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ABSTRACT 



We review recent progress in modehng the quark-gluon content of mesons and their low-energy 
interactions through the Global Color Model field theory. An emphasis is placed on techniques that are 
shared with the approach based on truncations of the Dyson-Schwinger equations of QCD. In contrast 
^ . to most other field theory models for QCD degrees of freedom in hadron physics, this approach directly 
^ ] deals with the derived intrinsic space-time extent of the meson modes in their role as field variables and 
. ^ ■ can accommodate confinement as well as dynamical breaking of chiral symmetry. Various theoretical 
^ techniques and approximations found useful in this approach are described. Selected applications 
'reviewed here include the properties and interactions of the Goldstone bosons, interaction vertex 
functions, low energy chiral observables, electromagnetic interactions and form factors, and transition 
form factors. Some initial considerations of Vector Meson Dominance and pion loop processes are 
discussed. 
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1 Introduction 



To investigate hadron physics from the perspective of nonperturbative quantum chromodynamics 
(QCD) and QCD-based models, a number of approaches are currently employed. These include the 
intensive computational method of lattice gauge field theory, the method of light cone quantization, 
QCD sum rules and chiral perturbation theory (ChPT). In one sense, the first two constitute a 
completely fresh start in that much of the empirical information on hadron dynamics that has been 
built up over many of the past decades through the refinement of model hadronic lagrangians is not 
easily joined. This becomes irrelevant when S-matrix elements for many-body hadronic processes can 
be routinely obtained from the relativistic gauge field theory of quarks and gluons. However it is 
apparent that for some time to come, model hadronic theories remain an important intermediate form 
into which new findings about non-perturbative QCD can be packaged for rapid deployment in hadron 
physics. One way to link QCD with such an intermediate form is to subject the quark-gluon action 
to a restructuring or change of field variables to expose hadronic or collective degrees of freedom and 
integrate out the quarks and gluons. At present, most work of this type that is able to confront a 
variety of hadronic observables deals with QCD-based model field theories. We review meson physics 
obtained from the Global Color Model (GCM) of QCD (Cahill and Roberts, 1985; Praschifka et al, 
1987a). Other overviews of this model that emphasize different aspects can be found in Cahill (1992) 
and in Cahill and Gunner (1997). 

What we refer to as a model hadronic field theory was in the past usually called an effective field 
theory. In either case this is in the old-fashioned sense that implies a semi-phenomenological model 
field theory meant for a limited low energy domain; but not implying an equivalent representation of 
an underlying theory through a systematic expansion and a counting scheme to order corrections. To 
distinguish the former from ChPT and other recent developments in hadron physics that are referred 
to as effective field theories in the technical sense, we shall use the terms model or approximate field 
theories. 

The GCM is a model field theory of quarks in which the interaction is between quark color currents 
mediated by an effective 2-point gluon function. The infrared form of the latter is the phenomenological 
element of the model, while it is constrained to match the known perturbative QCD behavior in 
the ultraviolet. The point coupling hmit gives the closely related and widely used Nambu-Jona- 
Lasinio (NJL) model (Nambu and Jona-Lasinio, 1961) that is simpler to work with. Special interest 
is taken here in the distinguishing features of the GCM that include a naturally accommodated 
confinement mechanism and the finite size of the qq meson modes. The derived interactions between 
these modes produce vertex form factors that are accountable to underlying quark-gluon dynamics. It 
is this aspect that may provide an interesting perspective to the thrust in nuclear theory that centers 
around hadronic field theory models. A key question is whether a consistent quark-gluon basis can be 
developed for important hadron interaction vertices that can support the present empirical description 
(especially of the form factors) thus correlating many hadronic parameters with far fewer sub-hadronic 
parameters. 

The Dyson- Schwinger equation (DSE) approach to non-perturbative QCD (Roberts and Williams, 
1994) is to truncate the infinite system of coupled equations that link the n-point functions of QCD at 
low order, use an Ansatz to replace the discarded information, and to seek a self-consistent solution. 
There is a close correspondance in form between the content of the GCM action at tree-level in 
the composite hadron fields and a particular low-order implementation of the DSE approach. For 
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this reason, a number of present calculations in each approach share common elements. A clearer 
distinction would become evident at the next level of treatment where the organization of the physics 
in each approach is different. However, they remain complementary; each providing insights into the 
role of nonperturbative QCD effects in hadron physics. We shall refer extensively to the mixing of 
elements of these approaches. 

In Sec. 1^ we introduce the Global Color Model and review the auxiliary field method for bosonization 
to an action for the meson modes of the theory. Also considered there is the decomposition to local field 
variables and the identification of the corresponding meson-quark vertex functions. That section ends 
with a discussion of issues related to extraction of physics from the formalism and calculations defined 
in Euclidean metric, and the relationship of the GCM and DSE approaches. Important elements for 
producing hadron observables from the GCM are the dressed quark propagators and forms used to 
facilitate calculations are summarized in Sec. ^ We examine in Sec. ^ the GCM description of the 
Goldstone boson sector, selected low energy chiral observables and methods used for other mesons. 
Electromagnetic coupling is developed in Sec. |^ where recent calculations of the pion and kaon charge 
form factors are reviewed. Meson interaction and transition form factors are discussed in Sec. |^. In 
Sec. selected issues beyond meson tree level are considered and these include vector meson dominance 
processes and pion loop contributions. Concluding remarks are presented in Sec. ||. 

2 Global Color Model 

The action for the GCM is defined in Euclidean metric as (Cahill and Roberts, 1985; Praschifka et 
al, 1987a) 

SccMiq, q]= J d'^x g(x)(7 ■ + m)q{x) + d^xd^y j^{x) g'^D^^{x - y) j^{y), (2.1) 

where m is a diagonal matrix of bare quark masses, and j^{x) = (i{x)^'~^uq{x) is the quark colored 
current. The employed Euclidean metric is such that a ■ b = a^Jj^, spacelike vectors satisfy > 0, 
and {7^,7i.} =25^^, with 7^ = 7^^. A model field theory of this form can summarize a number of 
diverse investigations in the modeling of hadron physics such as Cahill and Roberts (1985), McKay 
and Munczek (1985), Barducci et al. (1988), and McKay et al. (1989). This model may be viewed as 
a truncation of QCD and some consideration of that aspect can indicate the content and role that 
should be attributed to the phenomenological gluon two point function -D^jy. The connections between 
QCD and Eq. ( |2.1j ) have been considered by Cahill and Roberts (1985) and by McKay et al. (1989) 
and a brief summary follows. Functional integral methods are convenient for that purpose. In QCD 
(see, e.g. Marciano and Pagels, 1978) the generating functional for quark propagators is 

Z[f], ri] = N J DqDqDA exp (^-S[q, q,A\ + J (fx (r/g + gr/)^ , (2.2) 

where 

S% q,A]=J d'x [q{^ + m - ig^r)q + ^^^.^j , (2-3) 

and F^^ = d^A"^ — dyA'^^ + gf^'^'^A^^A^. We leave the gauge fixing term, the ghost field term and its 
integration measure to be understood. If we introduce 

exp iW[J]) = JdA exp (- J d'x {\f;,F;, - J^Al)) , (2.4) 
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the generating functional becomes 



Z[r], Tj] = N j DqDq exp (— + m)q + r]q + qi]) exp ( W 



(2.5) 



where the abbreviated notation imphes the usual spacetime integration. The functional W[J] is the 
generator of connected gluon n-point functions without quark-loop contributions. It may be written 
as 

W[J] = -j d'xd'y j;{x)D,,{x - y)J:{y) + Wn[J] (2.6) 
allowing the formal factorization 



Z[r],ri\ = exp Wr 



' 6 6' 



'GCM 



[ri,rj\. 



(2.7) 



Here the quark field variables in the argument of Wr have been replaced by the corresponding source 
derivatives. The truncated model identified this way is generated from 



ZccMirj, 1]]= N J DqDq exp(-5'GCA/[g, q] + VQ + QV), 



(2.8) 



where SgcmIqq], is given by Eq. (|2.1j ) in which the dependence upon the coupling constant g has been 
made explicit. The GCM corresponds to the retention of just the first term of Eq. ( ^.61 ), involving 
the gluon 2-point function excluding qq vacuum polarization contributions. The neglected term Wr 
involves the n-point functions, n > 3, from the pure gauge sector. Whether this can successfully 
model soft hadron physics can only be judged by consideration of the results that follow from a 
consistent phenomenological form for g^D^^[x — y). Certainly, any element of hadron physics that 
absolutely requires the gluonic coupling of three or more quark currents can be expected to resist 
description via the GCM; that finding in itself would constitute a valuable contribution. On the 
other hand, if a significant body of information can be correlated through a phenomenological gluon 
2-point function via accessible continuum calculations, this might highlight those aspects that deserve 
more concentrated attention through lattice-QCD methods (Rothe, 1992). Even though the model in 
Eq. ( |2.1|) looks Abelian in form (QED may be recast that way), the non-Abelian nature of the gauge 
field sector has not been totally ignored. Those non-Abelian aspects that contribute to the dressed 
pure-glue 2-point function can be implemented through the phenomenological element of the GCM. 



2.1 Bosonization to Mesons 

The auxiliary field method consists of a transformation of a field theory that introduces additional 
field variables through which important dynamical features can be revealed more readily than with the 
original fields. The method has found wide use in statistical mechanics as the Hubbard-Stratonovich 
transformation (Stratonovich, 1957; Hubbard, 1959). Early applications of the auxiliary field method 
to quark dynamics (Eguchi and Sugawara, 1974; Kikkawa, 1976; Kleinert, 1976a) were applied to 
four-fermion interaction models of the NJL type. Here the auxiliary fields are introduced so as to 
assume the role of hadronic collective modes generated by combinations of q{x) and q{x) at the same 
space-time point. The four-fermion interaction is exactly transformed into boson mass terms plus 
a local Yukawa coupling of boson fields to the quarks. The locality of the auxiliary boson fields is 
a consequence of the contact nature of the original four-fermion interaction. Since there remains 
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only a quadratic dependence upon the quark fields, they could then be integrated out leaving a 
bosonized action. Expansion of the action in quantum fluctuations about the classical values of these 
composite flelds was used to produce an effective action from which a number of important issues, 
e.g. renormalizabihty of four-fermion couphng (Eguchi, 1976) could be addressed. 

It was subsequently observed (Kleinert, 1976b; Shrauner, 1977) that the auxiliary field method can 
be generalized to account for the intrinsic nonlocality in the quark couplings generated by models 
that have an Abelian vector gauge field coupled to quarks. The vector gauge field can be integrated 
out to produce a finite range four-fermion interaction. The auxiliary boson fields needed to reduce 
this interaction to quadratic form are bilocal combinations of q{x) and q{y) whose free-field solutions 
turn out to be the ladder Bethe-Salpeter (BS) bound states. Expansion of the bilocal field in terms 
of these BS free-field eigenmodes produces (Kleinert, 1976b; Kugo, 1978) an effective action for the 
finite-size qq mesons of the model. 

In these works color was omitted and baryons were not considered. The achieved reformulations in 
terms of qq objects constitute a true hadronization for the case of a single Abelian vector gauge field 
coupled to any number of flavors of quarks. Such a field theory shares many properties in common 
with the strong interaction such as the conserved vector currents, global SU{Nf) (g) SU{Nf) chiral 
symmetry, and the consequent PCAC. Nevertheless, to obtain insight into approximate methods for 
QCD in the nonperturbative domain of hadron physics, a study was begun (Cahill and Roberts, 1985) 
on the bilocal field hadronization method for colored quarks that interact through finite range coupling 
of their colored vector currents. 

We shall consider only the Feynman-like gauge where g'^D^y[x — y) = 5fivD[x — y) and an Ansatz 
for D defines the model. This is done to permit the simplest form of the reordering identities that 
are to follow. It should not be taken to imply that the intent is to match with information from a 
corresponding gauge choice in QCD. Such a choice is not possible as the above Ansatz will depart 
from the perturbative form in both transverse and longitudinal parts while the latter will have no 
interaction or dressing in QCD. A reformulation of Sqcm in terms of fields that represent qq objects 
is achieved by the Fierz reordering identities 

l^jl'ki = , = (1, ^75, ^7., ^7.75) (2.9) 

for spin. 

It j 1 A*^ 

5ij5ki = -Fii-Ffcj > ^Nf=2 = 1 ^Nf=3 = (-^^ (2-10) 

for fiavor SU{2) singlet plus triplet (1/ -|- 3/), or SU{Z) singlet plus octet (1/ -|- 8/), and 

4 i 

KjKi = ^ii^kj ' = (3I' (2-11) 

for color SU (3) singlet plus octet (lc+8c). With the direct product set of matrices A^ = ^K°- C"^ 
where 9 — {a,b,c), the complete reordering needed to express the the quartic quark current-current 
coupling in terms of a quadratic coupling of qq objects is 
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and we use the summation convention for repeated indices. Hence in Eq. ( ^.11 ), this reordering causes 
q{x) to form a scalar product with q{y) and q{x) with q{y). An equivalent form that is useful later is 

^7^1^ M ^7^1^ = tr (A'M) (2.13) 

where M is any matrix. The Dirac matrix cr^j, does not appear in Eq. (|2.9| ) due to the use of the 
Feynman-like gauge. If one chooses a family of gauges parameterized by ^, a more complicated Fierz 
reordering involving a^u and ^-dependent weights of the spin matrices is obtained. Because of the 
phenomenological nature of the effective gluon propagator here, the simplifications of the Feynman-like 
gauge are appealing. 



With Eq. ( PJ^ ) the GCM action in Eq. can be written as 

S[q, q] = J d^x g(x)(7 ■ + m)q{x) d^xd^y j\x, yyD{x - y)j\x, y), (2.14) 

where the negative sign arises from the anticommuting property of the Grassmann variables q and 
q. There appear now bilocal currents J^{x,y) = q{x)A^q{y) that are hermitian, i.e. [J'^{x,y)]* = 
J'^{y,x), and carry the Lorentz and flavor indices of the candidate qq states of this model. However 
the color structure of the qq combinations consists of the singlet Ic as well as the octet 8c. The latter 
may be eliminated in favor of diquark correlations as discussed later. 



As a specific illustration of the content of the interaction in Eq. ( ^.141 ), consider, for Nj = 2, the point 



coupling limit where D{x — y) G5{x — y). The A in the color singlet channel are 

A^ = ^r" {l^,^75, ^7., ^7.75} , (2.15) 

where a = 0, 1,2,3 and = 1. It is convenient to separate out numerical factors ng and deal with 
the canonical set of matrices Kg defined by 

A^ 

i^e = — = r" {Id, «75, «7v, iluls} ■ (2.16) 
ng 



Then the interaction term from Eq. (|2.14| ) in the color singlet sector becomes 



= / d'x {Gs [{qr^q)' + (g*75r"g)'] + Gy [{qtl.r^q)' + (g^7.75r°g)^] } . (2.17) 

Here Gs = fG, and Gy = |G. This is one of the simplest forms of the Nambu-Jona-Lasinio model 
(NJL) (Nambu and Jona-Lasinio, 1961) as applied to quarks in the color singlet channel. For recent 
reviews of the NJL model as applied to hadron physics, see Vogl and Weise (1991) and Klevansky 
(1992). The coupling constant relation Gs = 2Gv evident here has its origin in the simple coupling of 
quark color currents. Applications are usually made from the point of view of an effective model that 
implements chiral symmetry and its dynamical breaking; the constants Gs and Gy are allowed to be 
independent and are fitted to data. The GCM retains the distributed coupling of the qq correlations 
via an effective gluon 2-point function D{x — y) and also maintains the relative Fierz weighting that 
accompanies that function to each meson channel. For large momentum g^, the magnitude of -D(g^) 
in the GCM is fixed by the requirement that perturbative QCD should be obeyed; this is not possible 
in the NJL since the constants Gs and Gy force the same coupling strength in the UV as in the IR. 



8 



The role asked of phenomenology is quite different in the two models and the outcome can also be 
quite different. 



We return now to Eq. ( |2.14D and before restructuring it, we note that the interaction term in Eq. 
is explicitly chirally symmetric as is its original form. However that is not true for a fixed value of 
the summation index 6. If one is to truncate that sum, it is necessary to retain a subset terms 
that mix under the chiral transformation to preserve the symmetry. Similarly for any symmetry e.g. 
an electromagnetic gauge transformation. The four-fermion term of the GCM action in Eq. ( p. 14 ), 



when viewed as the product of two bilocal currents, may now be recast through integration over 
auxiliary bilocal field variables that couple in a linear fashion to those bilocal currents. The functional 
path integral formulation is convenient for illustrating this. The generating functional for the quark 
propagators or n-point functions can be written as 

Z[f],ri] = N I DqDqexpS^-{q{-fd + m)q) + ^{J*DJ) + {m + QV)] , (2-18) 

where the bracket notation denotes the usual spacetime integrations and summation over discrete 
indices as is appropriate to the context. For example 

{J*DJ) = I d^xd^y J\y, x)D{x - y)j\x, y). (2.19) 

Since the normalization factor of the generating functional is irrelevant for the quantities of interest 
derivable from it, we consider now the restructuring that can be produced by multiplication of the 
generating functional in Eq. ( |2.18| ) by the constant 

[det{D-^)]--2 = I DBexp |-/ d'xd'y ^^^|^^^} = / ^xp . (2.20) 

Here one chooses auxiliary field variables B^{x,y) that have the same transformation properties as 
the associated bilocal currents J'^{y,x). In particular, if the flavor basis is hermitian, these fields 
are hermitian, that is, [B^{x,y)]* = B^{y,x). The integration measure is / DB = J DB'^DB^ ■ ■ ■. The 
result Eq. ( |2.2CI| ) is a direct generalization of the gaussian integral over a real, local field. After 
Eq. ( p.l^ ) is multiplied by Eq. ( p.20|) , the shift of integration variables B^{x,y) B^{x,y) + D{x — 



y)J^^{y,x), will eliminate the current-current (four-fermion) term \{J*DJ') in favor of a Yukawa 
coupling term {JB) = (qA^B^q) and the original quadratic term in B. 

The generating functional is now given by 

Z[f],r]] = N' J DqDqDBexp {qG-' [B]q) - ^{B*D-^B) + {m + gr/)) , (2.21) 

where the object similar to an inverse quark propagator that occurs here is 

G-\x, y) = (7 . 9, + m) 5{x -y)+ A'B\x, y). (2.22) 
We thus have Z[fj, rf\= N' j DqDqDB exp {—S[q, q, B] + {fjq + qrf)) where the new action is 

S%q,B] = J d'xdS |g(x) G-\x,y) q{y) + ^^^^^^|^}- (2-23) 
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The introduction of auxiliary boson fields thus allows the quark fields to be integrated out since they 
appear now only in a quadratic form. This integration is 

J DqDq exp + (m + qr])) = exp (TiLnG'^[B] + {f]Gri)) , (2.24) 

where the fermion determinant has been expressed as Det(G'^^) = exp TrLnG^^. This yields 

Z[f], r]] = N J DB exp {-S[B] + {f]Gr])) , (2.25) 

where the completely bosonized action is 

S[B] = -TiLnG-'lB] + / d'xd^y ^'(^' ^)^'(^^ ^) . (2.26) 

J 2D{x — y) 

This result is an exact functional change of field variables from quarks to qq bosons that preserves the 
physical content of the original action. 

The point coupling or NJL limit of this action is more familiar and can help illustrate the content of 
the more general distributed case. With D{x — y) — > G6{x — y), the necessary auxiliary fields take 
the special form 

B\x,y)^^Six-y), (2.27) 
He 

where ng are the purely numerical factors introduced by the Fierz reordering such that = ngKg 
(no summation), where Ke are Dirac matrix covariants (see Eq. ( |2.16| )). Thus for Nf = 2, Eq. ( p.23|) 
reduces to 

S^-^^[q,q,(j)] = J d^x i^q {-f ■ d + m + a + i-f^r ■ tt + i-f ■ oo + ■ ■ ■) q 



where Gg = UgG. After integration over the quark fields, the point coupling limit of Eq. ( p.26|) is 

2 2 

5^'^^[a,7r,o;,---] = -TrLn(^ + m + a + z75f-7f + ■■■)+ / d^x + ■ ■ ■ . (2.29) 

J 2Gs 

This is the familar bosonized form of the NJL model. 

With either the GCM or its point coupling limit, it is necessary to identify the classical ground state 
fields. From Eq. ( p.26| ) which covers both cases, the classical field configurations Bq{x, y) are defined by 
the saddle point condition 6S/6BQ{x,y) = and translation invariant solutions represent the classical 
vacuum about which quantum fluctuations may be developed. After functional differentiation of 
Eq. (|2.26|) this equation of motion is 

B^^{x -y) = D{x - y) tr [A^S{x - y)] , (2.30) 

where the associated propagator S depends self-consistently on Bq. In particular, 

S-\x -y) = {^-d, + m) 5{x - y) + S(x - y), (2.31) 



10 



where — ?/) = A.^BQ{x — y) represents the quark self-energy to the present level of treatment. From 
Eq. ( |2.30|) this self-energy satisfies 

S(a; -y) = D{x - y)A' tr [a'S{x - y)] = ^D{x - yh,S{x - yh^. (2.32) 

The second form is obtained by reversing the Fierz reordering using Eq. ( p.l3| ) and carrying out the 
color sum to produce = |- Eq- (|2.32|) is the Dyson-Schwinger equation (DSE) with a bare vertex, 
i.e. the ladder or rainbow approximation. This illustrates one of the efficiencies of the bilocal field 
method. Treatment of the bilocal fields B^{x — y) at the classical level generates a quark self-energy 
that sums an infinite subset of quantum loop Feynman diagrams in terms of the original quark and 
gluon fields. 

,4 

The momentum representation defined by /(r) = / ■(2^^^'^^ f^l^ allows the inverse propagator to be 
written as S~^{q) = if^A{q'^)+ B{q^) + m. The rainbow DSE for the Dirac scalar and vector amplitudes 
of S obtained from Eq. (p.32|) are 



A( 'i\ 1I 2 8 /• A ^ A{q^)q-p 

and 

It is useful to note that the only non-zero classical auxiliary fields are those for which the index 
6 corresponds to Ic for color, If for fiavor; for Dirac spin, there is only the scalar (which gives 
B{p^)) and vector 7^ (which gives A{p'^) — 1). 

In the point coupling or NJL limit, D[k) = G, and the rainbow DSE gives A = 1 and BijP') M 
where 

M='^Gtp- . (2.35) 

3 J (27r)4g2 + (]vf + m)2 ^ ' 

We have indicated the cutoff that is obviously needed here. Note that the (point coupling) NJL model 
leads to a free dressed quark propagator with mass M(A) = M(A) -|- m. The cutoff dependence cannot 
be removed and A becomes a parameter of the non-renormalizable NJL model. In contrast, the finite 
range property of the quark DSE, either Eq. (|2.33| ) and Eq. ( |2.34|) in the rainbow approximation 



needed at this stage in the GCM, or the exact DSE discussed later for the DSE approach, allows 
renormalization after regularization so that dependence upon a regulator parameter A is removed 
in favor of the renormalization scale /z where boundary conditions are imposed. The quark self- 
energy is necessarily momentum dependent and constraints from perturbative QCD (pQCD) can be 
implemented. There are a number of limitations that follow from the point coupling NJL limit and 
we will return to this topic in Sec. p^.3| after the contrasting features of the GCM approach have been 
developed. 

The bosonized action in Eq. ( |2.26D may be expanded about the classical configuration in terms of 
the quantum fiuctuation variables B^{x,y) = B^{x,y) — Bq{x — y) to define S[B\ = S[B] — S[Bo] 
as the action for the quantum modes. With functional generalization of the expansion ln(l + x) = 
X]n(~3;)"'"^^/n, and noting the fact that terms of first-order in fiuctuations about an extremum cancel. 
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one obtains 



S 



B 



This form of action, truncated to selected low mass mesons, was used to begin investigations nucle- 
ons as mean field solitons (Cahill and Roberts, 1985), and to initiate estimates of meson couplings 
(Praschifka et al, 1987a; 1987b). Also studies and estimates of the possible properties of diquark 
correlations in this framework were started (Praschifka et al, 1989). 

Although only color singlet components are present in the classical field configurations, the fiuctuation 
fields & in the present arrangement have, in general, both color octet and singlet components. The 
presence of 8c as well as Ic qq fields in the meson bosonization is not satisfying for several reasons. 
Color octet bosons are not physically realizable states and therefore may only have a role as virtual 
constituents of hadrons. However, as the gluon exchange interaction is repulsive in that channel, 
such fields may not be efficient as effective dynamical variables. Furthermore in order to treat color 
singlet baryons, diquark correlations in color anti-triplets are needed. A new change of field variables 
that leads to mesons and baryons was subsequently developed (Cahill, Praschifka and Burden, 1989; 
Reinhardt, 1990). The method is based on the other possible Fierz reordering of the original current- 
current interaction. That is, in Eq. ( p.l|) , this reordering causes q{x) to form a scalar product with 
q{y) and q{x) with q{y). The color part of this Fierz identity is 

KjKi = CiiC'^j + ClkC\j , = ^1 , Clf, = ^'^tik , (2.37) 

where euk is the antisymmetric Levi-Civita tensor. The first term introduces Ic qq combinations and 
the second term introduces 3c qq diquark as well as the corresponding 3c qq combinations. This leads 
to a bosonization in terms of meson and diquark bilocal fields and has been used to obtain a set of 
covariant Faddeev equations for baryons truncated to a quark-diquark format (Cahill, 1989; Buck, 
Alkofer and Reinhardt, 1992). For the purpose of discussing meson physics, the Ic sector obtained 
from the simpler qq meson bosonization is sufficient. 



The quadratic term of the action Eq. ( p.36| ) defines the propagator and the corresponding free field 



solutions indicate the physical interpretation required of the bilocal meson fields. For this purpose, 
it is convenient to use a momentum representation in which total meson momentum P is the Fourier 
transform conjugate of R = [x + y)/2, and the relative q — q momentum is q. For the equal mass 
case (e.g. pions with m„ = m^), q is the conjugate of r = x — y. Then the hermitian property 
[B^{r] R)]* = B^{—r; R) is equivalent to [B^{q; P)]* = B^{q; —P). In the cases with unequal masses 
for q and g, such as the kaon, the choice of relative momentum variable will be discussed as the 
need arises. The choice of non-hermitian fiavor basis needed to form charge eigenstates leads to 
[B^{q; P)]* = B^{q; —P) where B is the antiparticle field corresponding to B. 



The quadratic term of Eq. ( |2.36| ) can be written as 

'52[i3] = ^E/^^^''(?';^)*^.^^(?''?;^)^'(9;^)> (2-38) 

where the inverse propagator is 

-Do'liq', q; P) = 5e'eD-\q' - g) + (27r)^5(g' - q) tr [A^'5(g+)A^5(g_)l , (2.39) 
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with tr denoting a trace over spin, flavor and color and g± = g ± P/2. Note that S{q) is a diagonal 
matrix in flavor space and that the flavor component of the can be used to project onto the 



appropriate propagators for each flavor. The first term of the inverse propagator in Eq. ( 2.39 ) is 
the momentum representation of 1/D{r), while the second term is a quark loop with the momentum 
integral not yet carried out. The equation of motion for free mesons is 6S2[B]/6B'>' = which yields 

E / -^e'eW^ Pn)T:{q; Pn) = 0. (2.40) 

Here Pn indicates that solutions (labeled by f „) are possible only for special (on-mass-shell) values of 
the total momentum. A more recognizable form is produced by multiplication in the operator sense 
by D giving 

E / - [A''S{q^)A'S{q^)] T'^{q; P„) = -T':{q'; P„), (2.41) 

which is the ladder BS equation separated into components via a trace technique. To make this clear, 
use of the definition r„(g;P) = A^r^(g;P), and the reversal of the Fierz reordering via Eq. ( |2.13D , 
converts Eq. (|2.41|) to the standard ladder BS form 

f„(g'; Pn) = -J ^,D{q' - q)^l.S{q+)tn{q; P^)S{q^)^^,. (2.42) 

The ladder BS amplitude is f„(g;P) and it is a matrix of functions in the direct product space of 
color, flavor and spin, corresponding to mesons with masses given by P^ = —m"^. The solution label 
n denotes the mass as well as the quantum numbers J'^'" , G etc. In a given meson channel, e.g. 
TT, p, uj etc, a number of Dirac spin matrices are involved through the various A^ that are coupled 
in Eq. ( |2.41[ ) to form r„(g; P) = A^r^(g; P). The r^(q'; P) are not all invariants since in general they 



also contribute factors of momentum that are contracted with the Dirac matrices of the A^ to form 
true covariants whose coefficients are invariant amplitudes. With the a;-meson for example, the set 
of true covariants for T,y{q; P) includes not only the canonical Dirac matrix 7j,, but also the covariant 
qu^D involving the Dirac unit matrix, and q^'j ■ q, etc. 



2.2 Effective Local Meson Fields 



The free bifocal meson field solutions (Eq. (|J0D) are on-mass-shell ladder BS bound states. An 



expansion of the general interacting bilocal meson field variables in terms of the set of eigenmodes of 
that ladder BSE kernel serves to define a set of effective local meson field variables. Essentially the 
continuous internal qq degree of freedom is decomposed into this basis. The local fields will couple to 
quarks in a well-defined distributed fashion described by vertex amplitudes that are determined by 
the ladder BSE kernel. This technique is an important step to both clarify the physical content of the 
bilocal meson fields and to provide a link to familiar local field models and methods. This seems to 
have been first exploited by Kleinert (1976b) and Kugo (1978). Here we will follow a variation on the 
eigenmode expansion method as formulated by Cahill (1989). 



It is helpful to introduce an abbreviated notation so that the inverse meson propagator Eq. ( |2.39| ) can 
be expressed as 

V-\P) = + g{P), (2.43) 
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i4 

in terms of quantities that are operators in the space with the completeness relation 1 = J2e I j^;^ 
\0,q){9,q\ and normalization {9',q'\9,q) = (27r)^ 6g'o6{q' — q). In this operator form, is diagonal 
in 6', while Q{P) is diagonal in q. Then, the extension of the ladder BS Eq. ( ^.411 ) into an eigenvalue 
problem defined for any P^, namely 

E / -§^.D{q' - q)Ge>e{q; P)t%; P) = (a„(P^) - l)f ^'(g'; P), (2.44) 

can be summarized as 

DG{P)\tn{P)) = {an{P^) - l)|f„(P)). (2.45) 

Here r^(g;P) = (6', g|r„(P)). We shall refer to a„ as the eigenvalues and they are dimensionless 
and physical on-mass-shell solutions correspond to an{Pn = ~^n) — 0- That relation will be seen 
to provide the free equation of motion for the effective local meson fields that will emerge. The 
corresponding anti-meson solution r(g; P) = [C^^V{—q] P)CY has the same eigenvalue. In the limit 
of exact flavor symmetry, F = F. Due to the above-mentioned properties of the bilocal fields, the 
conjugate solution can be expressed as F„(g; — P) = (F„(P)|6', g). 

The symmetry property Qe'e{q; P) = Gee'iq', — P) allows the conjugate eigenproblem to be written 

{t4P)\giP)D = (f„(P)|(a„(p2) - 1). (2.46) 
With the convenient preliminary normalization, 

{tm{P)\D''\Tn{P)) = Sn.,n, (2.47) 

expansion in eigenmodes gives 

DQiP) = Y: |f„(P))(a„(P') - l){t„iP)\D-\ (2.48) 

n 

Thus the mode expansion of the bilocal inverse propagator is 

1^~\P) = T.D~'\tn{P))a4P'){tn{P)\D-\ (2.49) 

n 

Also of use is the operator for the free propagator of qq correlations that carry meson bound state 
quantum numbers. In this formulation it is the inverse of Eq. ( p.43|) and this gives 

V{P) =D- D g{P) V{P). (2.50) 

The more familiar form for this object can be recovered from this abbreviated or operator notation as 
follows. Form {6' , q'\T>{P)\6 , q)A^ summed over discrete indices, and use Eq. ( p. 12 ) to reverse the 
Fierz reordering. Then the first term of Eq. ( p.50| ) yields 

Ve>e{q', q; P) = D{q' - g) ® + . . . , (2.51) 

which is the one-gluon exchange term. The obtained index structure of the kernel term produces the 
ladder summation. Thus it is the ladder qq scattering operator or T-matrix operator that plays the 
role of the propagator of the bilocal meson fields as we have chosen to define them. Other choices are 
possible but all should agree near the mass-shell poles. 
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The eigenmode expansion, which is easily seen to be 

= E '^"^^^1^^"^^^' , (2.52) 

can be useful for analysis of how the composite meson fields propagate from, and couple to, given 
sources. Some of the issues that emerge can be pointed out with reference to Fig. and Eq. ( |2.52| ). 
The dependence upon the total t-channel momentum P is carried by the denominator a„(P^) (which 




Figure 2.1: Exchange of a qq correlation mode or meson between sources 
which might be considered nucleons. 



can be associated with an effective local propagator represented by the parallel qq lines in Fig. p.l]) , and 



also by the meson-gg vertex amplitudes r„(g;P) (which are represented by the cross-hatched circles 
in the Figure). These vertex amplitudes enter an integration over the quark amplitudes of the sources 
to produce effective meson sources which are momentum-dependent and constitute form factors. For 
a given mode, the distribution of momentum dependence and normalization between numerator and 
denominator of the representation in Eq. ( |2.52| ) is quite an arbitrary choice. This can be changed by 
a field redefinition which will also change the interactions accordingly if they are consistently derived. 
Only at the pole is there no ambiguity. Away from a pole the question of what is the propagator of a qq 
correlation with the quantum numbers of a given meson bound state cannot be answered uniquely by 
any separation of numerator from denominator in Eq. (|2.52|) or by any separation of one eigenmode 



in the sum from another. The above eigenmode representation is probably not a feasible route to 
calculations of the qq T-matrix. We view Eq. ( p.52|) as simply a formal device to make contact with 



effective local meson actions and models. The point coupling counterpart of this T-matrix in the NJL 
model has been used to investigate meson exchange mechanisms for the nucleon-nucleon interaction 
(Celenza et al, 1996). As we will see shortly, in that limit there is usually only one mode per meson 
in Eq. (p.52|) and the numerator in Eq. (|2.52|) becomes a constant. Thus the issues raised above do 



not arise because the dynamics is simpler. 



Substitution of the mode expansion for V ^(P) from Eq. ( p.49| ) into the quadratic or free term 
Eq. (|2.38|) of the bilocal field action yields 



(2vr) 



(2.53) 
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where the 6„(P) are effective local field variables defined by the projection (f„(P)|D ^\B), that is 

UP) = J:J^s ^nil'; -P)D-\q' - q)B%q; P). (2.54) 



That is, the bilocal fields B^{q;P) have been expanded in terms of the complete set of internal BS 
eigenstates according to 

B%q;P) = Y.r:iq;P)b^iP). (2.55) 

n 

The label n characterizes the ladder Bethe-Salpeter spectrum obtained from Eq. ( |2.44| ). The original 
set of bilocal fields B^ were finite in number and labelled by the spin-fiavor-color index 6. The 
eigenmode expansion has converted the continuous degree of freedom labelled by q into an infinite 
discrete classification included in the label n. Also included in n are the labels for the various spin- 
flavor representations that diagonalize the ladder Bethe-Salpeter kernel DQ. 

The form achieved in Eq. ( p. 53 ) is often loosely called a localization procedure. However it does 
not ignore meson substructure but rather produces a dynamically equivalent formulation in terms of 
local field variables that experience nonlocal coupling through the relevant BS amplitude. According 
to Eq. ( |2.36| ) the tree-level interactions of the mesons are generated by a quark loop. If the loop 
momentum is q the vertex describing the coupling of the local meson field bn{P) to quarks is given by 
Eq. (|2:55|) as 

qA'B'q = Y.q{q+)^n{q-,P)hn{P)q{q-). (2.56) 

n 

Note that it is the label n that characterizes the meson, and not the label 6', which identifies the 
various matrices making up the basis for the BS amplitude 

Physical normalization at the mass-shell is achieved by writing a„(P^) =(P^ + m\{P'^))Z~^ or equiv- 
alently a„(P^) =(P^ + m^)Z~^(P^), where m„ = niniP^ = ~^n) ^^e physical masses. Identifi- 
cation of the field renormalization constant 

Z„ = Z^{P' = -ml) = K{-ml)]-\ (2.57) 

produces the physical fields bn{P) = hniP) j \f2^. Then with a dynamical mass function defined by 

p2 + m„(p2) = a,(p2)Z„, (2.58) 

we have 

1 r d^P 

= ^ E y ^ UP){P" + m„(p2))6„(P). (2.59) 

Note that the dynamical mass function ■mn{P'^), which appears because the meson mode is not ele- 
mentary, has a zero first derivative at the mass-shell point. That is, the residue at the propagator 
pole is unity. This is equivalent to choosing a proper normalization for the BS amplitude. To see this, 
we note due to the choice of preliminary normalization in Eq. (|2.47|) , the eigenvalue can be written as 
a„(p2) = 1+ (f„(P)|e;(P)|f„(P)), that is 

a„(p2) = 1 + tr I ^ h{q; -P)5(g+)f„(g; P)5(g_). (2.60) 
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Application of the derivative d/dPu, followed by multiplication of both sides by Zn produces at the 
mass-shell 

Atr / ^ - 

J (27r) 



where we have defined 



(2.61) 

Tn{q;K) = ^J:^{q-K). (2.62) 

This is the physical BS amplitude since the relation in Eq. ( p.61| ) satisfied by this quantity is the 
canonical BS normalization condition in the form appropriate to the BS kernel being independent of 
total momentum. (Itzykson and Zuber, 1980). 

The minimum information that needs to be carried forward from the above analysis is that the 
bilocal fields may be factorized into ladder BS amplitudes and corresponding local field variables as 
in Eq. ( |2.56|) . An explicit and simple form for the resulting meson action can be presented if we make 



the approximation that only the canonical Dirac matrix covariant is employed for each meson. That 
is, each meson BS amplitude is described by just one scalar function Fp^, Fy, ■ ■ -. We can leave the 
normalization of the internal amplitudes to be revisited later and, for just the lowest mass modes, 
the meson action for effective local field variables vf, p^, tu, ■ ■ ■ can then be written, by substituting 
Eq. (|235|) into Eq. dpH), as 

oo / \n 

5[7r, p, u;, ■ ■ ■] = Tr ^ [S {i'-f^r ■ ntps + i'jf.uj^tv + i'j^r ■ p^ty H )]" 



n=2 



n 



J ^ 2D{x-y) ^ ' 

This dominant covariant approximation is for illustrative purposes as it allows a simpler view of the 
essential physics. It has been the common approximation in most work within the GCM until recently. 

In general, a meson BS amplitude involves a number of distinct matrix covariants. For example, the 
physical pion BS amplitude has the form F^^ = Z75-E7r(g; P) + ^757 ■ PFT,{q] P) + ■ ■ ■ where Et,, F^, ■ ■ • 
are scalar functions. The canonical Dirac covariant, which does not require momenta for its con- 
struction (e.g., 75 for the pion), is usually the dominant one as is confirmed by studies of the BSE 
such as that by Jain and Munczek (1993). The consequences of limiting the number of covariants 
allowed has not been subjected to systematic investigation, especially for various observables that can 
be generated. Some initial results in this line by Burden et al. (1997) indicate that for the light vr and 
K mesons the sub-dominant pseudovector amplitude can contribute up to 15% to the mass and up 
to 35% to the weak decay constant. For the next heaviest mesons, the contributions are indicated to 
be about 5% and to decrease rapidly with mass. The more general case of many covariants is already 
covered by the BSE eigenmode analysis covered earlier in this Section. Thus with general flavor and 
spin matrix BS amplitudes F, physically normalized, the action for lowest mass modes can be written 

d^P 



5[vr,p,c.,---] = ^E/^C(P) A;^(P)0„(P) 



+Tr 5] (f. ■ 7? + F>^ + . p-^ + . . .)f , (2.64) 

n=3 ^ 

where the local fields are 0„ = vr*, uj^, p^, ■ ■ -. The generalized inverse propagators are 

Kl{P) = Sum \P' + mliP')] + n„„(P) , (2.65) 
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where the off-diagonal elements that exist give the relevant mixed self-energies. It should be em- 
phasized that the meson fields at this level are bare or tree-level fields in the sense that quantum 
dressing effects from the cubic and higher order couplings among meson fields have not been applied. 
Nevertheless there is significant dynamical content in these bare fields as evidenced by the ladder 
Bethe-Salpeter structure they carry. From Eq. one may extract the distributed vertex function 

for a meson vertex of interest and its momentum structure will depend upon the internal structure of 
the mesons as well as the gluon substructure within the dressed quark propagators. All interaction 
terms allowed by symmetries are present in Eq. ( ^.641 ) and as an illustration, the low order form that 
is generated from use of a derivative expansion, and has been investigated (Praschifka et ai, 1987a), 



is 



d X 



]^[[d^Tif + m^vr^] + ^[p^^ ■ pf,y + mjp ■ p] + ]^[uj^^ u^^ + mlu^] (2.66) 



+ ^^e^j,^^tr(7f ■ rd^TT ■ tOi^tt ■ rd^n ■ tS^tt • r ) + ■ • 



Due to the quark loop structure, the chiral anomaly terms are properly embedded in this approach 
and will occur with the tensor e^,^o-r and a factor of i from the Euclidean metric. 



2.3 NJL Limit of Local Meson Fields 

Local meson fields also arise in the bosonized form of the point coupling NJL model. To explore the 
different dynamical content at work within the GCM, it is useful to summarize the point coupling 
limit of the analysis in the previous Section. With D{k) — > G and B^{q;P) — > (j)^{P)/nQ, Eq. (p.38|) 
for (S2[B] the quadratic term of the action, becomes 

= ^ E / 7t4 /*(^) ^o'liP) ^\P) ■ (2.67) 



(2vr) 



Here, the point limit has identified 



in terms of the bifocal inverse propagator defined in Eq. (|2.39|) . We have indicated that this integral 
must be regulated. One finds that Iq,1{P) = Sg/gGg^ + Jg/g{P) where the standard NJL model integral 
is 

Je,g{P) = tr f ^ [Kg, S{q^) Kg S{q^)] . (2.69) 

In general Jgig{P) is a matrix that is diagonal except for coupling between the channels characterized 
by the Dirac matrices Id and 75 and also and 71^75. The off-diagonal elements are often ignored 
in applications of the NJL model. In this case, the connection to the previous eigenmode analysis 
for the distributed coupling case is most transparent. Since the point coupling limit of Eq. ( p.44| ) 
shows that the eigenvectors r^(g;P) become independent of g, and the diagonal assumption makes 
them independent of P, we have oc 5nfi- That is, the modes are characterized only by the label 9 
which distinguishes among the flavor and Dirac matrices. There is only one mass state for each type 
of meson, and only a single covariant (the canonical one, i.e. 275 for vr, i'^i, for u;, etc) appearing in 
the BS amplitude. 
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When off-diagonal aspects of Ig,l{P) are included, as they should be for generality, diagonalization is 
needed for physical fields. This can be summarized as the point coupling limit of the distributed case 
described in Sec. p.2| . The relevant limit of Eq. ( p. 44 ) is 



EG,,J,,,(P)f^(P) = K(P^) - l)fr(P). (2.70) 



Thus using the preliminary normalization ^ = 6nm instead of the point limit of Eq. ( p.47| ), 

one finds the expansion 

lo'liP) =T.Ge'f'n{P) "n(P^) Ki-P)Ge\ (2.71) 



where 

a„(P^) = 1 + tr (01^" (-^) t'^^' ^(^+) ^(^-)] rn(^) (2-72) 

is the point limit counterpart of Eq. ( p.60|) from the GCM. The quadratic term of the meson action 
from Eq. ( p. 67 ) therefore becomes 

= ^ E / ^ 6:(P)«„(P2)6„(P), (2.73) 



which is identical in form to Eq. ( |2.53|) for the distributed case. However the local fields here come 
from the expansion 

<t>\P) = Y.^:iP)UP), (2.74) 

n 

in terms of the small number of basis eigenvectors F^. There are no more than 2 coupled terms since 
at most 2 Dirac matrix covariants are coupled by Jg'g{P) given in Eq. ( |2.69| ). When such off-diagonal 
coupling is ignored, the fields bn{P) are just rescaled 0^(P) fields. In either case, the analysis of 
Eq. ( |2.73|) to physically rescale the fields to produce unit residue at the mass pole proceeds in the 



same way as previously discussed. We note that the inverse of /^^ satisfies /(P) = G — G J{P) I{P) 
in a matrix notation and this is the qq T-matrix in the NJL model. It is the point limit counterpart 
of Eq. (|2.50| ) for the finite range ladder T-matrix operator "D^P) that arises in the GCM or truncated 
DSE approaches. 

Aspects in the description of meson physics can now be used to identify several consequences of the 
point coupling NJL limit by comparison to the GCM/DSE approach that retains the finite range 
nature of the interaction. For a successful NJL phenomenology, the cutoff A is typically rather low: 
A ~ 1 GeV. The lack of an interaction at higher momentum suggests that the physics will not 
be adequate for form factors of physical hadron processes except at low momenta and for coupling 
constants. The point coupling limit produces meson qq Bethe-Salpeter amplitudes that are momentum 
constants in the NJL. The quark loop integrals for meson interaction vertices will require a regulator 
parameter which is not necessarily the same as that for the self-energy. In contrast to this, the finite 
range or momentum dependent structure of the meson BS amplitudes in the GCM/DSE approaches 
(and in QCD) provide natural regularization of many of the quark loop integrals for meson vertex 
functions. On a more fundamental level, the presence of a momentum cutoff in the various DSEs of 
the NJL model destroys the translation invariance required for the variable shifts necessary to derive 
the various Ward-Takahashi identities (Ward, 1950; Takahashi, 1957) such as the one relating the 
photon-quark vertex to the quark propagator. An extension of the NJL model that modifies the point 
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interaction by matching to a one-gluon-exchange tail for momenta beyond a typical liadronic scale 
(~ Aqcd) has recently been put forward to address these issues (Langfeld, Kettner and Reinhardt, 
1996). 

The lack of a confinement capability in the NJL means that quark loops for physical meson processes 
will generate a spurious qq production threshold at s = 4M^. For a typical constituent quark mass of 
about 300 MeV, mesons heavier than the kaon would receive a spurious width. This problem can be 
avoided in the GCM/DSE approaches by taking advantage of the indications from studies of both the 
gluon and quark DSEs. This is that the gluon 2-point function shows sufficient enhancement in the IR 
to suggest confinement of both gluons and quarks through momentum-dependent dressed self-energies 
that do not produce a physical mass-shell. An extension of the NJL to add a phenomenological 
confinement interaction has been made recently by Celenza et al. (1995). There the confinement 
interaction is used selectively to dress the meson-gg vertices but the quark self-energy remains a 
constant. 

2.4 From Euclidean Metric to Physics 

Some important considerations in the type of approach considered here are the following. A typical 
first step is to use the effective meson action of Eq. ( |2.64| ) to produce, at meson tree-level, the meson 
self-energy amplitudes as integrals involving the dressed quark propagator S{p) and the ladder BS 
amplitudes F. Equivalently, the inverse propagators for the meson modes are calculated. An important 
structure that arises is a quark loop with two F vertices and one independent external momentum. A 
subsequent step is to calculate an interaction term such as the pnn vertex function or form factor. Here 
the calculation requires a quark loop with three F vertices and two independent external momenta. 
The approach here requires Euclidean variables for the momenta and thus the above-mentioned loop 
integrals for space-like values of external hadronic momenta are accommodated naturally. It will be 
necessary however to produce physical quantities, such as masses and coupling constants, defined at 
time-like values of external hadronic momenta. This has usually been done by an analytic continuation 
in the external momentum variables before evaluation of the loop integrals. The momenta occurring 
in the quark propagator parts of loop integrands are linear combinations of the internal integration 
momentum and the external momenta. Since one or more of the components of external momenta 
must be continued to complex values, the quark propagators must be evaluated in a domain of the 
complex momentum plane. 

This procedure is well-defined for the loop integrals mentioned above in those works we review where 
the dressed quark propagator is an entire function in order to implement quark confinement. The 
quark propagators will not contribute spurious production cuts to the loop integrals. The vertex 
functions of the loop integrals, being ladder qq BS amplitudes, have no production cuts either. For 
meson tree-level quantities, the needed continuation in the external momenta to define mass-shell 
quantities is not hindered by singularities. Beyond tree-level, there is, for example, a vrvr production 
cut in the p self-energy that first occurs at the one-meson-loop level of treatment. These quantities are 
explicitly constructed by a pion loop and the production cuts are automatically generated by the pion 
loop integration of functions that themselves have well-defined continuations previously established. 

The form of the GCM field theory model, and indeed that of the more general phenomenological 
approach based on the DSEs of QCD (Roberts and Williams, 1994), is such that use of the Euclidean 
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metric is necessary to draw upon the existing body of experience and results from previous studies 
of the dressed quark and gluon propagators. Such results with applications to hadron physics, either 
from lattice-QCD or from continuum studies through the DSEs, are invariably obtained in Euclidean 
metric for practical reasons. With the Minkowski metric, non-perturbative treatment of DSEs is 
extremely difficult due to the singularity structure and the indefinite norm. For the Bethe-Salpeter 
equation with scalars, a direct four-dimensional approach in Minkowski metric has only recently been 
developed and tested (Kusaka and WiUiams, 1995). There is no known treatment for fermions in that 
approach. 

The transcription of non-perturbative equations (such as the DSEs) from one metric to another is not 
necessarily equivalent to analytic continuation of the solution amplitudes in the (external) momenta. 
In the former case, establishment of the connection between solutions in each metric requires explicit 
integration along the closed contour generated by the Wick rotation (Wick, 1954). This is difficult 
except in very special cases. The singularity structure of the solutions in one metric dictates whether 
they are related to the solutions of the equations that difi^er only in form due to transcription into 
another metric. In the case of approaches that use an entire function quark propagator, there is an 
essential singularity at infinity and the Wick rotation cannot be used to justify change of metric by 
simple transcription of the form of equations. Analytic continuation of just the solution amplitudes 
(and not the total form of the non-perturbative equations) from the space-like domain near = 
to the time-like domain near = — m^, where m is not too large, requires analyticity information 
of the solution amplitude only in that neighborhood. This is the technique used in the GCM/DSE 
approaches to put external hadronic momenta on the relevant mass shell. 

One can take the point of view that the field theory in use is defined in Euclidean metric. The 
consistency of this point of view and the generation of physical quantities by continuation of the 
Euclidean n-point functions of a field theory, and of non-perturbative QCD models in particular, has 
been discussed in some detail by Roberts and WiUiams (1994). We note that this is the approach within 
lattice-QCD (Rothe, 1992) where the continuation of Euchdean n-point functions is achieved by taking 
the limit of large Euclidean time intervals in correlation functions of hadronic interpolating currents. 
Through the Cauchy integral analysis, this is equivalent to a continuation of the n-point functions 
in the Euclidean hadronic momentum component P4 from the real domain to the imaginary point 
corresponding to the nearest mass-shell. In the GCM/DSE approaches, specific model truncations 
are imposed that imply a limited set of quark-gluon interactions and dressings, but otherwise the 
Euclidean aspects and the definition of mass-shell quantities is equivalent to that of lattice-QCD. 
In the former case, those n-point functions that are presently used to provide the phenomenological 
input, such as the gluon 2-point function, should be viewed as interim parameterizations of dynamical 
quantities that are, in principle, well-defined from the Euclidean version of QCD and have clear 
counterparts in the lattice formulation. 

Given that the GCM/DSE approach involves analytic continuation of n-point functions (say /) in 
the external hadronic Euclidean momenta (say P^), there is more than one way to approach this. 
Initial investigations simply fit /(P^) on P^ > with known functions and then evaluated the fit at 
the mass-shell point where P^ < 0. Most of the work in recent years that we shall review here use 
the procedure that hadronic external momenta are to be continued before internal loop integrations 
are performed. As discussed at the beginning of this Section, this requires knowledge of various 
propagators and vertices in domains of the complex plane and the entire function parameterization of 
the confining dressed quark propagator has facilitated such work. A third method has recently been 



21 



put forward and tested by Burkardt, Frank and Mitchell (1997). Drawing upon the above-mentioned 
connections with lattice-QCD methods for extracting physical limits from the Euclidean metric, this 
method is to Fourier transform the external hadronic Euchdean energies P4 to Euchdean time and 
then take the large time interval limit. This removes the explicit need to determine or postulate the 
behavior of loop propagators and vertices in a complex momentum domain and, like lattice methods, 
requires only that the hadronic spectrum in a Lehmann representation of external hadronic correlation 
functions be known. In a simplified version of the impulse approximation or quark loop mechanism 
for the pion charge form factor, Burkardt, Frank and Mitchell (1997) have demonstrated that this 
method reproduces the results from the direct or second approach mentioned above. 



2.5 The DSE Approach 

Instead of considering the underpinnings of the present investigations to be a result of bosonization 
of the GCM action to produce an effective meson action, there is a more general viewpoint (Roberts 
and Williams, 1994). Selected truncation of the tower of coupled Dyson-Schwinger equations of QCD, 
together with use of a generalized impulse approximation, produces a formulation of meson vertex 
functions that are identical in form to what is generated in the effective meson action from the GCM. 
The content can be more general. In the generalized impulse approximation for form factors, the 
dressed quark propagator and the meson BS amplitudes are formally defined exact n-point functions. 
The interpretation of the role of phenomenology for those functions is then quite different from what it 
is for the same calculation taken from the tree-level GCM effective meson action where the dynamical 
content of those functions is ladder/rainbow DSE. Quantum meson loop dressing from the GCM would 
have to be applied to begin to recover a comparable interpretation. 

The distinction only becomes clearer when higher order terms are considered in each formulation. 
Progression through the loop expansion in the effective meson action derived from the GCM defines a 
specific ordering of quantum effects. The tree-level ladder mesons become dressed with ladder mesons, 
and then dressed again, etc. The point of view implemented through the DSE approach is that the 
quark-gluon content of important n-point functions for hadron physics is decomposed only in the 
way necessary for construction from other n-point functions that are active at the level of truncation 
being employed. There is potentially more freedom to develop an efficient ordering with the DSE 
approach as one moves beyond the impulse approximation for hadronic processes and considers less 
severe truncations of the DSE tower. Either route is difficult. Present efforts, for practical reasons, are 
focussed upon meson tree-level within the GCM, and meson and photon insertions in dressed quark 
loops within the DSE approach. At this level, the distinction between the GCM and the DSE is a 
subtle one of interpretation of dynamical content represented by the employed parameterization of the 
dressed quark propagators; this becomes important mainly as next order corrections are considered. 

We have mainly taken the effective meson action viewpoint that the GCM naturally leads to. One 
advantage of this is that it emphasizes the relevance to model hadronic field theories which summarize 
efficiently most of what is known and unknown about low energy hadron physics. It is becoming quite 
feasible to generate such models with many of the previously phenomenological coupling form factors 
now given a quark basis that, although approximate, captures the dominant influence on dynamics at 
the hadron size scale. For a quark-gluon model field theory to provide information on form factors for 
use in model hadronic theories, it is necessary to make contact with the same prescriptions traditionally 
used in such models for off-shell propagation of hadrons. This is really only defined by the form of 
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interacting hadronic lagrangian or action, and this is the reason for the emphasis upon the bosonization 
procedure described earher in this Section. At the very minimum, it provides a consistent identification 
of the propagators for the qq correlations with meson quantum numbers. The same information is 
available in the DSE approach, but consistent DSE truncations that amount to specification of a 
model action for mesons have not been developed so far. 



3 Quark and Gluon Dressed Propagators 

In QCD the fully-dressed and renormalized quark propagator is given by the Dyson-Schwinger equation 
(DSE) 

S-\p) = Z2[ii ■ V + mo(A)] + ^1 ^ / ^j-y^^-'^P - kh,S{k)Tl{k,p), (3.1) 

where mo (A) is the bare mass parameter and A is the regularization parameter which can be taken 
to be a cut-off for the integral which otherwise is divergent. Here the dressed gluon propagator 
D^y{q) and the dressed quark-gluon vertex rj^(fc,p) are the renormalized quantities and they satisfy 
their own DSEs required higher order n-point functions as input. The coupling constant g is also 
the renormalized one. The meaning of the renormalization constants is illustrated by the following 
relations 

S-\p) = Z, S],\p) , D-J{q) = Z, D-J^{q) , g = g^ , = Z^ , (3.2) 

where Zi, Z2, Z^ are the renormalization constants for the quark-gluon vertex function, the quark 
field and the gluon field respectively. Further details can be found in Itzykson and Zuber (1980). In 
quantum electrodynamics (QED), the Ward identity entails that Z\ = Z^- In modeling QCD via the 
quark DSE, it is common to model g^D^iylq) in terms of a running coupling constant and such an 
Abelian approximation leads to the simplification Z\ = Z^- 

At the chosen renormalization point /i, the constraint 

S'^{j))\pi=^fi=%-i -p^m (3.3) 

allows identification of ^2 (A, \x) and the renormalized current quark mass m(A). In general the cut-off 
A is different from the renormalization point /i. For sufficiently large /i and A, renormalized quantities 
for fixed /i become independent of A as A — 00. It is common to choose /i=A since for large enough 
A one can take Z^ as unity. Discussions of these issues can be found in Fomin et al. (1983) and 
references therein. For our purposes here, it is sufficient to note that Eq. (^TTj) has been subjected 
to many detailed investigations under a variety of approximations and renormalization schemes. A 
review can be found in Roberts and Williams (1994). Later works that provide useful details of 
numerical solutions under a subtractive renormalization scheme are Hawes and Williams (1995) and 
Frank and Roberts (1996). 

We have already seen that bosonization of the GCM produces a tree-level meson action that involves 
the dressed quark propagator from the rainbow approximation {r^{k,p) — >■ 7j,) DSE. Although there 
is a large body of early work that employed a phenomenologically successful form for g'^D^y[q), most 
of the more recent hadron physics investigations within the GCM format that we shall review in 
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the following Sections are based on use of parameterized model forms for dressed quark propagators. 
There are several reasons. 

Firstly it is difficult to prepare numerical solutions of the quark DSE that are sufficiently accurate and 
well-defined in the quark momentum domain required for the subsequent multi-dimensional integrals 
needed to produce hadron physics quantities. One of the complications is that while the solution 
is more readily, and hence usually, addressed in Euclidean metric, continuation of the solutions to 
the complex quark momentum domain is usually required for application to hadron physics due to 
mass-shell constraints on external hadronic momenta. When improvements over initial analyticity 
assumptions were sought, complex plane singularities generated by the rainbow form of DSE proved 
problematical (Stainsby and Cahill, 1992; Maris and Holties, 1992). The realization that dressing 
of the gluon-quark vertex and a realistically strong infrared-enhanced 2-point gluon function could 
produce a viable implementation of quark confinement (Burden, Roberts and Williams, 1992; Roberts, 
Williams and Krein, 1992), and could effectively remove propagator singularities, led to a focus on 
this aspect. Secondly, there is at present insufficient knowledge of the infra-red behavior of the n-point 
functions of QCD to prefer numerical solutions of the quark DSE (or other DSEs) over reasonable 
phenomenological models of the infra-red behavior for the purpose of investigating whether a wide 
selection of hadron physics can be reproduced. In fact it is worthwhile to take the point of view 
that well-selected hadronic observables might provide a valuable constraint on the allowable infrared 
behavior thus providing a reference point to link with, and be confronted by, more direct methods 
such as lattice QCD. 

To discuss the principal models of the dressed gluon 2-point function and the quark propagator that 
have recently been used within the general GCM/DSE format, we first briefiy recall the leading 
asymptotic behavior of the gluon and quark propagators in QCD. A more detailed treatment of this 
topic and the relation to renormalization procedures can be found in the review by Roberts and 
Williams (1994). In Euclidean metric the Landau gauge gluon propagator is 

where n(g^) is the gluon vacuum polarisation. If ghost contributions are unimportant for soft hadron 
physics, the effective coupling constant in Eq. (|3.4| ) satisfies the same renormalization group equation 
as the QCD running coupling constant a{q^) (Bar-Gadda, 1980). This leads to the approximate form 

9 D^,{q) = 16^, - 1 —^^ . (3.5) 

for the dressed gluon propagator in this gauge. This "Abelian" form accompanied by the QED-like 
relation Zi = Z2, is a common approximation from which many studies of the quark DSE have begun. 
It is reasonable to begin investigations of modeling hadron physics from QCD with this Abelian-like 
approximation. Ghost contributions in DSE studies have been studied in Landau gauge and shown 
not to modify the qualitative features of quark and gluon propagators (Brown and Pennington, 1988a; 
1988b; 1989). Quantitatively, ghosts provide a small (< 10%) effect. 

In Landau gauge, the two-loop renormalization group expression for the running coupling constant 
only receives small (~ 10%) corrections from higher orders for spacelike-g^ > 1 GeV^ and hence 
can be said to provide an accurate representation on this domain (Brown and Pennington, 1989). For 
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(3.4) 



< 1 GeV^, however, a(g^) can only be calculated nonperturbatively and it is not known. The recent 
status of lattice-QCD studies of the transverse gluon 2-point function is represented by Bernard et 
al. (1994). The current status of DSE studies is summarised in Roberts and Williams (1994). Present 
phenomeno logical quark-DSE studies rely on an Ansatz for a(g^ < 1 GeV^) motivated by these gluon- 
DSE studies. For spacelike ^ ^qcdj ^^e running coupling constant is given at leading-log or 
one-loop order by 

where d = 12/ (33 — ^Nf) is the anomalous dimension of the mass, Nf is the number of quark flavors 
and Aqcd ~ 0.20 GeV is the scale parameter of QCD. 



The general form of the dressed quark propagator is 
S{p) = -i^ .pav{p^) + o-5(p2) = 



1 



^7 ■ p A{p'^) + B{p^) + m 



(3.7) 



When 7^ 0, there is dynamical chiral symmetry breaking (DCSB). The behavior of the above am- 
plitudes in the perturbative or UV asymptotic region is well known from the QCD renormalization 
group and operator product expansion (OPE) and QCD sum rules ( Politzer, 1976,1982; Gasser and 
Leutwyler, 1982; Reinders et al., 1985). The behavior has been summarized in the course of explicit 
numerical solutions and model building (Williams, Krein and Roberts, 1991). For q^ ^ ^qcd the 
leading-log result for the running mass (M = {B + m)/A = (Ts/(Jv) is 



M(g2 



m 



d ' 



(3.^ 



where m is a scale invariant parameter to be eliminated in favor of the renormalized current quark 
mass m =M{fi'^) that has been chosen for the constraint Eq. ( |3.3| ). This holds when m 7^ 0, i.e. for 
explicit chiral symmetry breaking (ECSB). When there is no ECSB (m = 0), there is exact chiral 
symmetry, and conservation of the axial-vector current leads to (in Landau gauge) 



M{q^ 



K 



.00 ^2 



where k is a constant independent of fj, and is given by 

An'^d (qq) 



l-d ' 



(3.9) 



ln(/iVA^^^) 



(3.10) 



This implies the scaling behavior (qq) ~ ln(yU^/AQ^£,) for the quark condensate; hence m (qq) is a 
renormalization scale invariant quantity. The quark condensate (qq) is a measure of dynamical chiral 
symmetry breaking; in the absence of ECSB, the condensate for each flavor is 



(?9)u = - li"^, ^T^sc So{x) = -12 



(27r)4pM2(p2)+52(p2 



(3.11) 



where So{x) is the chiral limit quark propagator and where the trace tr^c is over spin and color. Quoted 
values for quark condensates in this paper will be at the scale fj? = 1 GeV^. 
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Early investigations with the GCM (Cahill and Roberts, 1985) employed the simple confining model 
quark propagator due to Munczek and Nemirovsky (1983). This is defined by the use of the ^-function 



5^^{2'nY 5^{q) 3aVl6 in the rainbow DSEs, Eq. (I^SD and Eq. (|j|). Equivalently, 



form -D^j.(g) - . 

Landau gauge can be used and the same results are obtained by replacement of 6^^ above by |(5^i/ — 
q^Qv/ff)- This Ansatz for the dressed gluon propagator models the infrared behavior of the quark- 
quark interaction in QCD via an integrable infrared singularity. The chiral limit DSEs reduce to non- 
linear algebraic equations and the solution that exhibits DCSB, also minimizes the vacuum energy 
density. The solution is 



A(/) 



2, 
1 + 



2a: 

r)2 



(a^ - 
0, 



p'>f 



(3.12) 



For < ^ this gives the simple mass function M^(p^ 



a 
4 



— p^. There is no mass-shell, i.e. 



+ M'^{p^) 7^ for any p^, and there is quark confinement. In model investigations, the above 
behavior was approximately corrected for the infiuence of a perturbative UV tail for the gluon 2- 
point function. Since recent realistic studies of the DSE for the gluon 2-point function (Brown and 
Pennington, 1989) have found strong infrared enhancement in qualitative agreement with a regularized 
singularity, the above Ansatz is more representative for low p^ than it might seem. 



A number of detailed studies of mesons, diquarks and the nucleon within the GCM (Praschifka et ai, 
1989; Burden et ai, 1989; Cahill et ai, 1989a; Hollenberg et al, 1992) were carried out with quark 
solutions of the rainbow DSE generated from g'^D^u{q) = 5^yD{q) with 

= ^ exp (-,VA) + . (3.13) 

with A = 0.19 GeV, A = 0.002 GeV^ x = 1-14 GeV, and e = 2.0. The first term simulates the in- 
frared enhancement and confinement; and the second term matches to the leading log renormalization 
group result in the UV. The parameter e can be varied in the range 1 — 250 without significant change 
in results. This was found to be phenomenologically successful for selected light mesons ( tc,p/lj, /i), 
various diquark correlations and the nucleon. Hadron mass-shell information was obtained for < 
(time-like) from extrapolated > calculations. Efforts to improve on this experienced difficulties 
with the numerical continuation of the resulting quark propagator into the complex plane due to com- 
plex conjugate pairs of logarithmic branch cuts (Stainsby and Cahill, 1992) that often are generated 
by rainbow-DSEs. This considerably complicates the extraction of hadron observables beyond the 
lightest Goldstone bosons. 

Important guidance for extending the range of hadronic observables accessible through modeling the 
quark propagator in the GCM/DSE approaches is provided by the work of Burden, Roberts and 
Williams (1992). This provides an analytic solution to the DSE in Eq. ( p.l|) with the following simple 
one-parameter Ansatz 

g'D,,{k) = (5,, - ^) 87i'D6\k) and r,{p,p) = -id,S-\p) , (3.14) 
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where D is a mass-scale parameter. This Ansatz ensures confinement, in the sense described below. 
The dressed-quark gluon vertex used here, and resembling the Ward identity constraint (Ward, 1950) 
for an Abelian theory, is a result of extensive analysis of its general form (Ball and Chiu, (1980); 
Curtis and Pennington, (1992)). The explicit solution of Eq. (|3.1|) obtained by Burden, Roberts and 
Williams (1992) is 



as{x) 



-— exp(-2x)Ji(4m|/) + — / f/a^i(K) UvO^W ' (^.15) 

imy y Jo \ ° J 

with m = m/ V2-D, x = y"^ = j {2D), and where Ji and Ki are Bessel functions, and 

avix) = — las{x) + ^-^asixU , (3.16) 
m \ 4:1/ ay I 



with CTsix) = \/2D (Jsip^) and cry(x) = 2 D av{p^)- In Eq- ( P-15| ), the first term is finite in the 



chiral limit and it is associated with dynamical chiral symmetry breaking. The parameter C is not 
determined by Eq. ( |3.1|) with Eq. ( 3.14 ). The second term of Eq. ( p.l5| ) is zero in the chiral limit and 



is associated with explicit chiral symmetry breaking. 

An important feature of this model, and its derivatives described below, is that both as and ay are 
entire functions in the complex-p^ plane. As a consequence the quark propagator does not have a 
Lehmann representation and can be interpreted as describing a confined particle (Roberts, Williams 
and Krein, 1992). This becomes evident when quark loop contributions to hadron couplings are 
considered. The pion electromagnetic vertex considered later in Eq. ( ^.311) is one example. The 



entire function property ensures the absence of free-quark production thresholds, under the reasonable 
assumptions that the pion BS amplitude is regular for spacelike-p^ and the dressed photon-quark vertex 
is regular for spacelike photon momenta (Roberts, 1996). There are other, more complicated, ways to 
realize the absence of free-quark production thresholds as discussed by Roberts and Williams (1994) 
but the effect is the same. Some of the phenomenological implications of a model with a simple 
realization of this confinement mechanism have been discussed by Efimov and Ivanov (1993). 



For the quark loops needed for hadron physics the entire function model given by Eqs. ( p. 15 ) and (p. 16 ) 



has a well-defined behavior in the complex plane. However, the involvement of Bessel functions 
and the integral representation complicate applications. With the expectation that the details of that 
solution are less important than the qualitative features, Mitchell et al. (1994) and Mitchell and Tandy 
(1997) have employed a simplified version in studies oi — uj mixing and the pTrvr vertex. This was 
obtained by expansion of Eqs. ( |3.15| ) and ( |3.16D to first order in ffi. This produces 



<Js{x) = C(m) e-'" + ^ (l - , (3.17) 

_ 1 _|_ p-2x 

ay{x) = C(m)me-2^ (3.18) 

2x'^ 

where recognition is given to the fact that the previous undetermined parameter C should depend 
upon rh if the current mass dependence of realistic DSE solutions (Williams, Krein and Roberts, 1991) 
is to be simulated. This is also necessary to produce a reasonable value for m^. It was verified by 
direct comparison that, for typically small current masses, these approximate expressions are accurate 
representations of the original expressions to better than 0.1% in the complex domain needed for 
a typical quark loop such as the self-energy of the p. 
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Figure 3.1: Confined dressed quark amplitudes A (right scale) and B 
and the dynamical mass function M (left scale) corresponding to the 
simplified propagator form in Eqs. (|3.17|) and (|3.18| ). Taken from Mitchell 
and Tandy (1997). 



The parameters A = y2D = 0.889 GeV, C(m) = 0.581, m = 16 MeV, produce the soft chiral 
physics quantities U = 90.1 MeV, m„ = 143 MeV, and (gg)^2 = (-173 MeV)^ at /i^ = 1 GeV^ The 
formulation of calculations such as these will be described in Sec. [4.1| . With the produced pions and 
vector mesons from this model, the pion loop mechanism gives rup — = —21 MeV (Mitchell and 
Tandy, 1997) which is in reasonable accord with the experimental value —12.0 ±0.8 MeV. The dressed 
quark amplitudes ^(p^), -B(p^) as well as the dynamical mass function M(p^) = {B{p^) +m)/A{p'^) 
that follow from Eqs. ( |3.17| ) and ( ^.18| ) are displayed in Fig. O. This model was found to be too 
crude to accurately represent typical DSE solutions at > 1 GeV^. Apart from the weak strength 
at such momenta, and other details, the behavior shown in Fig. ^TT| is typical of the entire function 
propagator models discussed here. 

The large spacelike-p^ behavior of the original entire function model solution, given in Eqs. ( ^.15| ) and 
(P.16|) , is found to be 



m 



m 



+ ... 



and ay {p^ 



D 



m 



+ ... 



(3.19) 



p^ p* 

The leading behavior of oy in this model agrees with that from QCD, but the leading behavior of as 
in this model is missing the QCD logarithmic correction due to the anomalous mass dimension. In 



particular, since (Tsij>^) — ^ M{jP')/p'^ the leading QCD behavior deduced from Eq. ( |3.8|) is 



m 



(^s{p ) ~ — 

p2 



HpyAlao) 



(3.20) 
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Thus the model defined by Eqs. ( ^.11 ) and ( p. 141 ) incorporates asymptotic freedom neglecting In(p^) 
terms. 



The main defect of this model solution can be seen by setting m = in Eq. ( p. 151 ), which yields 
as{x) = Ce~^^. This is a poor representation for finite momentum since without a bare mass and 
when chiral symmetry is dynamically broken, Eq. (^.91) implies that 



K 



p4 



1-c 



(3.21) 



This defect results from the fact that although the form of Df^i, in Eq. ( [3.14 ) generates confinement, 
it underestimates the strength of the coupling in QCD away from k"^ = 0. In the numerical studies of 
Eq. ( p.l| ) that have used a better approximation to D^v{k) (Williams, Krein and Roberts, 1991) there 
is no such defect. 



Table 3.1: Parameters for the model quark propagators in Eq. ( |3.23| ) 
and Eq. ( p.24 ). Taken from Burden, Roberts and Thomson (1996). Note 
A = 10~^ is not varied. The scale parameter is D = 0.160 GeV^. At 
the finite bare mass values, = = 0. The resulting vr and K 





u/fi-quark 


s- quark 


m 


0.00897 


0.224 


bo 


0.131 


0.105 


hi 


2.90 


2.90 


b2 


0.603 


0.740 


h 


0.185 


0.185 




0.121 


1.69 



To expedite investigations of a variety of hadron physics quantities, a useful parameterization of 
quark DSE experience has been developed from modifications to Eqs. ( |3.15|) and (|3.16| ) that simply 
restore the missing strength at intermediate-p^. This provides a better approximation to the numerical 
solutions, while retaining the confining characteristics present in Eqs. ( p.l5| ) and ( |3.16| ). The following 
approximating algebraic forms were introduced by Roberts (1996) in the u/d sector and extended by 
Burden, Roberts and Thomson, (1996) to include the s sector. The scalar and vector parts of the 
quark propagators are defined in terms of dimensionless functions: 



(3.22) 



with s = p^, X = s/X^, = 2D, where D is a mass-scale parameter. The explicit forms are 



™^ x + Thj 



4{x) = 



b{ 



X 



hlx 



bl + U 




(3.23) 
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and 



a 



^) = ^ ihr-TV, . (3.24) 



Here = rrif/X, and A = 10~^ is not a free parameter. In this work the u and d quarks are 
considered to be identical, except for their electric charge. The dressed-quark propagator described 
by Eqs. ( |3.23| ) and ( |3.24D is also an entire function in the finite complex p^-plane and therefore retains 
the same confined particle interpretation as the previous analytic model solution that motivates it. 
The exponential function form that ensures this is suggested by that analytic solution of the model 
DSE given in Eq. ( p.l5| ) and Eq. (|3.16|) . The behavior of Eqs. ( p.23| ) and ( p.24|) on the spacelike-p^ axis 
is such that, neglecting quantitatively unimportant ln[p^] corrections associated with the anomalous 
mass dimension in QCD, asymptotic freedom is manifest. In Eq. ( p.23| ) the third term allows for the 



representation of dynamical chiral symmetry breaking, and the second term represents explicit chiral 
symmetry breaking. 



In Burden, Roberts and Thomson, (1996), the five parameters {m„, b^, . . . , b^} in Eqs. and (WM) 

were varied in order to determine whether this model form could provide a good description of the pion 
observables: fn] m^; (gg); r^r; the tt-tt scattering lengths; and the electromagnetic pion form factor. 
The procedures followed were essentially those developed in an earlier work (Roberts, 1996) where more 



details can be found. A very good fit was found with the w-quark parameter values listed in Table 3.1 



The low energy vr — vr scattering partial wave amplitudes were not part of the fitted data, but a 
reasonable description of the qualitative features of the data for / = 0, 1, 2 and E < Am.,^ was obtained. 
Dyson-Schwinger equation studies (Williams, Krein and Roberts, 1991) indicate that while it is a good 
approximation to represent the u- and d-quarks by the same propagator, this is not true for the s-quark. 
For example; contemporary theoretical studies suggest that 2ms/ {f^u + ^d) ~ 17 — 25 (Particle Data 
Group, 1996) and (ss) ~ 0.5 — 0.8 (uu) (Narison, 1989) which is a nonperturbative difference. In 
Burden, Roberts and Thomson, (1996), with this in mind, the model forms in Eqs. (|3.23|) and (|3.24|) 
were employed in a study of the kaon observables: fx] (ss); rj^o; r/^±; and the electromagnetic form 
factors of the charged and neutral kaon. The sensitivity of these observables to rhs and (ss) was 
too weak for an independent determination and therefore rhs = 2bfhu and 6q = 0.86g, which ensures 
(ss) = 0.8 (mm), were chosen for consistency with other theoretical estimates. The parameter b\ was 
allowed to vary to provide a minimal residual difference between the u/d- and s-quark propagators and 
a very good fit to the kaon observables was obtained with the values listed in Table |3.1|. The resulting 



TT and K charged form factors are displayed in Fig. 5.3, Fig. 5.4 and Fig. 5.5. The predictions for other 



TT and K observables are displayed in Table |3.2| where a\ are the vr — vr scattering lengths. The notation 
^ave _ 1^^^ -^md)/2 has been used and the scale is fj? = 1 GeV^. The quoted experimental values 
of "^iQgva) ""^iGeV^' ('^'")iGev2 and (ss)iqcv2 are representative of current phenomenology including 
QCD sum rule analysis. The calculation of Qt^o^.^ is from Roberts (1996), and F^'^{Arn^) is from Alkofer 
and Roberts (1996). Predictions for other observables from use of this model, and related models, are 
discussed in following Sections. 



It is desirable to confirm that the above quark propagator modeling is indeed representative of what 
can be generated from the quark DSE with an appropriate model gluon dressed propagator. An initial 
step in that direction was taken recently by Frank and Roberts (1996) in showing that a one-parameter 
model confined-gluon propagator used in the rainbow-DSE produces a quark propagator with no pole 
on the real axis. The qualitative behavior for real is similar to that of the above entire function 
parameterizations except the there is a broad resonance-like peak centered at ^ —0.55 GeV^. This 



30 



Table 3.2: tt and K observables calculated by Burden, Roberts and 
Thomson (1996) using the parameters of Table p.l| . See text for ref- 
erences for the various calculations. Experimental values are extracted 
from Particle Data Group (1996) except that r^r is from Amendolia et 
al. (1986a), r^^ is from Amendolia et al. (1986b), and r^o is from Molzon 
et al. (1978). The experimental vr-vr scattering lengths, Oj, are obtained 
from Pocanic (1995) and Gasser (1995). 
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0.87 


0.88 ± 0.06 



suggests an interpretation of a " constituent- quark mass". The employed gluon propagator model in 
Landau gauge g^D^^{q) = {5^^ - q^,qu/q^) A{q'^) is specified by 



A{q' 



47r2rf 



47i^m^,S\q) + 



,(-gV[4m?]) 

q2 



(3.25) 



where d = 12/(33 — 2Nf), with Nf = 3 the number of light flavors. The first term provides an 
integrable, infrared singularity, which generates long-range effects associated with confinement, and 
the second ensures that the propagator has the correct large spacelike-g^ behavior, up to In(g^)- 
corrections. The reason for the relationship between the coefficients of the two terms in Eq. ( |3.25| ) is 
that 



A fx' 



d^q 



'■"A(g2 



d 



(3.26) 



shows that the effects of S'^lq) are completely cancelled at small x"^ where A(a;^) — > d/x"^ + 0(a;^) 
which is the form expected from QCD without logarithmic corrections. Hence rrit can be interpreted 
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as the mass scale that marks the transition from the perturbative to the nonperturbative regime. The 
single parameter rrit was varied and = 0.69 GeV provided a best fit to a range of calculated pion 
observables. The pion ladder BSE was solved with this model confirming that the chiral hmit quark 
amplitude Bo{q'^) provides an excellent approximation for the 75 part of the BS amplitude for a 
typical finite mass. We note that the renormalization point used by Frank and Roberts (1996) is 
/I ~ 9.5 GeV which is verified to be in the purely perturbative domain. The cut-off used for the quark 
DSE integral is A ~ 5500//. These values are typical for solution of the quark DSE for momenta scales 
needed for hadronic physics apphcations. The observables were verified to be independent of A to 
about 3%. Complete independence from A would occur if the bare vertex were replaced by a dressed 
one that preserves multiplicative renormalizability. 



4 Light Mesons 

The spectroscopy of light-quark mesons is made interesting because of the role played by dynamical 
chiral symmetry breaking, the natural scale of which is commensurate with other scales in this sector. 
It also explores quark and gluon confinement because most vector and axial- vector mesons have masses 
that are more than twice as large as typical constituent- quark masses. An important feature of the 
strong interaction spectrum is the exceptionally low mass of the pion. Furthermore, the pion mass must 
vanish in the chiral limit; i.e., when the bare or current mass of the quarks vanishes even though this 
does not entail a vanishing of the constituent-quark mass. These observations are an indication that 
the characteristics of the pion as a manifestation of the Goldstone theorem is a particular and crucial 
feature of the strong-interaction spectrum. These features can be captured in the form of relativistic 
potential models (Bernard, Brockmann and Weise, 1985; Le Yaouanc et al., 1985) that incorporate 
the consequences of the axial Ward identity from field theory. Without such a parameter-independent 
property embedded in the format, the Goldstone theorem requirements arc not well simulated through 
specific parameter choices (Sommerer et al., 1994; Gross and Milana, 1994). As a model field theory, 
the GCM can easily respect the Goldstone theorem independent of phenomenological parameters. 



4.1 The Goldstone Boson Sector 

In the chiral limit, the pseudoscalar octet are massless Nambu-Goldstone bosons associated with the 
dynamical (or spontaneous) breaking of chiral symmetry. In this situation, the quark propagator 
fully specifies the BS amplitude for these pseudoscalar mesons. At the level of tree level mesons 
in the GCM, this Goldstone theorem is generated exactly by the consistency between the ladder 
BS equation for the pseudoscalar amplitudes and the rainbow SDK for the employed dressed quark 
propagator. In particular, the chiral limit ladder BS equation for the amplitude Tb{q', P) becomes 
identical to the ladder Dyson- Schwinger equation for -60(5^) = -B(g^;m = 0), the scalar part of the 
dynamically generated quark self-energy (Delbourgo and Scadron, 1979). Thus the chiral limit pion, 
in this approach, is automatically both a massless Goldstone boson and a finite size qq bound state. 

Thus the only knowledge of the underlying effective gluon propagator that is needed in this limit is 
that which is already encoded in the dressed quark propagator. This suggests that, as one moves away 
from this limit to generate physical masses for the PS octet through small finite bare quark masses, 
explicit knowledge of the underlying effective gluon propagator is a secondary consideration. This is 
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certainly true to first order perturbation in the bare masses. For this reason, a considerable effort 
has emerged in recent years to develop an efficient parameterization of the dressed u, d and s quark 
propagators constrained only by pion and kaon physics. This effort has been very successful and we 
will recount some of the major results here. 

We first demonstrate that, in the chiral limit, and at the ladder-rainbow level, each member of the 
PS octet is massless and has the 75 part of BS amplitude given by the quark scalar self-energy. 
After projection onto the canonical covariant A'' = 275A^-\/2/3, where we take SU{3)f, the ladder BS 
Eq. ( |2.41|) can be written as 

r,(g'; P) = -lj - [^lAS{q+)^l,\S{q^)\ T,{q; P), (4.1) 

where P^ = — m^. Throughout this work the symbol tr denotes a trace over SU{'i) color, Dirac spin, 
and the relevant flavor indices. For a flavor symmetric S{q), and using iif{X\Xb) = 2, we have at 

^ = 

r.(g'; O) = ^ / ^,D{q' - g)tr. [S{-q)S{q)] T,{q; 0). (4.2) 
This is equivalent to the DSE given in Eq. (|2.34 ) for the scalar self-energy amplitude provided the 



chiral limit is taken; in this regard it is useful to note that So{—q)So{q) = {q'^Al + Bq)~^. Thus there 
is the massless solution Th{q; 0) oc Po(9^) for each PS meson. The above argument is unchanged if we 
consider the SU{3)f singlet with A° = Ij/y/S. However this (r]') is not a Goldstone boson and this 
circumstance (f/A(l) anomaly) must be addressed beyond the ladder BS and rainbow DSE level. 

The above argument covers only the pseudoscalar Dirac component of the pion BS amplitude. A 
full analysis should take into account the general form of the J^*" = 0^+ BS amplitude which is 
(LLewellyn-Smith, 1969) 

f^(g; P) = tl,T{E^{q; P) + 7 ■ PF^{q; P) + l ■ qq ■ PG^{q; P) + [7 ■ g, 7 ■ P]H^{q; P)} , (4.3) 

where the invariant amplitudes are functions of P^, and {q ■ P)"^. The last term does not arise in 
the present discussion because 0"^,^ is not part of the set of matrices A^ due to Feynman-like choice 
of gauge for the effective gluon 2-point function. The above result is that considering only the chiral 
limit Et,, one flnds ET,{q; 0) oc Boi^q"^). 

In general, without relying upon the explicit rainbow-ladder truncation, the Goldstone theorem for 
the pion, along with identiflcation of the proper normalization of the chiral limit BS amplitude, can 
be obtained from consideration of the axial vector vertex r^(g; P). The inhomogeneous BS equation 

for r^(g; P), driven by the bare term i'j^'j^^, has a pole when the homogeneous version has a solution. 
To be speciflc, the ladder BSE equation for this vertex is 

rliq'; P) = 11,1 J- - \ J -^,D{q' - g)7.5(g+)f J(g; P)S{q.)i^. (4.4) 

If we contract with P^, and note that P^r|^(g; P) is a pseudoscalar, we see that Eq. ( [4. 4] ) generates the 
inhomogeneous version of the ladder BSE for the pion BS amplitude. As we have discussed before, 
one may introduce eigenfunctions of that homogeneous kernel and its eigenmode expansion generates 
the pion propagator. The mass-shell BS amplitude provides the residue at the physical pole. Thus, 



33 



there is a pion pole in T^^{q; P) (Jackiw and Johnson, 1973), and near the mass shell of the chiral limit 
pion we have 

rliq; P) = - f.(g; P) (-^) (-^P,/.) + f f ^^(g; P), (4.5) 

where the last term is regular there, TT^{q] P) is the properly normalized pion BS amplitude, and the 
pion coupling to the point axial field is 

- ^P,fJ^, = tr / ^ [r;(g; -P)S{q^)i^,^J-^S{q^)] . (4.6) 

These two results are true in general; they are not limited to ladder approximation. Meanwhile, the 
isovector axial WTI identity 

- P.f^(g; P) = S-\q^)^J- + lJ-S-\q^) (4.7) 

also reveals the leading behavior of r^(g; P) at the pole (for the longitudinal component) in terms of 
the dressed quark propagator. In combination with Eq. ([4.5|) , this provides a means for connecting 
the chiral limit pion BS amplitude directly to the quark propagator (Delbourgo and Scadron, 1979). 
In particular, from Eq. ( [4.7|) we have 



P.fliq- P) = 75r ( -Po(g') + 11 ■ P^^ + il-qP- qA',{q^) + 0{P') \ . (4.^ 



From this and Eq. (4^), the amplitude E^^ for the 75 component of the properly normalized chiral 
limit pion BS amplitude is identified as 

EM:P = ^) = ^ (4.9) 



where the weak decay constant is given by Eq. ( |4.6| ). The other terms in Eq. ( ^.8| ) can arise from a 
combination of the pole term and the regular term. Further analysis is needed to be able to express 
the BS amplitudes P^- and directly in terms of the quark amplitudes A and B. The result in 
Eq. ([4.9| ) is essentially the Goldberger-Treiman relation for quarks (instead of nucleons in the original 
development). The rainbow-ladder approximation has the important property of preserving the above 
manifestations of Goldstone's theorem. This is due to the fact that contraction of Eq. ( [4.4| ) with P^ 
will produce exactly Eq. ( [4.7|) only if the propagators in Eq. (f4.4|) are DSE solutions in rainbow 
approximation. 



It is clear from the above that if the dynamically coupled pion amplitudes P^, P^ and G.^ are expressed 
in a normalization where the chiral limit of P,r is Po, then the normalization constant that must be 
divided into each of them for a physical amplitude is simply A^^r = /tt- This dual role for f^, (measuring 
the weak decay and providing the inverse of the field normalization constant in the above sense) is also 
evident in the 0{p^) chiral lagrangian where it is the only parameter. The dual role of is a useful 
device in model calculations. With the asssumption that the pion is dominated by the 75 amplitude, 
use of Eq. ( [4.9|) in the chiral limit BS normalization condition (see Eq. ( p.61|) ) 

= 1 (4.10) 

P'2=K2=0 



d 

ap2 



tr 



r^(g; -ir)S'o(g+)r^(g; K)So{q- 
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f: = 12 / 7^^o(g^) I'^Kas - aya's) + ayas]. (4.11) 



will produce a value for fj^ without explicitly computing the weak decay process. However, from the 
definition of given in Eq. ( [4.6|) , and again assuming that the pion is dominated by the 75 amplitude, 
one may also obtain in the chiral limit. 

The two resulting expressions for /^r do not generally agree and some caution is needed. In practice, 
the difference is usually found to be at most 10%. The reason they do not agree is that an approximate 
BS amplitude (just the 75 amplitude E^^) has been used; the equivalence of the two approaches to 
/tt comes from analysis of Eq. ( [4. 41) and Eq. ( [4.7|) which are compatible only if the full structure 
of the pion BS amplitude allowed by both of them is employed. The size of the contribution from 
the sub-dominant pion BS amplitudes F.^ and Gj, to various observables is a question that deserves 
further attention. Much work on modeling of the pion omits them. Many studies of the ladder BS 
equation (Aoki, Kugo and Mitchard, 1991; Jain and Munczek, 1993; Burden, Qian, Roberts, Tandy 
and Thomson, 1996)) confirm that the dominant amplitude in Eq. ( |4.3| ) is E.,^{q; P), especially as far as 
contributions to the mass are concerned. However and will contribute to /tt via Eq. ( |4.6| ) even 
in the chiral limit and the contribution can be as large as 35% (Burden, Qian, Roberts, Tandy and 
Thomson, 1996). For the special case oi A ^ 1, Eq. ( |4.8| ) indicates that the sub-dominant amplitudes 
of produce at most an 0{P^) contribution to Eq. ([4.6|) and therefore will not contribute to the 
chiral limit fj^. We also note that when A — > 1, Eq. ( 4.11| ) reduces to the formula obtained by Pagels 
and Stokar (1979). 

To go beyond the above symmetry-based observations and calculate the finite pion mass generated 
by a finite bare quark mass, a dynamical formulation must be specified. The Goldstone theorem at 
the nearby chiral point should be preserved by the dynamics and the ladder BSE and rainbow DSE 
combination therefore suggests itself. It is difficult to formulate general higher order descriptions 
of the BS kernel and the dressed gluon-quark vertex for the DSE that also preserve the Goldstone 
theorem. At order (7^ this has only recently been achieved (Bender, Roberts and V. Smekal, 1996) We 
will limit our considerations to the rainbow-ladder level. 

Consider now the PS octet modes for the physical situation of finite bare quark masses. We shall 
continue to make the dominant covariant assumption for simplicity and we shall proceed from the 
action given in Eq. ( |2.63|) . The free or quadratic part of this action for the PS octet can be written 



where 



S2[<P] = \j^. <PliP) ^b\P') MP), (4.12) 

A,-^(P^) = tr / ^ [^l,XlS{q^)^lr^bS{q-)] fg(g; P) + d'r^^^^. (4.13) 

Here the quark propagator S{q) can be considered a diagonal matrix in fiavor space. We can, at 
this stage, consider the TbiQ', P) to be ladder BS amplitudes that are not yet physically normed. 
Evidently A^^(P^) is an (unnormed) inverse propagator for the relevant local fields. Here g± = g ± ^ 
and Afe is the relevant flavor matrix. For example, with Nj = 3, we have A;, (A^ + i\^)/\/2 for tt"^, 
Ab ^ (A^ + for K+, and A* for r/; with Nf = 2, we have A^ ^ (r^ + iT^)/V2 for 7r+. 

We note that A^^ is a functional of f;, and that 6A^^[tb]/Stf, = reproduces Eq. ([4.1| ), the lad- 
der BS equation (projected onto the i'y^Xb channel) for the amplitude Tj,. Hence A^" (P^) must 
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have a zero at the corresponding invariant mass and the general momentum dependence can be 
written as A^^(P^) = (P^ + m^)Zj~^(P^). An equivalent observation is that the produced meson 
mass is a dynamical function mf,(P^) due to qq substructure. The field renormalization constant is 
Zb = Zbi^P"^ = —ml) and a physical normalization (unit residue at the mass-shell pole) is produced by 
absorbing 1 / \^Zf) into the fields. The physically normalized BS amplitude is then the dimensionless 
quantity Tb{p] P) = \/ZbTb{p] P)- This procedure satisfies the canonical Bethe-Salpeter normalization 
condition (Itzykson and Zuber, 1980). 

The functional A^^[f5] offers a very useful strategy for approximate solutions of the ladder BS equa- 
tion. One may use a variational trial function Tf){p) to minimize the node position in the resulting 
A^^(P^). This is equivalent to minimization of the mass functional mb[r6]. This technique was intro- 
duced by Cahill, Roberts and Praschifka (1987) where it is also assumed that A^^(P^), in the required 
time-like region P^ < 0, can be extrapolated linearly from knowledge of its value and derivative at 
P^ = 0. This is equivalent to ignoring momentum dependence of Zi,{P^). Then one has immedi- 
ately ml ^ ZbA^^{0), and Z^^ = A^^ '(0) where the prime indicates a derivative with respect to the 
argument. From Eq. ( [4.13|) one has 



ml[Tb] =ZbtTj J^,[^l,XlS{q)^J,hSiq)] f^(g) + ^Z, / d'r^, (4.14) 

and 

2b-' = ^tr / ^[^75Al5(g+)^75A,5(g_)] f^(g) 



(4.15) 

P2=0 



When f f,(g) — BQ{q^) for the pion, Eq. ( |4.15| ) indicates that y/Zbth^q) is physically normed, and from 
Eq. (|4.9| ) we see that Zj^^ is consistent with this present level of approximation. The assumptions 
underlying the above two expressions are best suited to small masses. Numerical investigations for a 
wide range of meson states were presented in Praschifka, Cahill and Roberts (1989). 

The above variational approach requires a specification of the model gluon propagator. To approx- 
imate the pion mass and normalization constant, an approach that has been used often is to set 
Tn{q', P) ~ B{q^) from the DSE solution at finite m. This has the correct chiral limit and allows the 
last term of Eq. ( [4.13|) involving D to be eliminated by requiring that the GMOR relation (Gell-Mann, 
Oakes and Renner, 1968) 

mlf^ ^ -2m{uu) (4.16) 
be satisfied. Since here A^^(O) = m^f^, this approach yields 

A-\P') = -2m{uu) + trj -0-^[Si-q)S{q) - S{-q+)S{q^)] B\q'), (4.17) 

The conventional expression for the chiral condensate is (uu) = —12 ]§^'-^si^)- ^he appropriate 
scale for such an approach is usually taken as /i = 1 GeVl Mitchell and Tandy (1997) have used 
this procedure to model the composite pion propagator in a study of p — mixing and mass splitting. 

A more natural approach, that does not force the GMOR relation, follows from using the chiral limit 
approximation for the BS amplitude, i.e. rTj{q;P) ~ Bo^q"^), while keeping the effect of a non-zero 
bare mass in the quark propagators. This gives the pion mass to first order perturbation from the bare 
quark mass. The last term of Eq. ( [4.13| ) may then be recast solely in terms of the chiral limit quark 
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propagator. The chiral limit DSE (see Eqs. ( |2.32D and ( |2.34] )) in the form B^ir) / D{r) = |trsS'o(r) 
leads to 



- / t/V 



D(r) 



2Nr_ tr. 



Then Eq. ( [4 .131) for the pion becomes 



2Nr, tr. 



[So{q)Bo{q') - S{-q^)S{q^) Bliq'^ 



(4.18) 



(4.19) 



where a subscript indicates a chiral limit quantity. Now A^^ is fully specified in terms of the quark 
propagator. Numerical evaluation of Eq. ( [4.19 ) may be used to identify the pion mass through the 
solution of A^^(P^ = —m^) = and the normalization constant may be evaluated from 



d 



dP' 



(4.20) 



P2 = 



These expressions require that the extrapolation of the quark propagator into the complex quark 
momentum plane be well-defined. This is the case with the confining entire function parameterizations 
that have been available in recent years. Before that development, a common approximation was to 
use A~^(P^) ^ /^(-P^ + ml) so that mlf^ ^ A~^(0) and ^ A~'^ '(0). Because of the small mass 
of the pion this is quite accurate and the resulting explicit pion mass formula (Cahill and Gunner, 
1995b; Frank and Roberts, 1996) is 



mlf:^8N,J^^Bo{q')[a'siq 



Baia^) 



(27 



B{q^) + m 



<^s{q% 



with 



J TT 



d 



dP 



2 2iVetr, 



(2vr) 



P2=0 



More explicit expressions are 

. f 2 ^ 24 / A B,{q^)[q^A,{q^)jV,{q^) + B,{B', + l)] 
m^Ur^ZmZ^J ^2^^^ q^Al + Bi 

to first order in m with /' = df /dm and subscript labelling chiral limit quantities; and 



(4.21) 



(4.22) 



(4.23) 



N r°° 

1-2 - . / 
J n 



iir"^ Jo 



J dss Bo{sY (ay - 2 [crsa'g + savcr'v] - sas - s'^ay 



(4.24) 



where / = //" — (/')^, and s = p^. 



We note that in this treatment, the finite value of the pion mass is tied directly to the explicit 
chiral symmetry breaking behavior of the quark propagator. It is determined by the response of the 
amplitudes B{q^) and crs{q^) to a small current mass. In realistic studies of the DSE, it has been 
found that the response of (Js{q^) at low momenta to a current mass is negative. That feature has 
been found to be an important element in generating pion mass via Eq. ( [4.211 ). Numerical studies 
of the BS and DSE equations (Frank and Roberts, 1996) reveal that, for typical finite current quark 
masses, TT,{q; P) ^ B^^q^) is indeed a good approximation. It is also found there that the pion mass 
formula in Eq. ([4.21[ ) accurately reproduces the pion mass found from the BSE solution. 
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We also note that the pion mass formula in Eq. ( [4.21[ ) is different from the GMOR relation in Eq. ( [4. 161 ) 



which is not exact in QCD. As a measurement of explicit chiral symmetry breaking, the leading term 



on the RHS of Eq. ( 4.21 ) is first-order in the current mass, it gives the response in the pion state 
(subject to the approximation r^(g;P) ^ Boi^q"^)), and it is therefore determined by vacuum quark 
properties. This is consistent with the general arguments of Gell-Mann, Oakes and Renner (1968). To 
the extent that one must make further approximations to obtain Eq. ( |4.16| ), the relation Eq. ( |4.21|) 



can be thought of as containing corrections defined by the above treatment within the GCM. Hence, 
as a summary of properties of the composite pion mode from either the GCM or the DSE approach, 
both truncated to 0*'' order in meson loops, the inverse propagator A^^(P^) given by Eq. ( [1.19| ) is 
quite accurate and practical near the mass-shell. One qualification is that only the amplitude for the 
dominant covariant 75 has been accounted for therein. As remarked earlier, there are recent indications 
that the common omission of the pseudovector pion component misrepresents some aspects. More 
work on this topic is warranted. 



The extension of Eq. ([4. 19]) to other members of the PS octet is straightforward although one would 



not expect that such a reliance upon first order effects from explicit chiral symmetry breaking would 
be as accurate. From Eq. ( [4.13|) for the we have 



A],\P') = 2N^tT,l^jB^,{q')S^{q) - B^^iq^f S-{-q^)S^{q.)] . (4.25) 

where is the propagator for flavor / and contains 7^ 0. The resulting mass formula is thus 

mlf^ ^ 2N^ trj ^,[B,{q') So{q) - B',{q') S-{-q)S%q)] . (4.26) 

Similarly for the pure flavor octet rj we have 

m^f^ ^ 2N^ tr. / ^ [B,{q') So{q) - \Bl{q') [S\-q)S-{q) + S\-q)S\q) + AS\-q)S\q)]] . 

(4.27) 

Again ^ A^-i '(0) for each case. The results to first order in the bare masses can be expressed 



2 €2 / I \ 2 c2 I I \ 2 £2 + '^d + 4ms) / , c-,o\ 

= i^u + md)p ; m^fj^ = (m„ + rusjp ; m^/^ = p, (4.28) 



where p can be identified from Eq. (^4.23|) and it should be noted that the constants are iden- 
tical in the chiral limit. With degenerate u and d masses, the Gell-Mann-Okubo mass formula 
ml + 3m^ = 4m|; follows immediately. 



4.2 Chiral Qbservables 

In principle, if a systematic approach is taken to produce from the GCM an effective action for all 
the Goldstone modes, and an expansion in momenta and masses is made, then the results will be 
predictions for the coefficients of the effective action employed in chiral perturbation theory (ChPT). 
This is because the GCM respects dynamically broken chiral symmetry and hence will share with 
ChPT all the consequences of this. The same relationship is in principle shared between QCD and 
ChPT. Given the difficulty of deriving the low energy constants of ChPT directly from QCD, they 
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are invariably determined as phenomenological parameters from a canonical set of data in the form of 
pseudoscalar meson decays, mass relationships, etc. Then predictions are made for other low energy 
experimental quantities. One may investigate whether a model field theory such as the GCM can 
provide some insight into the microscopic content of the low energy constants of ChPT. Although 
some threshold chiral observables including the vr — tt scattering lengths and the pion charge radius 
have been studied within the GCM, some caution should be exercised in comparing extracted low 
energy constants with values in use in ChPT. This is because constants that control the coupling of 
point interpolating fields with pion quantum numbers, as defined by the form of an effective field theory 
such as ChPT, imply an ordering of contributing mechanisms that does not necessarily correspond 
to the ordering associated with the coupling of distributed pionic qq correlations produced by an 
underlying theory at the quark- gluon level. 

The bosonization of the GCM to produce the meson action in Eq. ( ^.64[ ) provides an illustration. The 
fields are ladder qq objects, the couplings are expressed as distributed vertex functions, and this is 
only the tree-level result. In principle, one may integrate up through the loop expansion, and integrate 
out the non-Goldstone mesons, stopping when the effect on a given threshold observable has (hope- 
fully) converged. For an observable at given order in external momenta, the relative importance of a 
particular type of coupling or a particular order in loops can be expected to be influenced significantly 
by the treatment of the nonlocalities or finite size effects. These are characteristics of the non-point 
nature of the states. The meson loop integrals in the GCM are finite due to intrinsic substructure in 
the fields and hence in the vertex functions. To compare parameters in detail with ChPT, it matters 
whether one makes a derivative expansion of the GCM after quantum loop integrations or before. In 
the former case, the final importance of loops is a characteristic of the GCM; the latter case is not 
really consistent. Even a localization of the tree- level couplings means that the distributed nature of 
the dressed quark substructure of the bare mesons and interactions is ignored and this can change 
the relative contributions to observables. Unambiguous comparisons of the GCM and ChPT must be 
made through the calculated observables. Most of the investigations within the GCM have been at 
tree-level because loop effects are technically difficult due to nonlocalities. As we will discuss later, 
investigations of pion loop contributions within the GCM exist for a few selected observables and 
indicate small contributions, typically at the 10% level. Much more work remains to be done in that 
direction. 

In two recent works, low energy chiral observables were considered to fourth order in gradients and 
explicit chiral symmetry breaking masses. The usual chiral power counting is observed (Gasser and 
Leutwyler, 1984, 1985; Donoghue, Golowich and Holstein, 1992). The first derivation and analysis of 
the Goldstone modes and chiral coefficients to 0{p'^) from the GCM was made by Roberts, Cahill and 
Praschifka (1988a). The identification of the sums of quark loop integrals that contribute to masses, 
decay constants, three 0{p^) coefficients, and the coefficient of the anomalous Wess-Zumino (1971) 
five-pion term was made in that work. This is a complicated task and the interested reader is referred 
to that paper for details. The implications for vr — vr scattering from the terms through O(p^) were also 
drawn in that work and in particular the Weinberg (1966) result for the S-wave scattering lengths was 
well reproduced. The work of Roberts, Cahill, Sevior, and lannella (1994) extended the application 
for 7r — 7r scattering to include the implications from the three 0{p'^) coefficients considered and to 
employ a more realistic dressed quark propagator model. 
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To summarize the approach, we first note that in the chiral hmit of zero bare quark mass, when 



S-\x ~y;U)=j- d,Ao{x -y) + Boix - y)V 



(4.29) 



solves the DSE for = 1, it also does so for any constant unitary matrix V of the form (Praschifka, 
Roberts and Cahill (1987a)) 



V = exp 



(4.30) 



Here the 0^ are constants and Pr^l the standard right and left-handed projection operators (1 ± 75)/2. 
(One may include in Eq. (^4.30|) the A'' oc 1/ term for the singlet r/' but we restrict attention to 
the Goldstone flavor octet.) With a set of constants (jJ^ chosen to select one out of the degenerate 
vacuua, fluctuation fields will be incorporated via the coordinatization 0^ — > 0^(2^). Without loss of 
generality, the vacuum point about which to expand the action can be taken to he V = 1 = U . In 
terms of the chiral field U{R) defined by 



U{R) = exp 



(4.31) 



the action, given earlier by Eq. ( |2.26| ), may be written for the goldstone octet as 

S[U] = -TiLnS-^p] + ^ti j d^rd^R { (Sy(-r) + 5o(-r)f/^(i?)f/(i?)) 5(r; 1)} 



(4.32) 

where is the Dirac vector part of the quark self-energy. Although the second term of this action 
is irrelevant under the present circumstances where the constraint t/^t/ = 1 is operative, it is relevant 
in general. One example would be the generation, along these lines, of the composite qq version of 
the linear sigma model where WU = with x the chiral radius field variable (Cahill and Roberts, 
1985; Frank and Tandy, 1992). 



The physically relevant action 5[f/] = — 5[1] has a real part 



Sr[U] 



-TrLn 
2 



s-\u)] s{iys{i)s^\u)) . 



(4.33) 



With a perturbative introduction of the bare quark mass matrix M^2 = diag{mu,mj)^2, the low 
energy form through fourth order in gradients and masses for the terms relevant to vr — vr scattering 
was derived by Roberts, Cahill, Sevior, and lannella (1994) in the form 



Sr[U] = 
-N^Ki tr 



d^x I ^tr 



+ 



(qq). 



-tr 



+ N^K2 tr 



2If-U- U^) 



N,Ks - tr 



(4.34) 



Here the field equations of motion at second order, which could be used to eliminate the explicit 
appearance of the Ki term, have not been applied. The coefficients Ki, K2, and are given by 
convergent quark loop integrals which require only the chiral limit dressed quark propagator. In 
terms of the standard chiral coefficients Li_iq (Gasser and Leutwyler, 1985), the form in Eq. ([4.34|) 
corresponds to 2Li = L2 = L3 = Nc{Ki — K2 — K3), and Lg = —NcKi/4, with the remaining 

Li = 0. 
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The imaginary part of the action S[U] involves an odd number of fields and corresponds to all of the 
anomalous interactions. The lowest order contribution is fifth order in chiral fields and is 

tSi[U] = / {{U{x) - U^{x)) d,U^{x)dM{x)dpU\x)d^U{x)} , (4.35) 

where A is given by a convergent integral involving self-energy amplitudes Aq{p'^) and i?o(p^) and is 
given in Praschifka, Roberts and Cahill (1987a). It is remarkable that this quark loop integral gives 
exactly A = 1 independently of the detailed form of the quark propagator. Expansion of Eq. ( |4.35|) 



to lowest order in pion fields yields exactly the five-pion form of the Wess-Zumino term. This is a 
characteristic feature in that the form and strength of anomalies required by symmetries are properly 
embedded in the GCM. It is only necessary to properly maintain the symmetry constraints among 
the various dynamical quantities at the quark level. In the above case, the A = 1 result is critically 
dependent upon the symmetry result that -Bo(p^) is both the chiral limit 75 piece of the pion BS 
amplitude (as required by the axial Ward identity) and the quark scalar self-energy resulting from the 
dynamical breaking of chiral symmetry. Further studies of the anomalies generated by the GCM have 



been made by Roberts, Praschifka and Cahill (1989) where it is shown that extension of Eq. ( 4.35 ) to 



include vector mesons at the beginning of the derivation reproduces all of the anomalous interactions. 

The first two terms of Eq. ([4.34|) , the "kinetic" terms for the chiral field U{x), generate the familiar 
Weinberg result (Weinberg, 1966) for the invariant vr — vr scattering amplitude 

2 

g — ui 

Awis,t,u) = ^^ (4.36) 

J IT 

that provides the leading low-energy behavior. The scattering amplitude contributions from the 0(p^) 
action are linear in -R'l-s and were calculated from the dressed quark substructure of the pions by 
Roberts, Cahill, Sevior, and lannella (1994). That work employed model three-parameter forms of 
the quark self-energy amplitudes designed to reproduce the characteristic features that emerge from 
all studies of the quark propagator via the DSE. It was demonstrated that the infra-red structure of 
the dressed quark propagator in the GCM may be modeled to provide a sensible description, not only 
for /tt, m^r and which were the principal concerns to that stage, but also for the main features of 
the vr — vr scattering lengths aj. A reasonable description of the principal qualitative features of the 
TT — vr partial wave amplitudes for / < 2 and E < ArriT^ was also obtained. In this work, the leading 
asymptotic behavior of the propagator amphtudes corresponded to the outcomes from the operator 
product expansion and QCD renormalization (Politzer, 1976; 1982). In this and subsequent studies, 
it is observed that the logarithmic corrections inherent in such UV behavior, and associated with the 
anomalous dimension of the propagator, do not provide significant numerical contributions to physical 
quantities. 



In Table we show results for low energy pion observables including the vr — vr scattering lengths 
and associated chiral coefficients i^i-a from a tree-level calculation within the GCM that uses the 
parameterized confining quark propagator model discussed in Sec. ^. The original GCM calculation 
of the 71 — n scattering lengths (Roberts et al, 1994) employed a quite different propagator param- 
eterization; we choose the more recent calculation here for consistency with other elements of this 



review. The collection of quantities shown in Table |4.1| were fitted with the propagator model given 



in Eqs. ( |3.23| ) and ( p.24|) and with the parameters shown in Table [4.2| . Although the results for ob- 



servables are very similar to those that also appear in the more complete set shown as Table |3.2|, the 



relevant low energy constants are extracted only for this present calculation. The point meson limit 
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Table 4.1: A comparison between the low-energy vr observables calcu- 



lated at tree-level in the GCM using the parameters of Table [4.2| and 
their experimental values. The superscripts indicate the reference that 
experimental values are taken from: 1 - Particle Data Group (1996); 2 
- Amendolia et al. (1986a); 3 - Pocanic (1995) and Gasser (1995). The 
value for (qq) quoted as experimental is that typically used in QCD sum 
rule analysis. (Adapted from Roberts (1997)). 







T"^/ vnpri TTi pn f 


J vr 


0.0924 GeV 


0.0924 ±0.001^ 


_L 


U.Z4:i LreV 


U.zzU ± U.UJdU 




U.UU40 LreV 


n nns -i-n nnvi 
U.UUo itU.UU 1 




U.ioy LtC V 


v. loo 




0.547 fm 


0.663±0.0062 




0.505 (dimensionless) 


0.500 ±0.018^ 


a^ 
"o 


0.170 


0.26±0.05^ 


al 


-0.0481 


-0.028 ±0.012^ 


a\ 


0.0301 


0.038 ±0.002^ 


al 


0.00149 


0.0017 ±0.0003^ 


al 


-0.000214 


0.00013 ±0.0003^ 


Ki 


0.000702 




K2 


0.000867 






0.000875 





of this GCM calculation, which corresponds to the NJL model, produces Ki = l/967r^ = 0.00106 and 
the comparison with Table ^TT] indicates roughly that one can expect a dependence upon finite size 
effects associated with meson substructure at the 20% level. 

The values of the vr — vr scattering lengths shown in Table ^]l| are close to the current algebra or soft 
pion theorem results (Weinberg, 1966). The latter produces Oq = 0.16 for example. In ChPT up to 
one pion loop at C(p^), this value is corrected upwards by about 30% to become 0.20 ± 0.01 (Gasser 
and Leutwyler, 1983). Indications from two-loop results at 0{p^) are that the next correction is small 
(giving 0.217) and that it presently appears to be unlikely that ChPT will reach 0.26 (Ecker, 1996). 
The pion loop contribution to aj] has not been calculated in the GCM. Pion loop contributions in 
the GCM can in principle be obtained for each chiral observable independently and this is clearly a 
topic that deserves attention in the future. A meaningful comparison with one-loop ChPT clearly 
requires that loop contributions in the GCM be obtained for a number of observables at least equal 
to the number of independent low energy coefficients at a given order in ChPT. This is because a 
general observable is a linear combination of these coefficients. For example, for Nf = 2, one would 
need as many observables as needed to determine the Z,, z = 1 — 6 of Gasser and Leutwyler, (1984). 
The quantity a^ is a linear combination of li,l2,h that are determined by fitting D-wave scattering 
lengths and SU (3) mass relations. 

With Nf = 3, Frank and Meissner (1996) have studied the 0{p^) chiral coefficients Li_io that can 
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be produced at tree-level in the GCM. They were particularly interested in the relationship between 
the values of the coefficients and the chosen infra-red behavior of the phenomenological 2-point gluon 
function employed in the GCM. In order to facilitate a chiral invariant derivation, they adopt the 
common procedure of coupling the quarks to an external field iplx) = [s+i'j^p) which finally will be set 
to the bare quark mass matrix (Gasser and Leutwyler, 1983; 1985; Donoghue, Golowich and Holstein, 
1992). The x-dependence of iIj{x) helps prevent unphysical results for some of the coefficients and 
also enters into the equation of motion which is obtained at second order and employed to implement 
relations among the coefficients. They use 



(4.37) 



with r = X — y and R = {x + y)/2, and show that the produced tree-level structure of the GCM, to 
0(/), is 



d X 



^tr 
4 



- ^tr 
4 



Li (tr Ud^U)id^U^) 



- L-i tr 



{d^U){dM^)] - tr Ud^U){dM^) 



U tr 



{d^u){d^u^){d,u){d,u^)\ + l, tr \{d,u){d,u^){Ux^ + xu^) 



(4.38) 



Lg tr IxU^xU^ + Ux^Ux^ 



where x{^) = ~2(gg)'?/'(x)//^ is of order 2 (in mass) and the remaining trace is over flavor. 



Table 4.2: Parameters for the quark propagator model in Eqs. (|3.23| ) and 
( p.24| ) that produce the ffi shown in Table ^]l| for the low energy pion 
quantities. Note A = 10~^ is not varied. At the finite bare mass values. 

Cm. = 0. 



u/d- quark 


m 


0.008 


bo 


0.181 


hi 


2.89 


b2 


0.555 


bs 


0.184 


Cm=0 


0.119 


D 


0.159 GeV 



The chiral coefficients, classified according to Gasser and Leutwyler (1985), are produced in this work 
as convergent quark loop integrals. Of the complete set of C(p^) coefficients labelled Li_io, only those 
shown in Eq. ( 4.38 ) occur at tree-level in the GCM. The coefficients Lg and Liq refer to external vector 
and axial vector fields and are not considered here. Application of the SU{3) relation (Donoghue, 
Golowich and Holstein, 1992) between the terms fourth order in U produces L2 = 2Li. Apart from 
the special case of L7 which becomes finite when the Ua{^) anomaly is accounted for (Witten, 1979), 
these relationships agree with what is expected in the large Nc limit of QCD (Donoghue, Golowich and 
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Holstein, 1992). This is presumably related to the fact that the large Nc limit of QCD corresponds 
to a ladder description of mesons if the three and four gluon couplings are neglected as they are in 
the GCM (Frank and Meissner, 1996). A complicating element for considering a l/N^ expansion for 
the GCM is that a phenomenological two-point gluon function has an unknown dependence upon 
Nc- The neglect of higher mass states, explicitly excludes intermediate states of pure glue which are 
"iVc suppressed" . Since the GCM action at the quark level is consistent with explicit and dynamical 
breaking of chiral symmetry, the complete set of coefficients must arise if higher mass states and non- 
Goldstone meson modes were not simply discarded but integrated out, and if loop contributions were 
implemented. The minimum steps of this type needed to generate the complete set of coefficients in 
the GCM is not known. 



In principle, the four independent coefficients (which can be taken as Li, L3, L5, Lg ) produced by the 
GCM at tree-level in Eq. ( |4.38D depend on the details of the dressed quark propagator and hence 
implicitly on the form of the effective quark-quark interaction D. With g'^D{q^) = A7ra{q'^)/q'^, Frank 
and Meissner (1996) have investigated three different two-parameter models for a{s), namely 



ai[s) 
"2(5) 
"3(5) 



s/A 



+ 



Tid 



ird 



Tcd 



4A2 ln(s/A2 + e) 

sx^ 1 
+ 



s2 + A ln(s/A2 + e) 
1 + xe""/^ 



(4.39) 



ln(s/A2 + e) 



where s = and d = 12/(33 — 2Nf). Each of these forms incorporates the one-loop perturbative 
result for large s and extrapolates differently into the low-momentum region. The two low-momentum 
parameters, x ^i^nd A, are varied with the pion decay constant held fixed at fj^ = 86 MeV. This value 
was considered to be appropriate for a chiral limit calculation at zero-momentum rather than the 
physical value of 93 MeV appropriate to the physical mass-shell. The results however are not very 
sensitive to this small difference. By fixing f^^ the overall scale of D^SB is fixed. The remaining 
independent parameter is associated with the matching scale to the perturbative form. 



The first model, ai in Eq. ( ^.39| ), corresponds to that in Eq. ( p.l3| ) which, as discussed earlier, has 
been used in many previous investigations within the GCM (Praschifka et al, 1989). Here Frank and 
Meissner (1996) use slightly different parameters due to their chosen /,r value. The small A limit of 
this model simulates the infrared-dominant behavior which is known to model confinement (Burden 
et al, 1992). The infrared contribution to the second model, 02, generates a. singularity in the 
quark-quark interaction D{q^) in the limit as A — ^ 0. Such a singularity has also been considered 
previously as a model of confinement (Marciano and Pagels, 1978). The 1/g^ form falls much slower 
than the Gaussian in ai, and hence matches the UV form at a much higher scale. Finally the third 
model 0:3 was chosen here to be structurally different from ai and 0:2 in order to further illustrate the 
independence of the results to the details of the low-momentum parameterization. 

In all three cases the same pattern is observed: the coefficients Li and L3, which are important for 
TT-vr and K-K scattering, are nearly independent of the form of a{s) and therefore on the form of the 
quark-quark interaction, provided that the integrated strength of a{s) is fixed by On the other 
hand, the mass dependent coefficients, L5 and Lg, are moderately dependent on the actual form of the 
interaction. For example with L5, in order to reproduce the empirical value, forms of a{s) with a small 
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matching scale, i.e. which are relatively strong in the infrared region, are required. This observation 
agrees with the findings by Meissner (1994), where it is shown that infrared- dominant forms of the 
quark-quark interaction are required to achieve convergence of the chiral series in the strange quark 
sector. This also confirms the success of the "delta-function-plus-tail" models (obtained for example 
from ai in the limit A ^ 0) in describing chiral observables as most recently observed by Frank and 
Roberts (1996). 

Frank and Meissner (1996) also find that the results for the fourth-order coefficients are rather insensi- 
tive to the asymptotic UV tale of a{s); even omitting this tail completely gives no significant changes, 
again provided that /^r is fixed. Typical results are illustrated in Table |4.3| . The mixed quark-gluon 
condensate — fi's(^o"G'g)i cev^) as recently calculated by Meissner (1997), is also shown in Table |^ . 



Table 4.3: Low energy chiral coefficients, the quark condensate, and 
the quark-gluon mixed condensate from the GCM at meson tree-level. 
The interaction described by ai(s) in Eq. ( [4.39| ) is employed with 
A = 0.2 GeV, Nf = 3. Adapted from Frank and Meissner (1996) and 
Meissner (1997). Empirical values shown in the last row are taken from 
Meifiner (1993) and Bernard, Kaiser and Meifiner (1995) for the chiral co- 
efficients, and from the quenched lattice of Kremer and Schierholz (1987) 
for the condensates. 
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0.81 


-4.03 


1.66 


0.83 


0.020 


1.55 


160 
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0.82 


-4.39 


1.13 


0.84 


0.002 


1.45 


151 


432 


0.85 


-4.43 


0.97 


0.88 






225 


402 - 429 


0.7 ±0.5 


-3.6 ±1.3 


1.4 ±0.5 


0.9 ±0.3 



4.3 Other Mesons 

Beyond the light Goldstone bosons it is necessary to use the BSE and solve with an approximate 
kernel. Most common is the ladder approximation and in the GCM or DSE approach one must adopt 
an Ansatz for the effective gluon 2-point function D{q). It is not possible to use chiral symmetry to 



avoid this in the manner of Sec. [4.1| . To be consistent with what is arrived at naturally in the GCM, 
and consistent with preservation of the Goldstone theorem, the employed dressed quark propagators 
should be from rainbow approximation DSE using the same gluon function D{q). The most extensive 
and phenomenologically successful spectroscopic studies in the rainbow-ladder framework are those 
of Jain and Munczek (1993) in which the quark DSE is solved numerically for spacelike-p^ using a 
model gluon propagator. In Landau gauge the behavior of the gluon propagator is constrained by 
perturbation theory for /c^ > 1 — 2 GeV^ (Brown and Pennington, 1989) and one models the infrared 
behavior, which is presently unknown. Such studies have the ability to unify many observables via the 
few parameters that characterize the behavior of the model dressed-gluon propagator in the infrared. 

In this approach, solving the meson BSEs is complicated by the fact that these equations sample 
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the dressed-quark propagator off the spacehke-p^ axis. This difficulty was circumvented by Jain and 
Munczek (1993) through a derivative expansion of the dressed-quark propagator functions, ^(p^) 
and -B(p^), and estimating the error introduced thereby. This, however, obscures the discussion 
and exploration of the role of quark and gluon confinement. The successful use of parameterized 
entire function quark propagators for the light Goldstone bosons as described earlier has recently 
led to several investigations, roughly within the framework of the GCM and DSE approach, to see 
whether the constraints from pion and kaon physics contained therein can also constrain the properties 
of a wider class of light mesons. We describe briefly approaches by Cahill and Gunner (1995a), 
and by Burden et al. (1997), both of which are organized around use of a separable Ansatz for the 
ladder-rainbow BSE kernel. The purpose in choosing such an unusual representation is so that, when 
employed in the DSE it can easily be arranged to reproduce specified dressed quark propagators 
previously fitted to pion and kaon physics. The question is then asked whether subsequent use of this 
kernel Ansatz in the BSE can produce the dominant physics of certain other mesons such as the light 
vectors and axial vectors. 



Consider the ladder-rainbow BSE for a bound state of a quark of flavour /i and an antiquark of flavour 
/2 in the form 

r(p; P) = -\j - ?)7,5^n? + ^^)r(g; P)Sf\q - ^Ph, , (4.40) 

where ^ = 1 — and for simplicity a Feynman-like gauge has again been used. Here the quark and 
antiquark momenta are pi = q + C,P and p2 = —q + ^P respectively while the relative momentum 
that has been introduced is g = ^pi — ^p2- The conjugate coordinates to P, q are X = ^xi + and 
X = xi — X2- If we were dealing with the exact Bethe-Salpeter amplitude and equation, the definition 
of a relative momentum variable is a matter of convenience and ^ is arbitrary (Itzykson and Zuber, 
1980). However the ladder approximation does not respect this arbitrariness and physical observables 
depend on ^. This is an unavoidable deficiency of the ladder approximation and one must impose 
a physical constraint to fix ^. With equal masses mj^ = uif^ as in the case of the pion, charge 
conjugation symmetry requires ^ = 1/2. For unequal masses such as for the kaon where there is no 
charge conjugation symmetry, one can, for example, choose ^ so that the total charge of the neutral 
member is zero. This is discussed further in relation to the kaon charge form factor in Sec. ^71. 



One can write, without loss of generality, 

D{p -k) = Y. D^{p\ e) P - U^{p ■ k) , (4.41) 

n=0 ^ 

where p is the unit-magnitude direction- vector for p and {f/„(x)}^Q are the complete set of orthonor- 
mal Tschebyshev functions. Translational invariance is preserved if all contributing Tschebyshev 
moments are retained. Substitution of Eq. ( [4.41|) into the corresponding DSEs, given in Eq. (p.33|) 
and Eq. (|2.34| ), for quark propagator amplitudes shows that, because the Tschebyshev orthogonality 



weight function matches the angle integration measure that arises naturally in the DSE, only the mo- 
ment Dxij? ^ k"^) enters the DSE for v4(p^) and only the moment -Do(p^, k'^) enters the DSE for B{p'^). 
Hence using the simplest separable representation 

D{p -k) = G{p^) G{k^) + F(/) p ■ k F{e) , (4.42) 

both DSEs, Eq. (|j|) for A and Eq. (p|) for 5, are satisfied if 

G{s) = Ib{s) , F{s) = I {A{s) - 1) , (4.43) 
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where s = and a and h are constants that are determined by 



1 /-oo 

= dtt'[Ait)-l]avit) , (4.44) 

1 f°° 

= iT^ dttB{t)as{t) . (4.45) 
Sir^^ Jo 

This can be done for any quark flavor. Note that the UV behavior of the propagator amphtudes 
in general require some regulation to allow a and b to be determined. Therefore the applications 
are restricted suitably to soft physics. The separable Ansatz BSE kernel is now determined totally 
by the quark propagator and it can be seen that in the chiral limit the Goldstone theorem with a 
massless pseudoscalar is preserved. With given quark propagator amplitudes {as and ay, or A and B) 
constrained by pion and kaon physics, this BSE kernel Ansatz can be used to investigate the extent 
to which the major features of other light mesons are predicted. 



Table 4.4: Selected meson quantities obtained by Cahill and Gunner 
(1995a) from application of a simple form of the separable BSE Ansatz 
fitted to the observables shown. References for experimental values can 
be found in Table 14. 1|. 
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=135 MeV) 


496 MeV (fit) 


496 MeV 
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1230 MeV 




804 MeV 


782 MeV 




0.55 fm 


0.663±0.006 fm 


"0 


0.1634 


0.26± 0.05 


«o 


-0.0466 


-0.028± 0.012 


a\ 


0.0358 


0.038± 0.002 


al 


0.0017 


0.0017± 0.0003 


al 


-0.0005 





Consider the case of the pion BSE with given u and d quark propagators taken as being identical. 
Then use of Eq. ([4.42|) in Eq. (|4.40| ) with ^ = 1/2 for the BS amplitude allows the representation 
r^(p, P) = GijP') Xn{P) where A^(P) is a sum of a few amplitudes multiplied by matrix covariants. 
The integration in Eq. (|4.40| ) now involves only functions known from the propagator model. The 
distinct amplitudes in A^(P) can be separated by trace techniques. In general for mesons of other 
transformation character both separable form factors G{p'^) and Fip^) are involved. The result is that 
the BSE has been reduced to a low rank matrix eigenvalue problem of the form K T = (a(P^) — l)r, 
with the bound state mass identified from a;(P^ = — M^) = 0. The rank of the matrix is the number 
of distinct covariants that this ansatz allows to be coupled by the BSE. For the pion, the rank is 2 
corresponding to both 75 and 757 ■ P covariants. 

In the application of this framework by Cahill and Gunner (1995a), the chiral limit dressed quark 
propagator, in a certain parameterized entire function form, was used to determine the separable 
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form factors. Because the pion mass is one of the constraints that must be satisfied, it is necessary 
to have propagators that correspond to finite bare mass for the quarks. This was obtained by use 
of the separable kernel in the DSE to recalculate the quark such propagators which were then re- 
parameterized back into the entire function form. The simplified BSE was solved under the assumption 
that only the dominant or canonical matrix covariant is important. Some of the results obtained that 
way are displayed in Table Many other quantities have been calculated from this approach by 
Cahill and Gunner (1995a) including mass estimates for several types of di quark constituents for 
baryons, a nucleon mass without meson dressing. 



Table 4.5: Selected meson quantities obtained by Burden et al. (1997) 
from application of a simple form of the separable BSE Ansatz fitted to 
the observables shown. Top panel pertains to meson weak decay con- 
stants; bottom panel pertains to meson masses. The column labeled 
"Dom" refers to the results obtained using only the single dominant co- 
variant; for the other column all contributing covariants have been used. 
References for experimental values can be found in Table WX. 
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In the application of this separable Ansatz by Burden et al. (1997), a UV regulated form of the 
confining parameterized quark propagators described in Sec. |^ was used with finite quark bare masses 
to construct the kernel form factors. This procedure also preserves the Goldstone theorem. With the 
four quantities m^/^^ and f-j^/K fit, the predictions for other meson masses are shown in Table |0. 



Mass estimates for a variety of diquark constituents of baryons were also considered in that work. 
A very good description of the ground-state, SUf{3), isovector-pseudoscalar, vector and axial- vector 
meson spectrum was obtained. 



The result in Table for the rj meson corresponds to the BSE projected onto the SU{3) mixed 



48 



octet-singlet flavor covariant 

Fr, = ^ (A^ cos Op - A° sin Op) , (4.46) 
V z 



with A° = Y'2/3diag(l,l,l) and Op the mixing angle. The exact kernel of the BSE would lead 
to a prediction for 6 p. The rj' meson corresponds to the orthogonal projection produced by — >■ 
6p — 7r/2. One expects quark annihilation to timelike gluons, forbidden in the flavor octet channel, 
to be important in this mainly singlet channel. The ladder kernel can only provide coupling of 
flavor components such as uu and ss through the symmetry-breaking masses. This is insufficient to 
reproduce the large mass increase required by M^f^* = 958 MeV compared to M^^p* = 547 MeV. 
The ladder approximation is found to be adequate for the r] when the mixing angle is treated as an 
external parameter on which the mass and other properties depend. A small positive value for the 
mixing angle is favored by the ladder approximation here while the experimental value is 6p = —10° 
from Particle Data Group (1996) and has significant uncertainties. 

Scalar mesons were poorly described and this is to be expected in ladder approximation which is 
missing important contributions including quark annihilation to multi-gluon configurations and open 
channels like vrvr and K K which can provide a width and a significant change to the real mass. The 
scalars, and also the charge parity exotics (C = —1), are displayed in Table H]^ to point out these 



shortcomings. Within the limited reliability and applicability of the rainbow-ladder approximation to 
the quark-DSE/meson-BSE complex, these results imply that this crude separable Ansatz can convey 
to certain other meson channels, via the dressed quark propagator, the dominant ladder physics 
constrained by the Goldstone mechanism in the psuedoscalar channel. 

It was found by Burden et al. (1997) that in the isovector-pseudoscalar meson channel the sub- leading 
Dirac components of the BS amplitude provide quantitatively important contributions. In particular, 
the pion solution was obtained in the form 

r^(p; P) = Guip'') [175 Al - 75 7 ■ ^' A2] , (4.47) 

where is the time-like unit-vector associated with P^, i.e. [P^ = —1]. This is not the most general 
form of the pion amplitude (see Eq. (|4.3| ) given earlier), but it is the most general form allowed by 



the separable Ansatz. The amplitudes are displayed in Table ^l6| . The sub-leading (pseudovector) 
amplitude is more significant than it appears here since, as is evident from Table [4.5| , it contributes 
~ 15% to the mass and ~ 35% to the weak decay constant /t^. For the kaon, a similar situation is 
found to hold and there is an indication that for the rj — r]' complex the sub-leading covariants are also 
important (along with the processes missing from the ladder approximation as mentioned above). As 
is evident from Table the sub-leading covariants are found to be unimportant for mesons heavier 
than the pseudoscalar octet. 

In the vector meson case, the solutions of the BSE in Eq. ( |4.40| ) are labelled by a Lorentz index, and 
satisfy the transverse condition P^r^{p] P) = 0. Due to the Feynman-like gauge employed, the general 
form is a linear combination of the Euclidean covariants: 

Pu, ll, pIi ■ P, Pul ■ Pp-P, -i-,e^uXplp,P\Pp , (4.48) 

weighted by invariant amplitudes J-'i{p'^ , P"^ , p ■ P). The superscript T denotes a vector transverse 
to P^. The separable Ansatz Eq. ( [4.42| ) introduces further simplification since it cannot support 
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contributions to Tiy{p; P) that are bilinear in p. Hence, at the mass-shell, the produced vector meson 
amplitude is 

r^(p; P) = vlFuip'yXi + illGu{v')\2 + il.e^uXpl^.PxPpF^ip'yX^ , (4.49) 

where the Aj are produced as an eigenvector of a 3 x 3 matrix. For the axial- vector meson the transverse 
Euclidean covariants are 

pIi^P ■ P, 757J, P^757 ■ P, pIi^I ■ Pp-P, e^^xplf,P\Pp ■ (4.50) 

Again the terms bilinear in p do not contribute and, using the separable Ansatz, the axial-vector BSE 
reduces to a 2 x 2 matrix eigenvalue problem. The produced amplitude is 



^LiP^P) = h^llGu^p^YXi +iepyXplp,P\PpFu{p^)X2 
The solution amplitudes are displayed in Table ^]6. 



(4.51) 



Table 4.6: Selected Bethe-Salpeter amphtude solutions from the sepa- 
rable Ansatz of Eq. ( [4.42| ). The amplitudes are physically normalized 





jPC 


Ai A2 A3 




0-+ 


0.61 -0.045 


p/uj 




0.075 -0.33 0.049 




1++ 


0.056 -0.28 0.0 



These studies also show that separable Ansatze have a number of shortcomings. In the pseudoscalar 
channel one finds that the 75 and 757 ■ P components of the meson BS amplitude are characterised by 
the same function, B{p'^), which is not true in general. One also finds that the dominant components 
in the BS amplitudes of the vector and axial-vector mesons are characterised by the same functions 
that characterize these components of the pseudoscalar mesons, BijP'). More realistic studies indicate 
that the vector meson amplitudes are narrower in momentum space. The separable amplitudes should 
be used with caution in the calculation of meson interactions, and should eventually be replaced by 
more realistic calculations. 

The shortcomings notwithstanding, there are areas of study in hadronic physics for which the BS 
amplitudes provided by a separable Ansatz have significantly greater dynamical justification than the 
phenomenological amplitudes currently used. For example, as discussed in Sec. ||, hadronic coupling 
constants such as QpT,.,, and g.yT,p have been reproduced from the qq structure of the mesons in terms 
of a single dominant Dirac covariant if the amplitude is allowed some phenomenological freedom. At 
present, a more realistic treatment is facilitated by use of the separable Ansatz. 

Studies such as these, in combination with analysis (Bender et ai, 1996) of processes in the BSE 
kernel that are of higher order than the bare ladder term, provide a better understanding of the 
domain of applicability of the ladder BSE. The ladder kernel has the defect that it is purely attractive 
in both the color-singlet q-q and color-antitriplet q-q channels. This entails that it yields bound color- 
antitriplet diquarks. This is a peculiarity of ladder approximation. Measuring "order" by the number 
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of dressed-gluon lines in the Bethe-Salpeter kernel, ladder approximation is the lowest order kernel. 
Repulsive terms appear at every higher order. It has been shown by Bender et al. (1996) that in the 
isovector-pseudoscalar and vector meson channels, these repulsive terms are cancelled by attractive 
terms of the same order. This explains why ladder approximation is phenomenologically successful in 
these channels. That study also showed that in the color-antitriplet diquark channel the algebra of 
SUc{^) entails that the repulsive terms are stronger; they are not completely cancelled and eliminate 
the diquark bound states. 



5 Electromagnetic Coupling 

5.1 From Quark Currents to Composite Meson Currents 

To allow for electromagnetic (EM) coupling to the meson modes, the bosonization method can be 
generalized (Frank and Tandy, 1994) to account for a background electromagnetic field minimally 
coupled {di, di, — iQA^) to the bare quarks. Thus the action 

S[q, q,A]=J d'xq{x){^ - iQJ/. + m)q{x) + ^ | d^xd'y3l{x)D{x - y)3l{y), (5.1) 

where Q is the quark charge operator, is invariant under the gauge transformation 

q{x)^q'{x) = e^'3^(^)g(x), 
q{x)^q'{x) = q{x)e-'^'^^''\ 
A,{x)^A'^{x) = A,{x) + d,e{x). (5.2) 

The standard bosonization technique can be apphed in the presence of the external EM field to the 
expose qq meson fields and their electromagnetic interactions. Gauge invariance of the EM coupling 
at the quark level can be translated to the level of extended meson modes (Frank and Tandy, 1994). 
One obtains Z[Al\ — N J DE exp(— »S[;B; ^4]) where the gauge invariant boson action is 

S[B; A] = -TrLn g-' [B; A] + J d'xd'y^^^^^, (5.3) 

with = ^ — iQj^ + m + A'^B^. The saddle point configuration of the auxilary fields B'^, 

defined in the same manner as previously, now becomes a nonlinear functional iBo[A] of the background 
EM field. It produces a ladder dressing of the photon-quark vertices as well as the translationally 
invariant quark self-energy. It is convenient to define = A'^Bq [A] which plays the role of the 
quark self-energy in the background EM field. It satisfies 

4 

S(a;, y) = -D{x - y)-fugix, y)-^y. (5.4) 

where 

g-\A\^ ^-iQ^ + m + t[A\. (5.5) 

The structure of this formalism is exactly the same as when A^ = Q except that now the generalized 
quantities g and are not translationally invariant. The first-order dependence of g~^ upon Ai, 
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generates the EM vertex with dressed quarks via 



(27r)^ S{p-k-Q)T, 



p + k 



Apphcation to Eq. ( |5.4| ) immediately yields the ladder BSE for this vertex 



r.(g; Q) 



d'^k 
(2^ 



D{q-k)^,S{k+)T,{k-Q)S{k.)^^ 



(5.6) 



(5.7) 



where S = Q[Ai, = 0] and k± = k ± ^. Use of the rainbow SDE for the quark propagator S also 
confirms that this vertex satisfies the Ward-Takahashi identity (WTI) (Takahashi, 1957) 



Q.T,{q;Q)=Q{S-\q^)-S-\q^)}, 
and also the Ward identity (Ward, 1950) 

r.(g;g = 0) = 



Q ^S-\q). 
dqu 



(5.8) 



(5.9) 



These are properties the ladder vertex shares with the exact vertex. For this reason the ladder- 
rainbow pair of equations for the dressed quark propagator and EM vertex constituent a very useful 
approximation in that EM current conservation can be exactly preserved. If instead of the EM vertex, 
we the consider the ladder BSE for the Goldstone bosons, it was pointed out in Sec. ^]l] that the 
ladder-rainbow combination preserved chiral symmetry and the Goldstone theorem. It is therefore 
interesting and advantageous that, at meson tree level of the bosonized GCM, the produced quark 
substructure has exactly the ladder-rainbow content. 



The action may be expanded about the saddle point configuration Bo[A] to obtain 

S[B;A]=S[Bo[A];A]+S[B;A], 



(5.10) 



where B'^ = B'^ — Bq[A] are the new field variables for the propagating meson modes and the corre- 
sponding action for the meson sector containing electromagnetic interactions is 



S[B;A]=TtY: 



i-iy 



n=2 



n 



g[A] A^B^ 



-li 



(5.11) 



The EM coupling to the qq meson modes is contained in the first term of Eq. ( |5.11|) and can be 
unfolded through use of Q[A] = S — Ti,Ai,S + ■ ■ ■. The meson EM current is 



6S[B^; A] 



6AJx) 



6S[B'^;A] 



A„=0 



SA^{x) 



(5.12) 



A„=0 



The second equality in Eq. ( p.l2| ) expresses the fact that the saddle point action S [Bo[A]; A] does not 
contribute to the current. The functional dependence of S [Bo[A]; A] upon A^ that enters implicitly 
through Bo [A] makes no contribution since ^ = 0. The explicit A„ dependence generates a bare 
coupling of the photon to a vacuum quark loop which vanishes due to symmetric integration. The 
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saddle point action can contribute to higher order in A^, beginning with a vacuum polarization insertion 
UJ^KQ) for the photon propagator. In particular, the part of iS[i3o[A]; A] which is second order in A^, 
is 



where 

More explicitly this is 



5(2) [So [A]; A] = i / ^,M-Q) ^IZiQ) MQ), (5.13) 

UI:{Q) = -Tt[z^,QST^{Q)S]. (5.14) 

^IZiQ) = / M Siq+) r.(g; Q) S(g„)] . (5.15) 

This has the same form as the exact polarization tensor although the dressed vertex here is in ladder 
approximation as is appropriate for a treatment at zeroth level of the meson loop expansion. Due to 
the WTI in Eq. (|^), it follows that Q^ITJ^iQ) = producing the transversality required by gauge 
invariance. To go further in this direction and consider the renormalization of the photon propagator 
and field in conjunction with the required counter terms, the photon needs to be treated as a dynamical 
field variable instead of an external field. An important element is that the transverse amplitude of 
UJKQ) has a zero at = 0, that is Ilfj_^{Q) = {Q^S^^, — Qij,Qu)^{Q'^), thus allowing the physical 
renormalized photon to remain massless. This property applies to Eq. ( |5.14| ) as long as satisfies 
the Ward identity (Burden et ai, 1992a). 

One of the pleasing features of the action for electromagnetic coupling to interacting mesons that 
follows from Eq. (|5.1CI|) and Eq. ( ^.ll]) is that the current conservation required by gauge invariance is 



manifest in each term. One can then continue on and derive photon couplings to mesons in a way that 
maintains meson EM current conservation in the face of meson substructure or form factors that are 
accountable to the conserved currents at the quark level. This has been done for finite size pions by 
Frank and Tandy, (1994) where pion Ward-Takahashi identities are also derived. There it is also shown 
that when the pion action is restructured to an equivalent form governing effective local pion fields, 
then EM gauge invariance in the presence of the extended charge form factor will only be manifest if 
one recognizes the gauge transformation behavior of the effective pion inverse propagator that results 
from the substructure field content that has necessarily been absorbed into it. Such observations can 
greatly assist in the task of gauging effective lagrangians for hadrons that involve form factors related 
to substructure. 

An inconvenience of this approach is that there is no explicit term in Eq. (|5.11|) that is linear in a 
meson field such as the 7r° — > 77 process, for example. Also the physical processes involving the direct 
coupling of a photon to a vector meson (as needed for the decay p° — > e'^e~) are not evident. Such 
processes form the basis of the empirically successful Vector Meson Dominance (VMD) hypothesis 
and it is desirable to unfold them from Eq. (|5.11|) . In the point coupling limit, e.g. the NJL model. 



it is known that a point vector meson field redefinition via a simple shift p'^ = py + cAy reveals the 
VMD terms (Volkov and Ebert, 1982; Volkov, 1984; Ebert and Reinhardt, 1986) and in fact allows 
the vector meson fields to be chosen so that they mediate all electromagnetic coupling to hadrons. 
This idealization entails a mass-type (i.e. constant amplitude) mixing of the photon with point vector 
mesons, and is often taken to be the essential meaning of VMD. In this situation, an exact current-field 
identity (Sakurai, 1969) exists; that is, the EM current is a linear combination of the point vector 
meson fields. 



53 



Since the EM current is a local object, and in reality the vector mesons (both from the present 
GCM model and in QCD) are distributed or nonlocal objects, the VMD current-field identity can 
only hold approximately. Only the soft photon coupling to a nonperturbatively interacting qq vector 
correlation can be implemented through realistic vector meson degrees of freedom; there is a point g7g 
coupling {—i'jfj.Q) that cannot be properly simulated by a bound state object. Only in the context 
of effective lagrangian models with point-like vector meson fields can such fields take over the total 
role of the photon. For quark models (and QCD) where realistic finite size meson modes arise, a 
meson field redefinition that attempts to absorb the point q'-fq coupling does not make sense and leads 
to dynamical problems. However, one still expects to be able to implement resonant photon-hadron 
coupling via the composite vector meson fields leaving a non-resonant direct photon coupling to quarks 
while maintaining manifest gauge invariance. We will indicate how this can be done. 




Figure 5.1: The separated non-resonant and vector meson pole compo- 
nents of the dressed photon-quark vertex. 



The quantity S[y4] that enters Eq. (|5.5| ) contains ladder dynamics for the coupling of one or more 
photons with dressed quarks. For example, the 1-photon part must contain jp and 70; processes 
and the 2-photon part must contain a 71*^77 process. The former are contained within the dressed 
quark-photon vertex T,^ that was produced. This becomes evident for photon momenta near the 
vector meson mass shell where the homogeneous version of Eq. ( |5.7|) has a bound state solution. This 
generates a bound state pole in the inhomogeneous solution for the vertex F^. Thus, with respect to 
the invariant mass of a given vector meson, T^{q; Q) can be separated (non-uniquely) into a resonant 
or pole piece (which is transverse) and a background or non-resonant piece (which is both longitudinal 
and transverse). That is, = F"*^ + F^"'*^ and we shall consider the p pole to be the point of interest. 
We may formally make this separation using the eigenmode analysis of the kernel of the ladder BSE 
vertex Eq. ( [5.7| ) as discussed in Sec. |2^ . With the dressed part of FJJ'" representing the contributions 
from the eigenmodes other than the p, this yields 

r.(g; Q) = Fr (g; Q) - r^(g; Q) q2^^"^Iq2^ ^'^HQ)- (5-16) 
The p7 polarization tensor that is generated is given by 

IIZ{Q) = -Tr ni-Q) S il.Q S]. (5.17) 



This pole separation of the photon-quark vertex is illustrated in Fig. |5.1j. The employed p BSE 



amplitude is physically normalized and the associated dynamical mass function mp{Q ) is defined in 
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Sec. pr^ . Although a ladder BS structure is under discussion here, the exact IPI photon vertex can 
also be decomposed in this form. There will in general be an imaginary part to the mass function 
arising from physically accessible hadron channels (such as p ^ tttt) in a more complete analysis. The 
point of the present discussion does not require that such aspects be made explicit. The off-mass-shell 
behavior of individual terms in separations such as in Eq. (|5.16|) will differ in different approaches and 



adjustments can be made to accommodate convenient approximations for F"*^. 

With the non-resonant F"'' satisfying both the WTI and the Ward identity, then the above pole 
term of the vertex is zero in the soft photon limit. That is, with Il'j^iQ) = T^^{Q)Ilip {Q"^) , we have 
11^7(0) = 0. Thus, due to manifest EM current conservation, the photon decouples from the composite 
vector mesons at = 0. Similarly, with the separation in Eq. ( |5.16| ), the qq self-energy or polarization 
tensor for the photon, from Eq. (|5.14|) , can be written 



UlUQ) = -Tr M S Fr (Q) S] - UjfsiQ) '"ZlL, ^ZiQ)- (5.18) 



Since the first term can be shown to be zero in the soft photon limit if F"*^ satisfies the WTI and 
the Ward identities (Burden et ai, 1992a), the result 112^(0) = is again evident and linked to the 
massless condition of the photon. In general there can be many pole contributions to the spectral 
strength in the timelike domain that can be separated out. The common assumption is that selected 
vector mesons, such as p, u and 0, saturate the spectral strength. 

Although the pole term in the photon vertex admits an integral representation in terms of vector 
meson fields, there are already explicit vector meson fields elsewhere in the effective meson action of 
Eq. (|5.11|) . It is an awkward matter to consistently unfold photon-meson mixing from the previous 



formalism this way. It is more convenient to make the following field redefinition. We expand the 
bosonized action Eq. ( p.3|) about a new configuration ;B^[y4] of the auxiliary fields, that is not the 
saddle point Bq [A] as used before, but is chosen so that the resulting photon-quark vertex is the non- 
resonant dressed vertex F"*". Specifically, the new auxiliary field configuration yB^fA] about which we 
expand is defined so that the inverse quark propagator in the absence of the meson fluctuation fields 
is ^ 

S-'[A] = ^-iQ^ + m + A^B^[A] = ^ + m + S + F^'^A,. (5.19) 

Comparison with Eq. ( pl5|) indicates that the previous quantity E[^] is now replaced by a quantity that 
is linear in the photon field A^, is devoid of bound state poles, but can still allow current conservation. 
A practical Ansatz for F"'' is considered shortly. 

To proceed with the meson field redefinition, the bosonized action Eq. ( p.3| ) is expanded about B^[A] 
to produce 

S[B; A] = S[B[A];A] + S[B; A], (5.20) 

where B^ = B^ — B^ [A] are the new field variables for the propagating meson modes (although we 
use the same notation B'^ as before). The new action for the meson sector with electromagnetic 
interactions is 

S[B; A] = Si[B] A] + Tr ^ ^ — (^S[A] A^B^) + i (^B^*D~^B^) , (5.21) 

n=2 ^ 

where the term <Si [B; A] is first order in meson fields and is present because the expansion is not about 
an extremum. The second term is at least quadratic in mesons and contains conventional meson charge 
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form factors processes such as 77r7r. The processes involving one or more photons coupled to a single 
meson (including the so-called VMD processes) are exposed by the first term and its explicit form is 



Si[B;A] = -Tr [S[A] A'^B'''] + (B'''[A]*D-^B' 



(5.22) 



The 0{A^) terms of Si cancel out, i.e. <Si[B; A = 0] = 0, because then the expansion point is the 
A = saddle point configuration and the cancelling two terms in Eq. ( |5.22|) form the rainbow DSE. 
The 0{A^) terms of Si provide the sought-after photon-meson mixing terms. The (9(y4"') terms of 
Si with n > 2, such as 7r°77, are contained only in the quark loop term and are revealed by the 
expansion S[A] = S — SV^^Ai^S +ST'^^Ai,ST]l;^Ai,S + ■ ■ ■. We now use the composite meson field 
representation of Sec. to introduce the relevant BS amplitudes and associated local field variables, 
A't'B'f' r„(g;Q)0„(Q). The 0{A^) term of Eq. (§]^ then becomes closely related to the ladder 
BSE kernel, allowing the term to be recast as a quark loop by using the methods of Sec. p^ . 
Consideration of the relation between T^^ and the eigenmodes of the ladder BSE kernel indicates 



partial cancellation between the two terms of Eq. ( 5.22 ) due to orthogonality. The results can be 
illustrated by the explicit jp mixing term which is 



A^{-Q) UZiQ) Pu{Q) 



where the mixed self-energy tensor that is produced is 



Ul^^{Q) = -^T[^J,QST^,{Q) S] 



(5.23) 



(5.24) 



Here, at the p mass-shell, r^(g; Q) is the physically normalized BS amplitude, and for general Q it is the 
relevant eigenvector of the ladder BSE kernel defined in Sec. p.2| . The result in Eqs. (|5.23|) and (|5.24|) 
is consistent with the original analysis in Eqs. ( |5.16| ) and ( p.l7| ) and that p pole contribution to each 
photon-meson coupling would be recovered through integrating over the p field of the reformulated 
action. 



We summarize by retaining only a few low mass meson modes so that the reformulated meson action 
of Eq. ( p. 211 ) can be written in the more explicit form 



S[7l,p, 



Tr 



S[A] (f . • vf + ■ + ■ ■ Ol + ^ E / 7^ C(^) KLm;P) MP) 



n=3 



n 



s[A] (r^.vf + r^-p^ 



(5.25) 



where the local fields are (pn = vr*, p^, ■ ■ -. When Ay 0, this action reduces to Eq. (|2.64| ) that 
was obtained earlier. Here the generalized inverse propagator A~4([A]; P) is defined in an obvious way 
from the quadratic terms of Eq. (|5.21|) . We see in Eq. (|5.25|) that the photon coupling to hadronic 
matter consists of a resonant coupling mediated by single vector mesons which is part of the first 
term, as well as a direct non-resonant coupling to the conserved EM current of dressed quarks and 
represented by the other terms. A by-product of this reformulation is that the processes involving 
multiple photon coupling to a single meson are uncovered and are found in the first term. The latter 
include the chiral anomalies such as '7r''77 and 7777. This formulation is most useful when the F"'' 
satisfies the Ward-Takahashi (Takahashi, 1957) identity and the Ward identity (Ward, 1950); then 
EM current conservation is manifest for each term of the action in Eq. ( |5.25| ). 
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The solution of the DSE for the dressed quark-photon vertex is a difficuh problem. Only recently 
has the solution of the ladder version given in Eq. (^.Tp begun to be addressed (M.R. Frank, 1995). 
However, much progress has been made in constraining the form of r^(p; Q) and developing a realistic 
Ansatz (Ball and Chiu, 1980; Curtis and Pennington, 1992; Burden and Roberts, 1993; Dong, Munczek 
and Roberts, 1994). The bare vertex r^(p; Q) —>■ —i'^^Q is clearly an inadequate approximation when 
the quark propagator has momentum dependent dressing because it violates the Ward-Takahashi 
identity would lead to an EM current that is not conserved. An Ansatz with the desirable properties: 
a) satisfies the Ward-Takahashi identity; b) is free of kinematic singularities; c) reduces to the bare 
vertex in the free field limit as prescribed by perturbation theory; and d) has the same transformation 
properties as the bare vertex under charge conjugation and Lorentz transformations, is (Curtis and 
Pennington, 1992) 

T,{v; Q) = Q (ri^(p; Q) + rj(p; Q)) (5.26) 

where, with S^^{p) = 27 • pAlp^) + B{p'^) + m, the term advocated by Ball and Chiu (1980) is 

T';{p; Q) = -Z74(^(P+) + MP-)) + ■ p{MP-) - Mp+)) + {B{p_) - B{p+))] . (5.27) 

P ' ^ 

with p± = p ± Note that the undetermined piece P^ is transverse, P)^(5'; Q) = and vanishes at 
zero momentum Pj((3'; 0) = 0. It is also the case that with a bare quark propagator, which has A = 1 
and B = constant, P^ = 0. In this Ansatz the T'^^ piece is completely determined by the dressed quark 
propagator and it has both transverse and longitudinal components. A feature of this approach which 
follows from criterion c), is that with a quark propagator that incorporates asymptotic freedom the 
quark-photon vertex will reduce to the correct term at large spacelike-Q^, in the manner prescribed by 
perturbation theory in QCD. This Ansatz therefore provides a realistically constrained extrapolation 
of the quark-photon vertex to small spacelike-Q^. 



The representation for P^*" that has been used in most investigations of electromagnetic interactions of 
hadrons within the GCM or DSE format is the above Ball-Chiu Ansatz, that is, T^^{p; Q) ^ Qr^^{p; Q). 
The additional resonant or vector meson pole processes, that we have earlier identified within the 
GCM, provide a specific dynamical model representation for the P^ term in Eq. ( p.26| ). The reso- 
nant terms are obviously necessary in the timelike region sufficiently close to the poles; however in 
the spacelike region, the relative contribution of P"'" and pp°'^ to specific hadronic processes is very 
much dependent upon the off-mass-shell behavior of the particular separation employed. In another 
language, while the individual contributions depend upon choice of interpolating field variables, the 
summed contribution to an S-matrix element should not. Within the GCM and DSE approaches, it 
has been found that the spacelike pion charge form factor can be described quite successfully (obtain- 
ing 85 — 90% of the charge radius) through use of the non-resonant T^^ Ansatz. We shall comment 
further on this point later. There are several other choices available to model the non-resonant EM 
vertex (Roberts and Williams, 1994). However no inadequacy of the Ball-Chiu Ansatz has been found 
in practical calculations of form factors at spacelike momenta. 

In a recent analysis of effective local lagrangian models for electromagnetic coupling to point fields 
for vector mesons and other hadrons, O'Connell et al, (1995) have pointed out the merits of a form 
in which the VMD hypothesis is implemented in precisely the "direct plus resonant" manner arrived 
at above. A property of that particular local representation of VMD (called VMDl), is that the 
photon- meson mixing amplitude is momentum-dependent and vanishes in the soft photon limit. A 
point field redefinition via, e.g. p'y = Pv + cA^ can produce the equivalent, and more common form 
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(VMD2) of local effective lagrangian where the photon only couples to vector mesons and the relevant 
amplitude is a constant. This form is inconvenient as the zero value of the physical photon mass is 
only recovered as a cancellation between finite bare and dynamically produced contributions; also the 
reproduction of the total charge of hadrons requires fine tuning of parameters which is indicative of a 
formalism wherein current conservation is not manifest. 

The vanishing at zero momentum for the transverse vector self-energy amplitudes we have discussed is 
related to, but somewhat more general than, the node theorem for the self-energy amplitudes for point 
vector bosons coupled to conserved local currents that has been put forward recently (O'Connell et ai, 
1994). In the present analysis, it is the participation of the photon that precipitates the result, for it 
alone couples to a local conserved current. There is no aspect of the self-energies for finite size vector 
mesons that are involved in the dressed meson propagator of Eq. (|5.18|) that would produce a zero 
momentum node. In a similar fashion, with flavor asymmetric quark bare masses, the p — uj mixed 
self-energy amplitude contained in the coupled channel meson propagator that would now appear in 
Eq. ( |5.18| ) would not have a zero-momentum node. This is because the meson modes here, and in 
QCD, are extended objects. One would expect however, that a nearby node might occur if the meson 
size was suitably small. Indeed, a calculation of the p — oj amplitude at the quark loop level based on 
the GCM produces a node at about 0.1 GeV^ in a way related to meson size (Mitchell and Tandy, 
1997). The only model based at the quark level where the zero momentum node property would extend 
to the vector meson self-energies is evidently the NJL model which allows an exact VMD current-field 
identity to be realized. This follows from the observation that, in this point coupling approximation, 
the vector meson fields are proportional to the various flavor components of the photon; the relevant 
meson currents are local and conserved and the analysis reduces to the electromagnetic one. 

A recent study within a generalized NJL model utilizes just this correspondence between the vector 
mesons of the model and the flavor components of the photon to analyze VMD, p — uj mixing and 
the pion charge form factor from a quark basis organized around local current correlation functions 
(Shakin and Sun, 1997). In that work the NJL model is generalized in the sense that a regulated 
linear confining potential is added to dress quark vertices to eliminate spurious production thresholds 
from qq polarizations. 



5.2 Pion and Kaon Electromagnetic Form Factors 

The action component describing the '-jtxti interaction, in the form —J-A, can be written as 



S [77r7r] 



d'P,Q 
(27r)8 



Q 



n*_{P + ^)A,{Q)A,{P;Q)n4P-^) 



(5.28) 



which emphasizes the coupling of a photon with momentum Q to the electromagnetic current of a vr"*" 
with initial momentum P + ^ described by the field vr. = (tti — ^2) / a/2. An equivalent expression is 

d^P,Q 



S [77r7r] 



A,{Q) 3 ■ n{-P - |) X n{P - |) A,(P; Q). 



(5.29) 



The vertex A^ = (en)A^ + (— ed)A^ consists of contributions from the photon coupling to the quark 
and antiquark. With u and d quarks taken to be identical except for their electric charges e„ and e^, 
it is easily shown that AJ^ = A^ = A^. The second term of the action produced in Eq. ( |5.21| ) yields 



5 [7 



TTTT 



-Tr 



S T^A^ {S 7f ■ r.) 



(5.30) 
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One then obtains the 77r7r vertex at meson tree level in the GCM as the expression 



A^(P;g)= (5.31) 
2iVctr,| ^r.(g + |;-P+)^(g_, + |)r,(g+;Q)^(g+-|)r.(g-|;P_)5(g_) , 

where we have used P± = P ± ^ for the final and initial momenta of the pion and also q± = q± ^. 



The momentum assignments used for this vertex are illustrated in Fig. |5.2| . The trace over color 
and fiavor has been carried out leaving only the trace over Dirac spin. The quark charge operator 
Q has contibuted to the fiavor trace relations tr/r+Qr. = 2eu and tr/r_QT+ = 2ed. In Eq. (|5.31|) 
r^(g; P) denotes the pion Bethe-Salpeter amplitude in Dirac spin space, r^(g; P) = [C~^r,r(— ^Z; P)C]^ 
is for the corresponding antiparticle and here reduces to r^(g; P). The dressed quark-photon vertex 
associated with in and out quark momenta g ib ^ is denoted r^(g; Q). The above result is also what is 
obtained from the DSE approach to QCD truncated at the generalized impulse approximation. The 
generalization refers to the propagators and vertices being dressed. The dressing of each element is 
consistent in a way demanded by chiral symmetry and EM current conservation. 




Figure 5.2: Diagram for the generalized impulse approximation to the 
77r7r vertex. 



In the specific study (Roberts, 1996) that we summarize here, the Ball-Chiu Ansatz r^(g; Q) ^j^iq', Q) 



is used thus making the generalised impulse approximation given in Eq. ( p. 31 ) completely determined 



by the quark propagator without the need for new parameters. The pion charge form factor data 



was used, along with the other pion properties mentioned in Sec. ^?1| , to set the propagator parame- 
ters. (Various predictions from that parameterization are summarized in Sec. ^.) This approach thus 
entails that vector meson contributions are explicitly excluded and must be treated as an additional 
contribution. 

Using the charge conjugation properties: S{—k)'^= C^S{k)C; T^{—q;P)= C'''r^(g;P)C; and 
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r^(— g; Q) = —C"^ r^(g; Q) C, where C = 7274, one finds that the 'jim vertex has the symmetry prop- 
erties A^(P; Q) = — A^(— P; Q) = A^(P; —Q). Hence the most general form is 

A^(P; Q) = 2P^ P^(P2, Q^, (P . QY) + 2Q^P ■ Q H^{P\ Q\ (P • Qf). (5.32) 

For elastic scattering, (P— = (P+ ^)^, that is P • Q = 0. Thus the conservation of the electro- 
magnetic vr-current 

Q^A^(P;Q) = 0, (5.33) 

follows automatically in this formulation independently of the details of the parameterization of the 
dressed quark propagator. The physical form of the vertex is 

A^{P-Q) = 2P^F^{P\Q^) , (5.34) 
where the pion mass-shell condition is P^ = —rn^ — 

At = 0, the condition that the correct physical charge is obtained is equivalent to the physical 
normalization condition for the pion BS amplitude. To see this, the Ward identity for the photon- 
quark vertex r^(g; 0) = —dS^^{q)/dq^ can be used in Eq. ( ^.31|) to produce 



A^(P;0) = 2P^P^(P^0)= (5.35) 
d r d'^q — 



Comparison with the ladder BS normalization condition in Eq. ( p.61| ) shows that the former is 
equivalent to P7r(P^,0) = 1. In other words, generalised impulse approximation combined with a 
P-independent Bethe-Salpeter kernel provides a consistent approximation scheme. Also the condition 
of physical net electric charge may be viewed as the normalization condition for BS amplitudes. The 
expression for given in Eq. ( [4.22| ) can be obtained directly from the chiral limit of Eq. ( |5.35|) 



Since the quark propagator in use here is an entire function, the loop integral in Eq. (|5.31|) is free of 
"endpoint" and "pinch" singularities. This is a particular, sufficient manner in which to realise the 
requirement that Fig. |5.2| have no free-quark production thresholds. It also facilitates the numerical 
integration. 



The results obtained for Ft^{Q'^) are shown in Fig. ^]3| and Fig. |5]^ taken from Roberts (1996). The 



parameterization employed for the quark propagator in that work differs slightly from the later form 



that we have summarized in Sec. y. However Fig. |5.3| and Fig. are a faithful representation of 
the charge form factor resulting from the parameterization in Table |3.1| . We recall that the data 
for Ft^{Q'^) are included in the set of pion and kaon observables to which the propagator parameters 



are fit. Other pion and kaon observables arising from the same model are listed in Table |3.2|. Since 



at zero momentum the photon-quark vertex used here satisfies the Ward identity, the pion charge 
radius r^± generated by the generalized impulse approximation is totally determined by the quark 
propagator. The value of fnTn^ preferred by the fit was rather stable and having chosen to fit to the 
experimental /,r, the r^± value shown in Table p.2| , which is some 15% below the experimental value, 
is strongly preferred. Earlier GCM studies (Roberts et al, 1994), with quite different representations 
of the dressed quark propagator, also find essentially a 10 — 15% low value of r^^ from the generalized 
impulse approximation. Interestingly, the additional pion loop contribution to r^^ has been calculated 
within the GCM format (Alkofer, Bender and Roberts, 1995) and was found to have an upper limit 
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Figure 5.3: The pion electromagnetic form factor fit at low q^. Adapted 
from Roberts (1996). The data shown by crosses are from Amendolia et 
al. (1986a), the data shown by circles are from Bebek et al. (1976). 

of that size. The small role for pion loops is consistent with the lattice QCD studies of Leinweber and 
Cohen (1993). We discuss this point in more depth in Sec. |^. 

These results imply that with photon coupling to the dressed quark core of the pion as we have 
described, the p-meson contributes little to F-j^lQ"^) in the spacelike region, even though it is important 
for timelike-Q^ near the p pole. This is especially clear in the NJL model limit of the formulation 
described here where, for example, our dressed photon-quark vertex T^^{q] Q) —iln- The additional 
p pole term has been included in some NJL model studies (Blin, Hiller and Schaden, 1988; Bernard 
and MeiBner, 1988). In Sec. ^we present an analysis of whether such a reduced role for VMD processes 
is consistent with the GCM/DSE implementation of ladder structure for the pion in terms of dressed 
quarks. 

The asymptotic behavior of the generalized impulse approximation to -F^(Q^) was also investigated 
by Roberts (1996). The form of the quark propagator in Eqs. ( |3.23|) and ( p.24|) leads to all elements 
of the quark loop being consistent with the UV behavior of QCD up to logarithmic corrections. In 
particular, the quark propagator generates r^(g;0) oc -B(g^) — >■ ^ which, up to the Inp^ corrections 
associated with the anomalous dimension, reproduces the ultraviolet behavior of the BS amplitude 
in QCD (Miransky, 1990). It also requires, via the Ball-Chiu Ansatz, that the dressed photon-quark 
vertex evolve from soft to hard behavior in accord with such behavior of the quark propagators that 
attach to it as constrained by the Ward-Takahashi identity. From the non-factorized quark loop 
integral in Eq. ( [5. 31] ), the asymptotic behavior of Ft^{Q'^) is found to be determined by that of the 
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Figure 5.4: The pion electromagnetic form factor fit at high . Adapted 
from Roberts (1996). The data shown by crosses are from Brown et 
al. (1973), the data shown by diamonds are from Bebek et al. (1976), the 
data shown by circles are from Bebek et al. (1978). 



bound state BS amplitude in the form (Roberts, 1996) 



oc 



(2^ 



(5.36) 



where Xnil) is the chiral limit BS wavefunction Xtt{(1) = S{q)B{q'^)S{q). In the asymptotic analysis 
leading to this, the photon vertex tends to its bare limit, but the integrand can always find support 
in a domain that leaves several of the loop integral elements evaluated at soft values of their momen- 
tum arguments. Here the mass-shell pion wavefunctions dominate. With a constituent quark mass 
approximation for propagators, one has the approximate behavior Xn{q) ~ 75(q'^ + A^)~^ consistent 



with asymptotic freedom. Thus the integration in Eq. ( 5.36|) , after a now standard evaluation, gives 



FAQ' 



oc 



ing2 



(5.37) 



This is confirmed to be the leading behavior of the numerical evaluation of the full loop integral 
(Roberts, 1996). The \np^ corrections that arise due to the anomalous dimension of the propagator 
and Bethe-Salpeter amplitude in QCD, are found to only lead to the modification of Eq. (|5.37| ) in 
which ln[(5^] — > ln[Q^]'>', where I7I is 0{1). This asymptotic term is found to dominate only for 
> 10 GeV^ where it departs from the perturbative QCD factorization form (1/Q^). The latter 
is based on a factorization of the loop integral to separate out a hard factor leaving a soft integral. 
This requires assumptions about the favored loop kinematics guided by the singularity structure of 
the loop integral appropriate to non-confining quark propagators and is not applicable here. The 
result in Eq. ( |5.37| ) is due to full evaluation of the loop integral without factorization and it is found 
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that both bound state amphtudes contribute asymptotically. The resulting peak in Q'^F.„{Q'^) about 
~ 3 GeV^, which is also found for the form factors, is a signal of non-perturbative quark- 
antiquark recombination into the final state meson in this exclusive elastic scattering process. The 
same phenomena is found to occur with the dressed loop integral for the 7*7r° — > 7 transition form 
factor discussed in Sec. |6.1| . There the dressed vertex for the final real photon replaces the final pion 
amplitude; the real photon vertex remains soft. 

The extension of this approach to the kaon electromagnetic form factors has been studied by Burden, 
Roberts and Thomson (1996). In that work the generalized impulse approximation to the 'yKK vertex 
for the is written 

Af (P; Q) = (e.)A;i(P; Q) + (-e.)A^(P; Q) (5.38) 
where the u quark contributes 

A;:(P;Q)= (5.39) 
2iVc tr. / ^ r^(g + e|; -P+)5"(g+ + ^)T;{q^; Q)5"(g+ - |)r;,(g - e|; P^)S%q^) , 

and the s quark contributes 

A^(P;g)= (5.40) 
2iVctr./ ^r;,(g + e|;-P+)5^(g+ + |)r^(g+;g)5%+-|)r;,(g-^|;P_)5"(g_) . 

The K~ vertex is the negative of this. The momentum assignments are the same as used for the pion in 
Fig. ^.2| . We have used P± = P ± y for the final and initial momenta of the kaon. The conventions are 



the same as for the pion vertex, Eq. ( 5.31| ), except that now the required flavor dependence of the quark 



propagators and the photon-quark vertex are indicated as S-^ and Tj^ respectively, and the parameters 
$, and ^ = 1 — reflect the choice of qq relative momentum variable for the kaon BS amplitude Tk- 
That is, the kaon BSE is of the ladder form given in Eq. ( |4.4CI|) and the amplitude Txik; P) employs 



the choice of relative momentum variable k = ^pi — ^p2- The notation used here for the momentum 
dependence of the quark propagators is S^{q±) = ^"(g ± ^P) and S^{q±) = S^{q ± ^P). As with the 
pion, the physical elastic scattering vertex can be written A^^(P; Q) = 2P^Fk+{Q'^)- 

As discussed in Sec. the underlying ladder BSE that is assumed here does not preserve the 
independence with respect to the partitioning parameter ^ that would be a property of an exact 
treatment. The constraint imposed by Burden, Roberts and Thomson (1996) to fix ^ is to require 
that the neutral kaon have exactly zero charge, i.e. Fko{Q'^ = 0) = 0. This works in the following 
way. The generalized impulse approximation to the vertex is given by 

Af (P; Q) = (e,)A;(P; Q) + (-e.)A^(P; Q) . (5.41) 

The u and d quarks are taken as identical except for the electric charge and hence A^ = A^ = A^. If the 
u and s quark properties, except for the charge, are set equal in Eqs. (|5.39|) and (|5.40|) the generalised 



impulse approximation to the ■ymr vertex, Eq. (|5.31| ), is recovered provided that Tk — > F^r, fx fn 
and ^ = 1/2, which is required by charge conjugation symmetry. This provides a consistency check. 
The same limit applied to Eq. (|5.41|) verifies the charge conjugation symmetry requirement that there 
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Figure 5.5: Kaon charge form factor from Burden, Roberts and Thomson 
(1996). The available data is shown and is from Amendolia et al. (1986b), 
and Dally et al. (1980). 



is no photon coupling to the (self-conjugate) vr^. This formulation ensures EM current conservation 
because the quark-photon vertex satisfies the Ward identity and produces Fk±{Q^ = 0) = 1. This 
happens because at = the generalized impulse approximation for the charged kaon vertex reduces 
to the normalization condition for the BS amplitude in a similar way that led to Eq. ( |5.35|) for the 
pion. However for the non-self-conjugate K^, there is not complete cancellation in Eq. ( |5.41| ) at = 
unless ^ is suitably different from 1/2 to counteract the different dynamics for the d and s quarks. 
The presence of this dependence reflects a deficiency of the ladder approximation to the BSE. 



This formulation was applied by Burden, Roberts and Thomson (1996) using the Ball-Chiu Ansatz 
for r^. The u/d quark propagator parameters were fixed by pion physics and the only free parameters 
were the two s quark parameters C^^^q and 63 described in Sec. |[ These parameters are determined by 
requiring that the model reproduce, as well as possible, the experimental values for the dimensionless 
quantities fx/U = 1-22 ± 0.02, rK±/r^± = 0.88 ± 0.07 and m^//^ = 4.37 ± 0.05. The ansatz 

r^(g; P2 = _rnl) ^ 275 ^^7°^^'^ , (5.42) 



K 



is used with the associated kaon mass being obtained by solving A^^(P^) = 0, where the inverse 
propagator is as given in Eq. ([4.25|) . The normalization constant is obtained from fj^ = '{~'^'k) 



which also ensures the correct normalisation of the charged kaon form factor. 



The obtained set of quark propagator parameters are those shown earlier in Table p.l| . The fit selects 
differences between the u and s quark propagators. The calculated electromagnetic form factor is 
reproduced in Fig. |5.5|. The comparison between calculated and experimental pion and kaon quantities 
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from tha work is presented in Table For the neutral kaon r'^o < 0, and Fko{Q^) is found to be 



similar in form and magnitude to the charge form factor of the neutron. The other quantities shown 



in Table |3.2| are obtained consistently from the same quark propagator model. 



The difference between the calculated and measured values of the charge radii and scattering lengths 
in Table \i.'2\ is a measure of the importance of final-state, pseudoscalar rescattering interactions and 
photon-vector-meson mixing, which are not included in generalised-impulse approximation. The 
above calculation suggests that such effects contribute less than ~ 15% and become unimportant for 
> 1 GeV^. The fact that the calculated values of fx/ fir and tk/tt^ agree with the experimental 
values of these ratios suggests that such effects are no more important for the kaon than for the pion. 

The results for F{Q'^) found by Burden, Roberts and Thomson (1996) for K+ and at < 2 GeV^ 
are not sensitive to details of the propagator parametrization. Their qualitative features are similar 
to those from a recent quark loop investigation based on a separable interaction (Buck et al, 1995). 
However specific magnitudes are different. As is the case of the pion, the large behavior is a 
faster fall-off than 1/Q^ due to both non-perturbative bound state amplitudes playing a role in this 
exclusive scattering process. The behavior of -FV±(Q^) at > 2 ~ 3 GeV^ is influenced by details 
of the Ansatz for the kaon Bethe-Salpeter amplitude, Eq. ( |5.42| ), that are not presently constrained 



by data. This emphasizes that measurement of the electromagnetic form factors is a probe of the 
bound state structure of the meson and of non-perturbatively generated differences between the u and 
s quark propagators. 



6 Meson Interactions 

6.1 The 7r'^77 Form Factor 

The pion charge form factor for space-like momenta is one of the simplest but non-trivial testing 
grounds for applications of QCD to hadronic properties. A closely related quantity that has received 
less attention is the 7*7r° — » 7 transition form factor (Ametller, Bijnens, Bramon and Cornet, 1992; 
Jaus, 1991). Here the photon momentum dependence maps out a particular off-shell extension of the 
axial anomaly (Adler, 1969; Bell and Jackiw, 1969). Presently available data for this transition form 
factor in the space-like region < 2.5 GeV^ is from the CELLO collaboration (Behrend et al., 1991) 
at the PETRA storage ring where the process e~^e~ e'^e^n^ was measured with geometry requiring 
one of the two intermediate photons to be almost real. There is renewed interest in this transition 
form factor due to the prospect of obtaining higher precision data over a broader momentum range via 
virtual Compton scattering from a proton target at CEBAF (Afanasev, 1994; Afanasev, Gomez and 
Nanda, 1994). The theoretical results (Frank et al, 1995) we summarize here can be viewed either as 
an application of the GCM model at the 0^^ level of the meson loop expansion or as an application of 
the DSE approach truncated to generalized impulse approximation. The improvements over previous 
quark loop studies (Ito et al, 1992; Anikin et al, 1994) are: quark confinement (thus eliminating 
spurious quark production thresholds), dressing of the photon-quark vertex, and the dynamical rela- 
tion between the pion Bethe-Salpeter amplitude and the quark propagator. The latter elements are 
crucial for producing the correct mass-shell axial anomaly independently of model details. 



Using the action given in Eq. ( p.21| ), the relevant interaction is generated by the first term iSi[7r; A]. 
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Figure 6.1: The quark triangle diagram for the generahzed impulse ap- 
proximation to the 7*7r°7 vertex. 

Expansion of the quark loop term shown as the first part of Eq. ( ^.22|) yields 

5[7r°77] = -Tt[S T^A, S V,A, S V^tx\ (6.1) 

That is 

r Pd^O 

5[vr°77] = j -l^M-P - Q)A.{Q)^'{P)A,.{P, Q), (6.2) 
where the vertex function is given by the integral 

A..(P, Q) = / Sik-P- Q)T,ik - I - |; -P - Q)S{k) 

xT^ik - |; Q)Sik - Q)Kik -^-Q-P). (6.3) 

The momentum assignments are shown in the quark triangle diagram of Fig. |6.1| . The general 
form of the vertex allowed by CPT symmetry is K^y{P, Q) = —i:^(-^uai3PaQi3 g{Q^, P^, P ■ Q) where 
^4123 = 1, a is the fine-structure constant, is the pion decay constant, and g is the off-mass-shell 
invariant amplitude. With the one photon mass-shell condition {P + QY = 0, the invariant amplitude, 
denoted by g{Q^, P^), is the object of the calculation. For a physical pion the shape of the 7*7r° 7 
transition form factor is given by glQ"^, —m^). The chiral limit for the physical 7r° 77 decay 
amplitude is fixed at by the axial anomaly (Itzykson and Zuber, 1980) which gives an excellent 
account of the 7.7 eV width and requires g{0,0) = 1/2. This follows only from gauge invariance and 
chiral symmetry in quantum field theory and provides a stringent check upon model calculations. 

In the work of Frank et al. (1995) the dressed photon-quark vertex Fj, is taken to be the Ball-Chiu 
Ansatz T^J' from Eq. ( |5.27| ) and the pion BS amplitude is taken to be the 75 part B{p'^ ,m) / f^^ as 
suggested by Eq. (|4.9| ). The 7r°77 vertex function in Eq. ( |6.3|) is now completely specified in terms 
of the quark propagator. An exactly parallel situation holds for the spacelike pion charge form factor 
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produced from the same approach (Roberts, 1996) and discussed in Sec. pT^ . The parameters employed 
there are used without adjustment for 7177. At = 0, numerical evaluation (Frank et ai, 1995) yields 
g^Oyy = g{0, —m\) = 0.496, in agreement with the previous application of this model (Roberts, 1996) 
and in good agreement with the experimental value 0.504 ± 0.019. The chiral limit of this approach 
has been shown (Roberts, 1996) to correctly incorporate the exact result g{0,0) = 1/2 produced by 
the axial anomaly independent of the form and details of the quark propagator. This is only possible 
when all elements of the loop are dressed consistently as dictated by symmetries. Here this means 
the chiral limit 75 term of the pion BS amplitude is fixed by chiral symmetry through Eq. ( |4.9| ) in 
terms of the quark scalar self-energy amplitude, and the physical soft photon-quark vertex is fixed by 
EM gauge invariance through the Ward identity in terms of the quark propagator. It is advantageous 
to account for the distributed nature of pions and dressed quarks, not only because it is a realistic 
outcome of QCD, but also because important symmetries are implemented this way. 
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Figure 6.2: The 7*7r'^7 transition form factor from Frank et al. (1995) 
along with the CELLO data from Behrend et al. (1991). 



The obtained form factor F^Q"^) = gi^Q"^, —ml)/g{0, —ml) at the pion mass-shell is displayed as 
Q'^F{Q'^) in Fig. |6]^ by the solid line. A "radius" or interaction size, defined via r^^o^ = —6 -^'(Q^) |q2=o) 
yields 0.47 fm while a monopole fit to the data yields (Behrend et al., 1991) 0.65 ± 0.03 fm. The insert 
compares the low result with a recent result from a QCD sum rule approach (Radyuskin, 1994), 
and a monopole form (Brodsky and Lepage, 1981) that interpolates from the leading asymptotic 
behavior F^Q"^) —>■ Sn'^f^/Q^ argued from the pQCD factorization approach (Lepage and Brodsky, 
1980). In the latter two approaches there is ambiguity due to: A) the unknown momentum scale at 
which perturbative behavior should set in; and B) assumptions for the pion wavefunction and how it 
should evolve with the momentum scale (Radyuskin, 1994). Within the present GCM/DSE approach, 
both the photon coupling and the pion wavefunction evolve with in a way determined by the evo- 
lution of the dressed quark propagators. This produces, in a single expression, both the ultra-violet 
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behavior required by pQCD and the infra-red hmit dictated by the axial anomaly. We note from the 
parameterization of the quark propagator given in Eq. ( p.23|) and Eq. (|3.24|) that the employed pion 
Bethe-Salpeter amplitude B{p'^,m)/ f.,^ has the correct leading power law behavior mX^/p^fT^ which 
implements the hard gluon contribution that dominates pQCD. 



The dotted straight line in Fig. B]^ is the pQCD factorization (Lepage and Brodsky, 1980) limit 
Q2f(Q2) ^ SvrV^ = 0.67 GeV^ Although the GCM/DSE non-factorized calculation reaches this 
value near 

Q2 = 3 GeV^ there is a slow decrease with higher consistent with a logarithmic 
correction. An excellent fit to the numerical results for 3.3 GeV^ < < 10 GeV^ is provided by 
F(Q2) = A [1.0+ B Q2 in(C Q'^)]-\ where A = 1.021, B = 0.461 /m^ = 0.777 GeV'^ and C = 
l.lQ/m?p = 1.45 GeV~^. (Neglected in this calculation is the anomalous dimension of the quark 
propagator; this would modify the power of the In-correction but would be a numerically small effect.) 

The logarithmic correction to the anticipated 1/Q^ asymptotic behavior can be attributed to the 
persistent nonperturbative nature of the coupling to the final state soft photon in this exclusive 
process (Frank et ai, 1995). Numerically it is found that, if a bare coupling were to be used for 
both photons as is implicit in the pQCD factorization approach, F{Q^) would eventually approach 
Sn'^f^/Q'^. The turn-over in Q'^F{Q'^) near 3 GeV^ predicted in Fig. ^]2]is barely within the Q'^ limit 
of 4 GeV^ anticipated for measurements at CEBAF if a 6 GeV electron beam becomes available. This 
turn-over and the logarithmic corrections generated by the loop integral are features also found in the 
parallel approach to the pion charge form factor (Roberts, 1996). For a study of the behavior of the 
vertex function off the pion mass shell, as needed for the anticipated CEBAF experiment, see Frank 
et al. (1995). 

This 7r77 approach has been recently extended without modification to the 77r* — > tttt form factor 
(Alkofer and Roberts, 1996). The value F^'^(4m^) of the amplitude at s = 4m^ from that work is 
shown in Table E?^. This process is also due to the chiral anomaly and it is found, again, that the exact 



chiral limit current algebra result at the soft point F^'^(O) = eNd {12tx'^ f^) is produced independently 
of the details of consistently dressed elements of the loop. Together with the parallel result (Praschifka, 
Roberts and Cahill, 1987a) for the anomalous five-pion Wess-Zumino term, Eq. ( |4.35| ), these exact 
results from symmetries are produced at a truncated level of the theory. The correction terms, 
including meson loops, can apparently be absorbed into the details of the quark propagator dressing 
to which the chiral anomalies are evidently blind. 



6.2 The piTTi Form Factor 



From the second term of the action in Eq. ( p.21| ), we identify S [pirn] = — Tr S V^- {S T^, ■ vf)" 
which yields 

where the vertex is 

A.(^, Q) = j tSi^pI^ + f ; Q)E.{k + |; -P - ^)E^{k - f ; P - |)T,(fc, P, Q) (6.5) 



(27r) 

T^{k,P,Q) = 2N,tis 



with 



(6.6) 
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These expressions have been simphfied through use of the approximations: only the i'^^TET^lq; P) part 
of the pion BS amphtude T.,, has been used, and the p BS amphtude has been represented by 
i^l^TVplq] Q) where 7J is transverse to the p momentum Q. The pnir vertex could have been obtained 
directly from the 77r7r vertex discussed in Sec. ^.21 . Since only the isovector part (T3/2) of Q from the 
photon-quark vertex contributes there, it is clear that the result for p^Trvr can be obtained by replacing 
for the photon with i'yJ^Vp{q; Q) for the p°. 

From symmetry properties, it is not difficult to show that Ap{P,Q) = — A^(— P, Q) = A^(P, — Q), 
which requires the general form 

A,{P, Q) = -P, Fp^^{P\ Q\ {P ■ Qf) -Q,P-Q H^^P^ Q\ {P ■ Qf). (6.7) 

With both pions on the mass-shell, (P — y)^ = (P + y)^ = —m^. Equivalently, P ■ Q = and 
P^ = — — ^ so that only the first term of Eq. ( |6.7| ) survives. The corresponding form factor 
FpTTTriQ"^) contains the coupling constant as its mass-shell value, i.e. Fp^^T^^Q"^ = —m?^) = gp-,,-,,. If the 
form factor is held at this triple mass-shell value independent of momenta, then Eq. ( |6.4|) , when 
converted to coordinate representation, becomes 

Spt [ptttt] = -gp^^ j d^x pp{x) ■ ti{x) x dpTi{x) , (6.8) 

which is the point coupling limit expressed in a standard form. There are, of course, other forms 
of action for this interaction that are on-shell equivalent. The general form of Eq. (|6.4| ) is just one 
of them (in coordinate space it is a complicated non-local form). The momentum space form of 
Eq. (|6.4| ) is an efficient way to account for the non-locality or off-shell structure of the interaction 
vertex produced by the quark-gluon substructure of the tree-level quantities within the GCM. One 
should bear in mind that the ofF-mass-shell behavior here in the meson momenta depends on what 
has been identified as the propagator for the qq meson modes as discussed in Sec. pT^ . It is always 
subject to a field redefinition and one must use a consistent definition throughout all elements used 
to construct an observable or S-matrix element. 



Previous studies of the pvrTr vertex within the GCM (Praschifka et ai, 1987b; Roberts et ai, 1988b; 
HoUenburg et ai, 1992) made the approximation that gpT^^, could be estimated at the point P = Q = 0. 
This was necessary to avoid ambiguities associated with numerical extrapolations of dressed quark 
propagators and vertex amplitudes in the integral Eq. ( |6.5| ). The parameterized analytic representation 
of the quark propagator described in Sec. ^ allows the correct mass-shell constraints to be maintained 
for (yfpTTTT and like quantities. Since the propagator parameters have so far been constrained only by 
pion physics, there is no guarantee that continuation of S{p) away from the real axis by an amount 
proportional to will be correct. The exploratory calculations involving vector mesons described 
here constitute a first investigation of this matter. 

The piTTT vertex has recently been studied within the GCM/DSE approach by Tandy (1996) and also 
by Mitchell and Tandy (1997). Those works use £'^(g;P) = P(g^,m)//^ and the parameterization 
r^(p^) oc . The strength is set by the canonical normalization condition given in Eq. ( p.61|) . 



The range a is then adjusted to reproduce the empirical value g^"^^ = 6.05. This approach has proved 
to be phenomenologically successful for studies of other processes such as p — mixing (Mitchell et 
ai, 1994) and the p — > e"*" e~ decay (Pichowshy and Lee, 1996) from essentially the same theoretical 
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Figure 6.3: The pvrvr form factor versus p momentum from Mitchell and 
Tandy, (1997). 



framework. The produced p ^ tttt decay width, given by 

2 

-p 9 p-KTV "^p 

is 151 MeV. The calculated form factor Fpj^T^{Q'^) is shown in Fig. |6.3| for timelike and spacelike 
momenta in the vicinity of the mass-shell at —0.6 GeV^ indicated by a dot. The main conclusion from 
this calculation is that the previous approximation of using zero momentum to extract a coupling 
constant (Praschifka et ai, 1987b) can underestimate the value by almost a factor of 2. 

These investigations are exploratory in the sense that an independent calculation of the p BS amplitude 
Tp{p) has not been employed. There are indications that an improved parameterization of V^(p^) to 
include a more realistic power-law fall-off at high momenta provides a simultaneous accounting for 
both p TTTT and p — e"*" e~ decays (Roberts, 1997). Covariants other than the canonical one 7J for 
the p BS amplitude can contribute to the pnn vertex and this has recently been investigated together 
with the effect of the sub-dominant pion covariants (Qian and Tandy, 1997). Truncation to the 
dominant p amplitude is found to be an accurate approximation while the effect of the pseudoscalar 
pion component can be quite important. 

The strong suppression displayed in Fig. p.3| as increases from the mass-shell point to the spacelike 
region indicates that the effective coupling strength gp„^(Q^) appropriate to the p contribution to the 
pion charge form factor and radius is significantly smaller than what is assumed in the standard VMD 



1 - 



4m! 



0/ ^ 



(6.9) 
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approach. This is consistent with findings discussed in Sec. |^. The strong dependence in Fig. |^ 
is very similar to that found recently (Frank and Roberts, 1996) in a study that used a p amplitude 
from solution of the BSE. In that work the time-like region of the form factor could not be calculated 
directly because the numerical solution of the p BSE was performed only on the real > axis. The 
results there suggest that the common procedure (e.g. see Eq. (56) of Frank and Roberts (1996)) of 
interpolation of calculated form factors in the meson momentum from the Euclidean domain to the 
mass-shell point can be reliable if the mass is not too large. 



6.3 The 77r/9 Form Factor 

The empirical p BS amplitude set by (7^,^77 as above enables a parameter-free prediction for the 77rp 
vertex. Apart from being a consistency check in this manner, the 7'7rp interaction together with the 
77r7r, 777r, and the pTrvr processes, provide important guidance for extending the present approach for 
nonperturbative QCD modeling of meson physics beyond phenomena dictated by chiral symmetry. 
Within nuclear physics, the associated isoscalar 7*7rp meson-exchange current contributes significantly 
to electron scattering from light nuclei at presently available momenta. In particular, our understand- 
ing of the deuteron EM structure functions for ^ 2 — 4 GeV^ is presently hindered by uncertainties 
in the behavior of this form factor (Ito and Gross, 1993). 

.6 I \ \ \ \ \ \ I \ \ \ \ \ \ \ \ \ \ \ \ 




Expansion of the second term of the action in Eq. (|5.21|) to first order in the EM field yields the 77rp 
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interaction as the pair of contributions 



5[77rp] = -Ti[S T^A^S ■ tt S f^ p.] - Tt[S T^A^ S f^ ^ ■ vf]. (6.10) 
With a vertex function A^j^(P, Q) defined by 

5[7^p] = - / A^Q) 7f(-P - |) ■ p,{P - |) A,.(P, Q), (6.11) 

one may combine the two terms in Eq. ( |6.10| ) to obtain the integral 



X 2iV,tr,[5(A:+ - |)r.(fc+; Q)5(A;+ + |)z755(fc_)z7j] • (6.12) 
Here k± = k ^ ^. The approximations made here for the vr and p BS amphtudes are identical to 



those used for the pTrvr interaction in Sec. |6.2| . Use of the Ward-Takahashi identity from Eq. (|5.8|) in 
Eq. ( |6.12| ) shows that in this formulation the 77rp EM current is conserved, that is K^y{P, Q) = 0. 
The general form of the vertex function can be shown by symmetries to be 

A^,(P, Q) = -I— e^,^p P^Qfs gp^^ f{Q\ P\ P ■ Q) , (6.13) 
nip 

as is expected for a coupling arising from the chiral anomaly. We have used the standard definition 
of the coupling constant so that, at the triple mass-shell point, the form factor / = 1. 

The coupling constant and form factor resulting from the GCM or DSE approach truncated to the 
level shown above have been obtained by Mitchell (1995) and by Tandy (1996). The form of the 
quark propagator parameterization is that discussed in Sec. || and the Ball-Chiu Ansatz is used for 
Fj,. Those works obtain (yf^^p = 0.5 The experimental p^ ir'^'-f partial width (67±7 keV) determines 
the empirical value g^^p = 0.54 ±0.03. The 77rp form factor obtained with on-mass-shell vr and p, and 
weighted by Q^, is shown by the solid curve in Figure Also shown is the form factor suggested 
by vector meson dominance (VDM) as initially used to treat such meson exchange effects in electron 
scattering analysis (Hummel and Tjon, 1990), as well as the result from a free constituent quark 
loop with bare photon coupling (Ito and Gross, 1993; Ito et ai, 1992). Both quark-based approaches 
produce a much softer form factor than does the VDM assumption. Above 50 fm~^, which is readily 
accessible in electron scattering, the differences are serious. Since this calculation for 77rp is equivalent 
to the O*'^ term in the meson loop expansion. An estimate of meson-loop corrections to this result 
would be interesting. 



7 Beyond Meson Tree Level 

The essence of the GCM approach is that the tree-level meson fields are ladder qq bound states whose 
structure and interactions are governed by what we might call their dressed quark core. Further 
quantum treatment will generate meson loop effects from integration over these fields. That is the 
dressing of ladder mesons by ladder mesons. This extends the bound state structure of the mesons to 
add meson cloud effects to the original dressed quark core. In the more general DSE approach, such 
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processes are implemented by use of less severe truncations of the equations of motion; also exposed 
that way in the DSE approach will be contributions from explicit gluon n-point functions that are 
not already accounted for by dressed quarks and mesons. In either case, it is clear that the approach 
just outlined is out of practical reach if one has to proceed too far into the loop expansion or with 
higher ra-point functions to capture the dominant low energy physics. A crucial question then is the 
size of meson loop contributions to the various processes of interest, such as those observables we have 
described in previous Sections. In the GCM/DSE approach, the composite and extended structure 
of the tree-level meson fields makes such investigations difficult. On the other hand, the distributed 
vertex functions for the coupling tend to suppress the numerical importance of the meson dressing 
loop integrals, and it might be spectulated that for this approach the meson loops play much less of 
a role than for point coupling models or for effective field theories built from point coupling of point 
interpolating meson fields. We shall review the few investigations that have tackled this question in 
the GCM/DSE approach. Besides the pion loop contribution to the pion charge form factor, we also 
include here a discussion of vector meson mode contributions to this quantity. 



7.1 Role of Vector Mesons in the Pion Form Factor 

The relative unimportance of the p for the space-like pion charge form factor in the present quark- 
gluon approach is consistent with considerations of vector-meson-photon mixing in QCD. It can be 
understood in terms of Eq. ( ^.16| ) which is the dressed photon-quark vertex Vy separated into a p pole 



term and a non-resonant or direct background term. With both terms of that vertex employed for 
coupling to the pion in the manner discussed in Sec. |5.2|, one obtains 



^ + ^^gl^igp . (7,1) 

The first term is from the non-resonant photon vertex term and is the generalized impulse approxi- 
mation (GIA) that we have described in Sec. f).2\ The second term is the p resonant term that arises 
in the GCM or DSE approaches. The relative importance of the two terms away from the p mass- 
shell depends strongly upon their momentum dependence. The off-mass-shell momentum behavior of 
individual terms in photon vertex separations such as T^, = F"'" + F^"'^, given in Eq. ( p.l6| ), and that 
underly Eq. (|7. 1|) , will differ in different approaches. The total result should not. 

In the naive limit of a point coupling model with structureless hadrons, the elements of Eq. (|7. 1|) 
become F^^'^^Q'^) 1, Fp^.^lQ'^) Qp^,^, and nj7(Q^) —Q^/gv- The latter Q"^ factor arises from 
the description of the mixing as proportional to / d'^xpp^i,{x) F^i,{x). This produces the VMD empirical 
form 

fM) = 1 - • (7.2) 

where the first (contact) term accounts for the pion charge and the charge radius is required to 
come totally from the second (p propagator) term. This implies ~ ^gpmr/i'fn'jtgv), which under 
the universality assumption of vector coupling gy ^ gpT„v, produces ~ 0.4 fm^ which compares very 
well with the experimental value 0.44 fm . In this sense it is often considered that the spacelike pion 
charge form factor is the tail of the p resonance. A recent re-examination of the VMD phenomenology 
for this and related processes with two choices of point hadron lagrangian can be found in the work 
of O'Connell et ai, (1995). We note that the phenomenological form (VMDl) found preferable there 
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obeys the massless photon constraint n^''^(0) = and produces just the point hmit in Eq. (|7^). We also 
note that the more standard form of VMD, called VMD2 by O'Connell et ai, (1995), uses universal 
vector coupling and also the assumption that Ufp^Q"^) ~ Tn^p/gv for all Q^] then the equivalence of 
Eq. (|7.2| ) to m?p/ (m^ + Q^) is often used to invoke the extreme form of VMD, namely that the photon 
couples only through vector mesons. 

From the perspective of the GCM or DSE approaches, the reality of a space-time distributed quark- 
gluon substructure of the pions, for example, demands the distributed form for the form factor shown 
in Eq. ( |7.1| ) rather than the point form in Eq. ( |7.2| ). These considerations certainly apply also to QCD. 
The question of the relative importance of the two terms in Eq. (|7.1| ) then arises. This will receive 
different answers in different formulations for it depends on how the separation is defined off-mass-shell. 
This corresponds to the observation that there is no unique way to separate a pole contribution from 
a background contribution except precisely at the pole. The sum of the two terms is the meaningfuU 
quantity. With adoption of the Ball-Chiu Ansatz, this representation of the non-resonant term of 
the photon-quark vertex shares with the BSE vertex those aspects of momentum dependence that 
are constrained by symmetries including the Ward-Takahashi identity. Since the vector mesons are 
described by the same kernel, it should not be surprizing that at least some extrapolation of vector 
qq correlations into the spacelike region is accounted for through the Ball-Chiu vertex. A similarity 
in strength between the Ball-Chiu vertex and the spacelike extension of low mass vector meson poles 
has recently been found (Frank, 1995) through a numerical study of the ladder BSE for the photon 
vertex. It is quite a robust result that about 85 — 90% of the experimental value of is produced 
by the resulting non-resonant or generalized impulse term in the work of Roberts (1996); the value 
of /^r^r is quite insensitive to variation in the representation of the dressed quark propagator if other 
constaints such as m.,^ and {qq) are maintained. 

This leads to considerations of why the p contribution should be small in such an approach. There is 
a constraint at = 0. Electromagnetic gauge invariance requires the p-7 mixing amplitude Ilfp{Q'^), 
and hence the pole term of -F^((5^) in Eq. ( [7. 1| ) , to vanish at = 0. Asymptotic freedom ensures 
Utp^Q'^) vanishes at large spacelike-Q^; hence a small contribution to F^^IQ"^) from the 7p mixing 
term over a large range of space-like can be consistent. A more detailed investigation requires a 
consistent definition of the momentum dependence of the separation in Eq. ( ^.16|) as related to the 



p BS amplitude and the defined non-resonant photon vertex. An initial investigation has utilized a 
mass-shell BS amplitude for the p-meson which is well described as a qq bound state within the ladder 
BS approach; such a dressed-quark core gives ~ 90% of its mass (HoUenburg, Roberts and McKellar, 
1992; Burden, Qian, Roberts, Tandy and Thomson, 1996). A calculation (Roberts, 1996) based on 
this confirms that Il'^^{Q'^) is a slowly varying function for spacelike Q^. The momentum suppression 
of the pole term of Eq. ( [TID due to -Fp7r7r(<5^) is quite strong as discussed in Sec. and displayed in 
Fig. |673| . We may use the representations FpT.T.^Q'^) = gp-^Trfp-K-niQ'^) and Ili^ [Q"^) = —Q^fp'y{Q'^)/gv, 
where the two new functions / and / depend on the details of the dynamics. The p pole contribution 
to the pion charge radius from Eq. (|7.1| ) can then be written in the form 



{r'/r = rl- (rf ^)^ = 1.2 /,..(0)/,,(0) , (7.3) 



p 



where we have used the empirical result gpTnr/gv ~ 1-2. Since from Fig. |O| /o^^(0) ~ 0.5, one estimates 
that a value /p7(0) ~ 0.4 would account for all of the difference between the non-resonant (r^^"^)^ in 
Table |3]^ (0.31 fm^) and the experimental value (0.44 fm^). A reasonable value fp^{0) ~ 0.3 would 
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leave room for a contribution as well from the pion loop process of the size estimated by Alkofer et 
a/. (1995) and described below. 



7.2 Pion Loop Shift of the p Mass 

At the tree level of the GCM, the ladder description of the p and u is degenerate, nip = m^. The 
pion loop correction to the p self-energy produces a mass shift and also a width due to the open 27r 
channel. Within the GCM, this mass shift was originally studied by Hollenburg et a/. (1992). We 
review the later work on the shift and width by Mitchell and Tandy (1997) and highlight some issues 
that generally arise in the GCM or DSE approaches. 

Mitchell and Tandy (1997) treat the transverse p and u modes and use fields that contain the transverse 
projector T^,(P) = 6^, - P^PJ P\ That is, p^(P) = T^,{P)p,{P) and u^iP) = Tp,{P)uj,{P). In a 
matrix notation where denotes (p^,c<j^), the tree-level action from Eq. ( p.63|) , up to second order 
in the fields, can be written as 

Here, in meson channel space, is diagonal and the only non-zero elements of H''' are the off-diagonal 
ones that provide p^ — uj coupling via the u — d current mass difference. The superscript q identifies 
the quark-loop contribution to distinguish from subsequent pion loop quantities. The composite qq 
nature of the mesons is reflected in the fact that the diagonal inverse propagators contain a dynamical 
mass function. The general form 

A-i(p2) = {P^ + ml)Z^\P^) = (P^ + mliP^)) (7.5) 



as calculated from the formalism in Sec. p.2| identifies the degenerate mass my of the tree-level p 
and u states. The momentum dependence of Zi(P^) is due to the qq substructure. Without p° — u; 

mixing, a physical normalization (unit residue at the mass-shell pole) would be produced by absorbing 

_ 1 _ 1 

at least the on-mass-shell value of Zi ^ into the fields. The function Zi ^ (P^) may be absorbed into 
the fields producing unmixed propagators have the standard point meson form. As in Sec. the 

associated BS amplitudes become correspondingly rescaled as ry(g;P) = ry(g;P) Zf^(P^), so that 
at the mass-shell the standard normalization condition (Itzykson and Zuber, 1980) is satisfied. 



The relevant low-order terms in the effective action for vr, p, uo from the bosonization result in Eq. (p.64|) 



are 

S [vr, p, uj] = §2 [p, uj] + §2 [it] + Slpmr] + S[Lunn] H . (7.6) 

Here the quadratic terms 5*2 [p,uj] and 5*2 [tt] represent respectively the p — u sector containing the 
quark-loop mixing mechanism and the free pion sector. To integrate out the pion fields, it is convenient 
to first combine the last three terms of Eq. ( [7.6|) so that it becomes 



S[n,p,uj] = S[p,uj] + ^ I -j^ix.,{P')D-^{P\P)'n,{P). (7.7) 



1 fd^P,P' 
The three terms in 

DrMp\ P) = ArMp\ P) + V,,{P', P) + W,,{P' , P) (7.8) 
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correspond to the tree-level pion inverse propagator given by 

Ar;(P', P) = {2^)%^6\P' + P)K\P\ (7.9) 
the pvTTr interaction term given by 

V,,{P\P) = 2ie,,,pl{-P' - P) A^(-^^^; -P' - P), (7.10) 
and the utttt interaction term given by 

W.,{P\P) = 2ie,,,u,{-P' - P) A;{-^^^; -P' - P). (7.11) 

In this work the pion inverse propagator of the GCM is reduced to the standard point form by 

absorbing the calculated Ztt ^ (P^) into the pion field. Thus the pion momentum dependence of the 
pTTTT and tUTTTT vcrtex functions contain the compensating aspects of the qq substructure of the pion. 
The pTTTT vertex function is described in Sec. |6.2| . The corresponding cutttt vertex function, which is 
generated by a u — d current quark mass difference, is calculated and included by Mitchell and Tandy 
(1997) to extract the pion loop contribution to p — uj mixing. Integration over the pion fields allows 
a new effective action iS[p, cj] to be identified from 

Z = N j DpDujD7iexp(^-S[7r,p,uj]) = N' J DpDco exp (^-S[p,uj]) (7.12) 

by using the functional integral result 

/ ^vrexp (^-i J ^j^MP')DT^\P\P)n,{P)^ = exp [-^TrLn D~'^ . (7.13) 

The field-independent term exp(— iTrLn A~^) from the result may be absorbed into the normalization 
constant to arrive at 

S [p, iu] = S2 [p, uj] + -TrLn(l + A{V + W)). (7.14) 
The second term here gives the vector meson coupling to all orders produced by a single pion loop. 

The quadratic part of the second term of Eq. ( |7.14| ) adds a contribution to the tree-level p° — matrix 
inverse propagator. The net transverse result is 

j.-Up2. _( P' + ml + W^{P^) W-{P^) \ 

^ ^^>-\ n^-(p2) p2 + m2, + n--(p2) j' ^^■^^> 

where the net mixing amplitude is n^'^(P^) = n''(P^) + n'^(P^). The additional terms produced by 
the pion loop are the diagonal self-energies 11^'^ (P^) and 11'^'^ (P^), and the mixing term n'^(P^). From 
Eq. (|7.14|) , the self-energy contribution I\.pp{P'^) is given by the transverse component of 

^7^{P) = -4 / ^A.(g+)A.(g_)A^(g; -P)K{t. P), (7.16) 

where q± = q ± P/2, and the summation over pion isospin labels has been carried out. This contribu- 
tion is the same for each isospin component of the p. The u self-energy contribution is obtained from 



76 



0.10 

0.05 
-0.00 
-0.05 

-0.10 

> 

& -0.15 
^ -0.20 

a 

^ -0.25 
-0.30 
-0.35 
-0.40 

-1.0 -0.5 0.0 0.5 1.0 

Figure 7.1: The real and imaginary parts of the pion loop contribution 
to the p self-energy. Taken from Mitchell and Tandy (1997). 

Eq. ( [7.16[ ) by the replacement of by A"^. It is quadratic in the small symmetry breaking mechanism 
and is ignored. 

In the absence of mixing, the real part of 11'''' (P^) generates a mass shift for the isospin eigenstate pj. 
This mass shift now represents the p — uj mass splitting. For timelike momenta such that < — 4m^, 
both IIPP{P'^) and n'^(P^) have imaginary parts associated with the decay p — > 27r. This mechanism 
should produce most of the p width. The mass modifications due to mixing are second order in a very 
small quantity and are ignored. The shifted rup is therefore identified from the position of the zero in 
the real part of the relevant inverse propagator in Eq. ( [7.15| ). That is, 

ml = ml + ReUPP{-ml), (7.17) 

where now = my. With the definition f (P^) = — Im 11'''' (P^)/mp, the inverse propagator in the p 
channel can be written as 

V;\P') = (P2 + - tmpTiP')) Z,\P'), (7.18) 

where r(P^) = f (P^) Z2{P'^). The function Z2{P'^) arises from the momentum dependence of the 
real part of the pion loop self-energy. The on-mass-shell value, which is the p field renormalization 
constant due to the 27r content, is given by 

= 1 + Ren'""(-mJ), (7.19) 

with the prime superscript denoting differentiation with respect to the argument. One absorbs 
_ i_ 

Z2 ^ (P^) into the field p^{P) so that the resulting propagator has the conventional Breit-Wigner 
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form. The physical width is given by 



= -m-^ Z^i-ml) lmnPP{-ml). (7.20) 
Numerical results for the p mass shift and width are obtained from Eq. ( 7.16| ) in the specific form 
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Figure 7.2: Graphical determination of the p — uj mass difference due to 
the pion loop. The influence of the qq composite nature of the vr and p 
propagators is also shown. Taken from Mitchell and Tandy (1997). 



(Mitchell and Tandy, 1997) 



n^^(P^) = --f%AP' 



(g2 - (g . Pf/P^)e 



-2g2/A2(p2) 



(27r) 



(g + f )2 + ml 



le 



{q - f )^ + ml- ie 



(7.21) 



where a parameterized form of the calculated pmr vertex function is used for convenience. The nu- 
merical results for IieIl'^'^{P^) and Imn''''(P^) are shown in Fig. [7.1| . The self-consistent solution for 



produced by Eq. (|7.17| ) is displayed in Fig. ^]2| in the following way. The quantity 



IS 



plotted as a long dash line and its intercept with the solid line representing Ren^''(P^) identifies the 
mass shell point. With = my = 782 MeV one finds rrip = 761 MeV giving rrip — rriuj = — 21 MeV. 
The experimental value is —12.0 ± 0.8 MeV. The calculated width is Tp = 156 MeV while the experi- 
mental value is 151 MeV. Fig. also illustrates the contribution to the mass splitting result made 



by the momentum dependence of the meson self-energies generated by the qq substructure. This is 
most marked for the pion and enters the pion loop integral in Eq. ( |7.16| ) through the factors -Z^(g^) 
which are contributed by each pion propagator. (In the calculational procedure these factors have 
been moved into the vertex functions for convenience.) If each function Z^^ is held constant at its 
mass-shell value (as would be the case for a structureless pion), the full calculation in Fig. |7]^ (solid 
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line) becomes the dot-dashed hne. The intercept then indicates that the mass shift would be essen- 
tially zero. In contrast to this, the short dashed line indicates that the influence of the corresponding 
quantity Zi{P^), from the qq substructure of the p, is quite negligible here. This is because the range 
of variation in is very small. 




Figure 7.3: The pion charge form factor. The first term is the general- 
ized impulse approximation, and the second term is the one pion loop 
correction. Double lines represent the ladder qq pions of the GCM/DSE 
approach. Quark lines are dressed propagators and circles are dressed 



7.3 Pion Loop Part of the Pion Charge Radius 

A somewhat more complicated task is presented by the pion loop contribution to the pion charge 
radius. This is illustrated in Fig. [7.3| . The analysis has been carried out within the GCM framework 
by Alkofer, Bender and Roberts (1995) in the following way. Starting from the tree level GCM meson 
action of Eq. ( ^.25 ), which includes coupling to a background EM field, the first-order contributions in 



the pion loop expansion may be developed for various vertex functions. For the 'yinr vertex function, 
Alkofer et al, (1995) obtain 

A,{P; Q) = Kf\P- Q) + AJ'°^(P; Q) , (7.22) 

where the first term is the generalized impulse approximation result (i.e., at meson tree level) given 
by Eq. ( ^.31| ) and discussed in Sec. |5]^. The pion loop contribution is given by 



AJ"'^(P; Q) = -\j 0^J'''\Q^ P-K,K + P) A.(ir+) A.(A-_) e'"^ Kf^{K- Q) , (7.23) 

where K± = K ± Q /2 and T^-^'^^ipi, P2] Ps, P4) is the n — tt scattering kernel which may be written 
in terms of a single invariant amplitude A{pi,p2;p3,P4). A threshold expansion for this amplitude in 
low powers of momentum is not adequate here since the integral in Eq. ( [7. 23] ) samples a large domain 
of spacelike momenta. The contribution to the required A{pi,p2;p3,P4) that yields both the leading 
order chiral limit result at low s, the term (Weinberg, 1966) given in Eq. ( |4.36| ), and the leading 
contribution within the GCM at tree level for large s, is obtained from the action in Eq. ( [4.33| ) in the 
form 

Sb\U] = J d^xd'y ^tr [d,U{x)d,Uiyy] f{x - y) , (7.24) 
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which avoids the derivative expansion by retaining the nonlocahty. Here f{x — y) is the fourier 
transform of the momentum space quark loop integral with two chiral insertions. It is of the same 



form that contributes to the pion inverse propagator as discussed in Sec. ^A. 



For the integral in Eq. (|7.23|) , the onshell approximation A'^^'^{K; Q) ^ 2K^F^^^{Q'^) is made. The 
result is expressed in the form 



F,{Q'') = F^^\Q'') l + L{Q 



(7.25) 



where L is a result from the pion loop integral. Since the condition ^,^(0) = 1 is equivalent to the proper 
normalization condition for the pion BS amplitude, the above form emphasises that the pion loop 
mechanism generates additional internal structure for the pion. A renormalization of the modified BS 
amplitude is necessary. The form z75_B(g^)//^ used at each nqq vertex is now effectively supplemented 
by an additional pion cloud term. The normalization constant previously determined at the quark 
loop level, must be renormalized to adjust for the pion cloud now being included in the pion. From 
these considerations, Alkofer et ai, (1995) report results in the form 



r. 



loop \ 2 



(7.26) 



where the first-order pion loop contribution 



loop I 



is at most 10 — 15% of the generalized impulse 



approximation (meson tree-level) result. The loop term has the expected divergence as a chiral 
logarithm in the m^r limit, but it has the stated small effect at the physical point and is seen to 
not dominate until m^r ~ 10 MeV. The pion loop correction to is found to be at most 2% of the 
impulse or tree- level result. 



In ChPT with Nf 
is 



3, the result given by Gasser and Leutwyler (1985) up to one pseudoscalar loop 



12 a 



fl 327r2 /2 



2 In 



+ In 



m 



K 



+ 3 



(7.27) 



where Lg is the 9*'' low energy constant at O(p^) 



after the divergent part from the loop has been 



absorbed and /z^ is the loop regularization scale. Any reasonable choice of /x^ such as m^, leads to the 
"chiral logarithms" in Eq. (|7.27|) contributing typically 15% of the experimental value 0.44 fm^. The 



dominant Lg term is the counterpart in ChPT of 



MIA\2 



obtained from the distributed quark-gluon 



substructure of the 77r7r vertex in the GCM/DSE approaches. However the correspondence is only 
approximate; the distinction between tree-level and loops in the two approaches is not the same due 
to the very different character of the pseudoscalar fields. This has been discussed in the early part of 
Sec. O 



8 Concluding Remarks 

We have summarized the Global Color Model as it applies to a variety of meson physics observables. 
A constant theme that emerges is that the tree level interactions among qq Goldstone mesons and 
photons are well-described in terms of only the dressed quark propagators. This efficiency is due 
to the connection between the pseudoscalar Bethe-Salpeter amplitudes and the quark propagator 
brought about by chiral symmetry and also by the connection between the photon-quark vertex and 
the propagator brought about by electromagnetic gauge invariance. Only tentative steps have been 
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taken so far in the task of relating the phenomenological aspects of the confining quark propagators 
used at this stage to an underlying gluon 2-point function. This step is required if the GCM is 
to fulfil the initially designed role of a field theory model based upon an effective gluon 2-point 
function. This need is emphasized by investigations of mesons other than the Goldstone bosons 
where more approximations must necessarily be made. However, even with an increased degree of 
phenomenological input, valuable information is being obtained. 

We have seen that various couphng constants obtained at meson tree-level are within 10% of exper- 
iment. At that stage the meson modes from the GCM have ladder Bethe-Salpeter content and the 
meson coupling is via a dressed quark loop. An interesting question is whether this will be a persistent 
theme that leaves meson loop corrections relatively unimportant in this approach. This is possible 
because the subsequent meson loop dressing would have to overcome distributed coupling caused by 
the finite size effects of the tree-level mesons in this approach. There are only several isolated pieces 
of evidence. Pion-loop corrections to the pion EM form factor were seen to contribute at the level 
of < 15% to the charge radius (Alkofer et ai, 1995). Similarly, the pion-loop contribution to the 
downward shift in the p mass relative to in the GCM has been found to be < 5% (HoUenburg 
et al., 1992) and in a more recent study 2% (Mitchell and Tandy, 1997); this latter result is in good 
agreement with experiment. It would certainly be a convenience if the low- mass mesons and their 
soft interactions were amenable to a representation in which their dressed quark content captures the 
dominant quantum loop effects. However this remains to be determined and in this regard future work 
within the GCM is called for to obtain the pion loop contributions to the tt — tt scattering lengths espe- 
cially Cq. In ChPT at 0{p^), the only distributed substructure enters through loops of the (auxihary) 
point pion fields, while at tree level in the GCM there is significant substructure through all orders in 
momenta. A systematic comparison through a complete set of threshold observables at 0{p'^) would 
be quite useful. One of the features of the GCM is that constraints from asymptotic freedom at large 
momenta are automatically included and the complete momentum dependence of form factors and 
scattering processes can be examined. Prom this perspective, one may speculate that with a model 
field theory with a few parameters such as the GCM, high accuracy in threshold behavior might be 
incompatible with quality descriptions of many processes over a large momentum domain. More work 
is needed to explore a large variety of observables. 

There has been significant effort in the past devoted to investigations of baryons (particularly the 
nucleon) within a covariant Paddeev formalism truncated to a quark-diquark approximation derived 
within the GCM. Probably the first covariant Paddeev calculation of the nucleon (Burden, Cahill and 
Praschifka, 1989) was carried out within the GCM framework. Other works (Cahill, Roberts and 
Praschifka, 1987; 1989a; Praschifka et ai, 1989) have made important contributions to understanding 
the properties and role of diquarks as constituents of baryons. Progress in this endeavour is compli- 
cated by the rich singularity structure of the covariant Paddeev equations. A more accessible, but 
necessarily more approximate approach to the nucleon has been via chiral meson-quark mean field 
models of the non-topological soliton variety. The GCM has provided a definitive underpinning to 
such methods by producing the soliton equations of motion for these degrees of freedom as dressed 
composite modes based at the quark level (Cahill and Roberts, 1985; Prank, Tandy and Pai, 1991). 
As a confining quark condensate model of the nucleon it has provided very efficient and sensible nu- 
merical results (Prank and Tandy, 1992). Space has not permitted an analysis of such baryon work to 
be included here. 
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